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CII. Gamma Radiation from the Decay of mee PUsa, * "CM, and *8Cm,, 


By J. O. Newton, B. Rosz and J. Minstep 
Atomic Energy Research Establishment, Harwell + 


[Received February 10, 1956] 


ABSTRACT 


Accurate energy measurements on the gamma rays from the decay of 
*38Pu and ?42Cm have been made with proportional counters. Values of 
43-49-+- 0-08, and 99-8+-0-4 and 153-1--0-6 kev for 238Pu, and 44:03-L0-06, 
101:80+-0-17 and 157-61-L0-3 kev for 242Cm were obtained. The ratios 
of these energies for each nuclide agree well with those predicted in the 
strong coupling limit of the Unified model. Gamma rays having energies 
of 210--1-5, 228+2 and 277+2 kev and K—X radiation were found from 
the decay of 248Cm. These gamma rays are shown to be magnetic dipole 
with a 30%, mixture of electric quadrupole ; the magnetic dipole com- 
ponent is 104 times slower than the single particle estimate. 


$1. INTRODUCTION 

Tuis work was originally undertaken with the aim of investigating the 
modes of decay of the alpha particle emitters 78*Pu and 24Cm. At 
that time (early 1953) there was little evidence regarding the second and 
higher excited states of the heavy even—even nuclei and their systematic 
variation with mass number was not realized. It was apparent from our 
early measurements that the gamma ray spectra from the decay of 
238Pu and 242Cm were similar. Later it was found that the ratios of the 
energies of these gamma rays fitted very well with the values to be expected 
if they arose from the decay of nuclear rotational states as predicted by 
Bohr and Mottelson (1953) from the unified model, and as demonstrated, 
with rather less precision, by the work of Asaro and Perlman (1954) and 
Asaro, Thompson and. Perlman (1953) on the alpha-particle fine structure 
of the even—even elements. 

In view of the rather precise predictions of the unified model regarding 
the relative spacings of energy levels it seemed desirable to make the 
energy measurements as accurate as possible. Krypton and xenon filled 
proportional counters of the type previously described by West (1953) 
were used. The main advantage of these counters is their resolution, 
approximately 5°%, full width at half height for gamma ray energies of 
100 kev. Their intrinsic efficiency for detection is also quite good, being 
about 1-5° for radiation of 100 kev in a 2 in. diameter counter filled to 
atmospheric pressure with xenon. 


+ Communicated by the Authors. 
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A preliminary account of the work on 238Py and 2420m was presented 
at the Birmingham Nuclear Physics conference in July 1953 (Rose and 
Newton 1953) and the later results have also been quoted in other publica- 
tions (Bohr 1954, Newton 1954). The energy values given in this paper 
differ slightly from those previously reported because the energy of the 
gamma ray from the decay of *4*Am, which was used as a standard for 
these measurements, has recently been revised. It was originally 
reported as 59-78-+0-02 kev (Browne 1952, Asaro, Reynolds and Perlman 
1952) but is now given as 59:568-+0-017 kev (Day 1955) and 59-62+ 
0:06 kev (Jaffe e¢ al. 1955). 


§2. PREPARATION OF 74Cm anp 73%Pu SOURCES 


Curium 242 was prepared by the neutron irradiation of americium 241 : 
Am 241 (n,y) Am 242m oe Cm 242. 

Two samples of curium were used in the present work, one from an 
irradiation in BEPO at Harwell, and the other from a much more intense 
irradiation in the NRX reactor at Chalk River. The integrated neutron 
doses for these irradiations were approximately 21018 and 10° 
neutrons/cm2, respectively. It is unlikely that an appreciable amount of 
Cm 243 was formed in the BEPO irradiation, but a fractional percentage 
by weight of this isotope could have been formed in the NRX irradiation, 
by further neutron capture in Cm 242. 

In each case, the curium was separated from the irradiated americium 
by the following chemical steps : 

(a) Precipitation as fluoride (with a little residual trivalent americium) 
from a solution in which most of the americium had been oxidized to the 
fluoride-soluble hexavalent state. The oxidizing agent used was 
ammonium persulphate with silver catalyst, and the technique was 
similar to that described by Asphrey ef al. (1951). This step also separated 
the curium from most fission products except the rare earths. 

(6) Separation from rare earths by elution from cation exchange 
column with 12 N hydrochloric acid (Street and Seaborg 1950). 

(c) Separation of curium from residual americium by elution from a 
cation exchange column with ammonium citrate solution at 85°c. This 
was for several years the standard separation technique for the higher 
actinide elements, and it is described (inter alia) by Thompson e¢ al. 
(1950). A two-stage separation was necessary to achieve the desired 
purity. The early curium-rich fractions from the first elution were — 
combined and a second elution from a smaller column was carried out. 

Plutonium 238, the a-decay daughter of 162 day curium 242 was 
extracted from purified curium fractions which had been standing for 
several months, by a further ion exchange step. The plutonium fraction 
was purified chemically by conventional solvent extraction and precipita- 
tion reactions on an ultramicro scale. The Pu 238 used in the present 


work was shown by pulse analysis to contain negligible amounts of Cm 242 
and other «-emitters. 
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Sources on aluminium foils were prepared by the tetra-ethylene glycol 
spreading technique (Glover and Borrell 1955) and were assayed by direct 
low geometry «-counting. 

The larger sources, in sealed tubes, were prepared by evaporating a 
known volume of solution to dryness or to a small volume. The solution 
had previously been assayed by taking aliquots with a micro-pipette and 
a-counting in the low geometry counter. 


§3. EXPERIMENTAL ARRANGEMENT 


A diagram of the apparatus is shown in fig. 1. The glass proportional 
counter A, of the Maze type, was ~20 in. long and had an inside diameter 


Fig. 1 


Diagram -of the apparatus. 


of 2in. and a 0-005 in. tungsten wire. The counter A was placed in the 
lead castle B in order to reduce the general background. Lead K 

radiation (~75 kev) is produced by the process of fluorescent oer te 
from gamma radiation having an energy greater than the K a ing 
energy of lead. For this reason the castle was usually lined with 1 mm 


SU 2 
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cadmium sheet C, which has a high attenuation for 75 kev radiation. 
The fluorescent radiation from the cadmium has an energy of 26 kev. 
This, being in general below the energy region of interest, was not 
troublesome. 

The sources D were placed as shown in a box made from 2 mm cadmium 
sheet and lined inside at the bottom with lead. This box had the following 
functions : (a) To provide some collimation so that absorbers E could be 
used, (b) to give a reproducible position for the sources, (c) to minimize 
the amount of backscattered radiation detected in the counter. The 
backscattered Compton radiation is troublesome at these low energies 
because it differs in energy from the primary radiation by only a small 
amount, e.g. by 11-4 kev for primary energy of 60 kev. 

The counter was sometimes used in a longitudinal magnetic field. This 
field, which could be varied from zero to 8000 gauss, was produced by the 
magnet of the large electromagnetic separator at Harwell. 


$4 
4.1. Measurement of the Gamma Ray Energies 


The proportional counter cannot give absolute measurements of energy. 
It is therefore necessary to refer all the measurements to those of standard 
calibration energies. These are most conveniently provided by gamma 
rays of known energy from radioactive sources. It is desirable that the 
error in the measurement of the standard gamma ray should be less than 
tke error in the relative measurements with the proportional counter. 
It is also necessary that the standard gamma ray be not associated with 
other gamma rays of comparable intensity and close to it in energy and 
that the lifetime of the source is long compared with the time taken for 
the measurements. The 59-57-+0-02 kev gamma ray from ?44Am decay 
(Day 1955) fulfills all these requirements very well, but we were not able 
to find any suitable sources giving gamma rays with energies between this 
and that of 191-5+-0-5 kev from the decay of 114"In (Mihelich 1952). The 
accuracy of our measurements therefore depends seriously on the linearity 
of our arrangement. The 59-57 kev gamma ray was used as a standard 
and the 191-5 kev radiation as a check on the linearity. 

A typical spectrum of the 244Am gamma ray taken in a xenon filled 
counter (fig. 2) consists of three peaks. That at the highest voltage 
corresponds to the full energy of the gamma ray and is called the normal 
peak, whilst the other two are ‘escape’ peaks, corresponding to the 
differences between the full energy and those of the Kx and Kf x-rays 
of xenon ; the energies of these x-rays can be estimated to be 29:66 and 
33°61 kev from the table of x-ray absorption edges given by Hill et al. 
(1952). In this particular case the K8 escape peak contains also a 
contribution from the 26-4 kev gamma ray of 241Am. 

For gamma ray energies greater than about 60 kev the Ka and Kp 
peaks are no longer resolved. In this case the difference in energy between 
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the maxima of the full energy and escape peaks is a function of the 
resolution of the counter for the particular energies concerned. A correc- 
tion for this is easy to make since the ratio of intensities of the Kx and 
Kf peaks is constant and the variation of counter resolution with energy 
can be found empirically. The corrections were 0-01 kev and 0:64 key 
for escape peaks having energies of 60 kev and 160 kev respectively. 


3000 
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Pulse height distribution from a 59-6 kev gamma ray taken in a xenon filled 
proportional counter. 


Both the normal and escape peaks were used in the energy measure- 
ments. As the escape peaks are not single lines and therefore do not have _ 
a Gaussian shape, it is not possible to use the method of the Ogive plot 
for determining their position. We therefore defined a peak position, for 
the purpose of these measurements, as the intersection with the curve of 


its median line. 
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Energy measurements of the type described here suffer from errors due 
to the following: (a) non-linearity in the electronic apparatus, (b) non- 
linearity in the proportional counter, (¢) instability in the counter power 
supplies, in the amplifier and the kicksorter, and (@) statistical fluctuations, 

In the original measurements a thirty-channel pulse amplitude 
analyser (A.E.R.E. type 1091A) was used and set up to be linear, with 
an accuracy of one part in 104, by means of a calibrated pulse generator 
(A.E.R.E. type 1091A). The linearity of the amplifier was checked as 
follows. A condenser was charged from the positive d.c. output of the 
type 1091A pulse generator and discharged twice per second by means of 
two oscillating mercury switches. The sharp negative voltage step so 
obtained was fed through a small capacity on to the first grid of the 
amplifier, with the counter still connected. The output pulses from the 
amplifier were arranged to have a spread of about +-0-5 volts by adjusting 
the amplifier noise, and their amplitude was controlled by varying the 
voltage derived from the pulse generator ; they were observed with the 
pulse amplitude analyser. The accuracy of the measurements was two 
parts in 10% and no significant non-linearity was observed in the range of 
output voltage from 16 to 56 volts. The error in this measurement was 
taken into account in assessing the errors of the final results. 

For the more recent measurements on ?42Cm a 100 channel kicksorter 
(type 1363A) and'a type 1405A waveform generator were available. The 
latter produced pulses with adjustable height (defined to 0-1°%), with 
adjustable length and with a repetition rate of 200 c/s; they were fed 
into the amplifier put. It was possible to achieve an accuracy higher 
than that of the earlier measurements with this arrangement. 

The linearity of a proportional gas counter is dependent on two factors 
(cf. West 1953, Chs. [Land IV). The first of these is the magnitude of the 
electron avalanche in the counter ; if this exceeds a certain size then 
saturation sets in. The size of the avalanche is proportional to the 
product of the recoil electron energy and the gas mutliplication. The 
latter can be controlled by the counter voltage, which was therefore kept 
well below the point where saturation set in for these measurements. The 
second factor upon which the linearity is dependent is the variation of the 
mean energy per liberated electron with energy. It did not prove possible 
to show the presence of any significant non-linearity for recoil electrons 
having energies between 11 kev and 60 key ; the accuracy of the measure- 
ments varied from 7 parts in 10° at 11-5 kev to 2 parts in 10° for 60 kev. 
Sources of 1°8Gd and ?41Am respectively were used in these measurements. 
The linearity in the region of energy greater than 60 kev was checked by 
making a measurement of the energy of the gamma ray from the decay 
of "In; values of 191-30-7 kev in the earlier measurements and 
191-3-+-0-4 kev in the later measurements were obtained. These are in 
good agreement with the best reported value of 191-5-40-5 kev (Mihelich 
1952) obtained from measurements on the conversion lines. 

The stability in channel width of the pulse amplitude analyser during a 
measurement was better than 1%, and the stability of a given channel 
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threshold about -+0-025 volts. Since the pulse heights normally used 
were between 30 and 60 volts, errors due to drifts in the pulse amplitude 
analyser were negligible compared with those from other sources. 
Instability in the counter power supplies and in the amplifier caused an 
overall long term variation of about +1°%. The line to be measured and 
the calibration line were always observed simultaneously so that errors 
due to small changes in the overall gain of the system were eliminated. 
When the gamma ray under investigation had an energy greater than 
100 kev it was not compared directly with the standard 241Am line but 
with another line close to it in energy ; this substandard line was then 
compared directly with the ?44Am line or with another intermediate sub- 
standard. 


4.2. Hstimation of the Errors of the Energy Measurements 


The error in peak position was estimated in the following way. First 
a ‘ visual error ’ was defined as one quarter of the difference between two 
visually estimated extreme values of the peak position. This error was 
determined by two observers to reduce the subjectivity of the determina- 
tion, and in general good agreement was found. It is related to the 
statistical error in the measurements, for much more extreme curves can 
be drawn through a set of poimts with poor statistical accuracy than 
through a similar but more accurate set. 

In order to relate the visual error to the standard error of the measure - 
ments, the standard deviation from the mean of fourteen determinations 
of the energy of a 44 kev gamma ray was compared with that calculated 
from the visual errors. The difference between the two results is not 
significant. The standard error of the energy was determined by suitably 
combining the visual error with that of the linearity measurement. When 
the old value of 59-78+0-02 kev (Browne 1952, Asaro, Reynolds and 
Perlman 1952) was used as standard there was a ‘ barely significant ’ 
difference between the energy values obtained by comparison with the 
normal and escape peaks of the ?44Am gamma ray and in this case the 
visual error appeared to be a factor of 1-4 times too small. This difference 
disappeared when the more recent value of 59-57 kev (Day 1955) was used. 


§ 5. MHASUREMENT OF THE INTENSITIES 


With the experimental arrangements used here, it is possible to deduce 
the relative intensities of gamma rays having energies less than 50 kev 
from the areas under their respective peaks and a knowledge of their 
absorption coefficients in the counter gas and in the absorbers. 

The mass absorption coefficients for photoelectric absorption in krypton 
and xenon were obtained as follows. Values for the total mass absorption 
coefficients were obtained from the tables given by Compton and Allison 
(1935); they had to be obtained by interpolation for xenon. After 
correction for absorption by the Compton process they were plotted 
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on a log-log scale against the gamma ray energies. The points lay on 
straight lines so that the photoelectric absorption coefficients could be 
read off for any energy. 

When the ranges of the recoil electrons become comparable with the 
counter radius a fraction of the pulses is lost from the peak by wall 
collisions. In general these ‘lost’ pulses have an energy much below 
the peak value (West et al. 1952). The loss of electrons can be reduced by 
placing the counter in a longitudinal magnetic field (West and Rothwell 
1950). The electron tracks are then curled up around the field direction. 

If an infinitely large longitudinal magnetic field was applied to a long 
proportional counter, irradiated with gamma rays in its central region 
only, then there would be no loss of electrons provided that they could not 
escape out of the ends. Hence, if the law relating the counter efficiency 
to the magnetic field were known, it would be possible to calculate the 
efficiency for any field. A rigorous theoretical calculation of this law 
would be difficult since the effect of scattering, which is considerable 
for low energy electrons, would have to be taken into account. Elementary 
considerations, neglecting the effect of scattering, suggest that the 
efficiency « (H) of the counter for magnetic field H, and energy HL, should 
be given approximately by an expression of the form e(H)= 
[<( 00)—e(0) ](1—f/H)? where f is a constant for a given H and where H 
is so large that the radius of curvature of the electron is smaller than the 

counter radius. From our experiments it appeared that an expression of 
this form was in fact obeyed. In figs. 3A and 4 A are shown the spectra 
of the 192 kev gamma ray from 14JIn taken in a krypton filled counter and 
of the 411 kev gamma ray from 1°’Au taken in a xenon counter with and 
without a magnetic field, whilst in figs. 3B and 4B are plotted the 
square roots of the numbers of counts in the peaks against the inverse of 
the magnetic fields. Similar results were obtained for gamma rays of 
other energies. 

It seems plausible that this semi-empirical relation should be obeyed 
to infinite field. Relative intensities were therefore estimated by 
obtaining a few points on the linear portion of the curve and extrapolating 
to infinite field. 

The results were compared with those obtained with a sodium iodide 
scintillation detector and good agreement was found. The procedure 
with the scintillation counter is direct since there is no significant logs for 
low gamma ray energies provided that the escape peak is included ; 
however the resolution is worse than that of a proportional counter by 
approximately a factor of five. The rather complex procedure for 
obtaining relative intensities with a proportional counter is therefore 
justified when a complicated spectrum is being analysed. 

The numbers of low energy (~44 kev) gamma rays per alpha-particle 
disintegration were determined as follows. The intensity of the gamma 
rays from a source of each isotope was measured relative to that of the 
59-57 kev gamma rays from an *41Am source by counting under the same 
geometrical conditions and correcting for counter efficiencies. By 
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measuring the disintegration rates of the sources and assuming the value 
of 0-40+-0-025 for the number of 59-57 kev gamma rays per alpha- 
particle (Beling et al. 1952) the number of 44 kev gamma rays per dis- 
integration was obtained. The 242Cm and ?41Am sources were thin and 
their disintegration rates for alpha-particle emission were measured in a 
standard low geometry counter. The *°*Pu was in solution and contained 
in a thin walled polythene vessel. Its source strength was found by alpha- 
counting a known weight of solution. 


§ 6. EXPERIMENTAL RESULTS 
6.1. The Gamma Radiation from the Decay of 8®Pug, 


The energies and relative intensities of the gamma ray lines are shown 
in table 1. 


Table 1. Energies and relative intensities of gamma rays from the decay 


Energy of line in kev 43-49 +0-08 99-82 +0-4 153-1 +0-6 
Relative intensity 1-00 (2-7540-25) 10-1 | (3-0+0-4) 10-? 


Two sources of 238Pu, which had been put through different purification 
procedures, were used in these measurements. They both gave the lines 
reported in the table with the same relative intensities. Some weak 
high energy radiation (~750 key) was also observed but this differed in 
intensity in the two sources and may have been due to fission product 
contamination. It was not possible to establish definitely the presence 
of any other radiation, apart from the L radiation of U. In particular an 
upper limit to the intensity of uranium-k X -radiation relative to the 43 kev 
radiation was 3x 10-?. 

The absolute intensity of the 43 kev radiation was found to be 
(3-8-+0-4) x 10-4 gamma rays per alpha-particle. 

These results are consistent with those reported by Asaro and Perlman 
(1954). 


6.2. Gamma Radiation from the Decay of ?4#Cm and *Cm 


As explained in § 2, sources of 4°?Cm in general contain some *Cm ; 
the ratio of 242Cm to 248Cm depends on the irradiation of the source and 
the time after preparation. Sources having different irradiation histories 
were studied and the decay of one source containing an appreciable 
fraction of 24Cm gamma activity was followed over a period of 780 
days. From these measurements and from known data on the beta decay 
of 22°Np, which leads to the same daughter nucleus as does *“*Cm it was 
possible to assign the various observed lines to the decay of ?4°Cm, 
243('m and 241Am, which occurred as a contaminant in some of the sources. 
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Fig. 5 
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Fig. 6 
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Pulse height spectra in the energy range 140-320 kev from a source containing 
242Cm and 24%Cm. 
(a) Taken in a xenon filled proportional counter with a longitudinal magnetic 
field of 7500 oersteds. 
(b) Taken in a sodium iodide crystal. 
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The spectra from the gamma rays with energies between 60 and 170 kev 
taken in xenon and krypton filled proportional counters are shown and 
compared in fig. 5; the normal and escape peaks of the various lines 
are joined by bars. The advantage of using two counters with different 
gas fillings in investigation of a complex spectrum Is well illustrated ; 
the different escape energies of the two gases make it easy to tell which are 
the norma] and which are the escape peaks. : 

In fig. 6A is shown the spectrum of the gamma rays with energies 
greater than 150 kev taken in a xenon filled counter placed in a longi- 
tudinal magnetic field of 7500 oersteds, whilst in fig. 6B is shown the 
spectrum of the gamma rays taken in a sodium iodide crystal. 

The peak at 180 kev in fig. 6A is likely to be the escape peak of a 210 kev 
gamma ray rather than the normal peak of a 180 kev gamma ray for the 
following reasons, (a) the peak at 200 kev is too large to be accounted for 
solely by the escape peak of a 228 kev gamma ray, (b) in figs. 6B and 
8B the intensity of the peak at 225 kev relative to that at 275 kev is 
1:00-0:05. This is consistent with the peak at 225 kev arising from 
210 and 228 kev lines but not consistent with it arising solely from the 
228 kev line, if the relative intensities of the 228 and 277 kev lines are 
as given in table 3, (c) a gamma ray of energy 210 kev is observed in the 
decay of ?2°Np (Graham and Bell 1951) and would also be expected to 
appear in the 248Cm decay (see § 7.2). 

The variation of the intensities of the various lines relative to the 
59-6 kev 241Am line is shown in fig. 7. 744Am has a half life of 470 years 
(Harvey 1952) so that it does not decay appreciably in the time involved 
here. It can be seen that the 44 kev and 158 kev lines decay with a half 
life consistent with that of 162-5 days for 74*Cm (Glover and Milsted 1954) ; 
the 102 kev line decays with a longer apparent half life and the 117 key 
and 210-277 kev lines remain constant within the errors of the measure- 
ments. The half life of 7*%Cm has been reported to be 35 years (Asaro 
1953) so that radiation from it will remain approximately constant during 
the period of our observations. It therefore seems likely that the 117 kev 
and the 210-277 kev radiations arise from from 243Cm decay. 

Gamma ray spectra, taken with a sodium iodide crystal, for two sources : 
(a) with negligible **Cm gamma activity, and (b) containing a large 
fraction of *Cm gamma activity are compared in fig. 8. It can be seen 
that the 102 key line must be present in both the 242Cm and 243Cm spectra. 
This explains the longer apparent half life for the 102 kev line shown in 
gee te 

In table 2 are shown the energies and relative intensities of the gamma 
rays which are attributed to the decay of 242Cm. In addition to these the 
L X-radiation of Pu was observed. The number of 44:03 kev gamma 
rays per alpha-particle was found to be (3-9-L0-7) x 10-4. 

The energies and relative intensities of the lines which are attributed 
to the decay of ?Cm are shown in table 3. The energies of the three - 
higher energy gamma rays are in good agreement with those of 210, 227 
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The intensities of the various lines from a Cm source, relative to that of the 
59-6 kev line from 241Am, as a function of time. 
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Fig. 8 
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and 276 kev obtained by Graham and Bell (1951) from a study of the 
conversion lines from the decay of 2°Np. The two lower energy lines 


Table 2. Energies and relative intensities of the gamma ray lines from 
the decay of 24#Cmg, 


Energy of line kev 44-03 +-06 101-80+-17 157-61+:3 
Relative intensity 1-00 (9+1-1)10-2 (6-0 +:5)10-2 


have energies in good agreement with those of 102-29 and 117-96 key 
predicted for the K« and Kf x-rays of plutonium from the tables of Hill 
et al. (1952). The intensity of the plutonium-K—X-radiation relative 
to that of the three gamma rays was measured with a sodium iodide 
erystal and found to be 2-34-L0-25. 


Table 3. Energies and relative intensities of the gamma rays from 
248C0m decay 


Energy of line |102-15+-2 | 117-30+0-4 | 210+1-5 223-2 a LD 
kev 
Relative 3-6 ---4 1-4 +02 | 0:5 £-15 | 0-65+0-15 1-00 


intensity 


$7. INTERPRETATION 
7.1. The Gamma Radiation from the Decay of 2*Pu and 247Cm 


The gamma ray spectra arising from the decays of #3%Pu and ?#Cm 
are very similar. The assignment of the gamma ray lines to transitions 
between excited states in 724U and *8%Pu has been discussed by Asaro 
and Perlman (1954) Asaro, Thomson and Perlman (1953) and Rose and 
Newton (1953). The proposed decay schemes are shown in fig. 9. The 
two lower energy gamma rays in each case almost certainly correspond to 
transitions between the three lowest states, since their energies correspond, 
within the experimental errors, to the energy differences of these states 
as determined by the alpha spectra. Asaro and Perlman (1954) have 
shown that 153 kev and 100 kev gamma rays from *°*Pu decay are in 
coincidence. This gives strong support to the proposed decay scheme. 
Support is also given by the fact that the ratios of the energies of the 
excited states and the cascade gamma ray transitions are just as expected 
from the strong coupling limit of the Unified Nuclear Model of Bohr and 
Mottelson (1953). 

According to this model, in even-even nuclei well removed from closed 
shells there should exist rotational levels. Their energies, to first order, 
are given by the expression 


2 
B(L)= T_Ii+ 1), J=0, 2, 4even parity 
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where J is the spin of the level and 9 is an effective moment of inertia 
related to the deformation of the nucleus from the spherical shape. The 
rotational levels should decay by cascade #2 transitions to the ground 
state. A correction E to the energies, analogous to that due to the 
rotation—vibrational interaction in molecular spectra is given by 


HI) —2 : nN TAT 1)? 
Cie (a) (5) ofa) 


where wy, is a measure of the vibrational frequencies of the nucleus. 


Fig. 9 
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The decay schemes of ?°8Pu and 242Cm. 


In table 4 are shown the ratios of the energies #, and H, of the second 
and third excited states of °4U and 288Pu to those H, of their first excited 
states. These ratios agree remarkably well with those expected from the 
first order theory. The small deviations can be taken into account by 
assuming values for w,,, of 1:1--0-1 Mev and 1-41-+0-08 mev for 234U 
and ***Pu respectively. These values are to be compared with a value of 


about 1 Mev expected on the hydrodynamic approximati 
ximat 
Mottelson 1953, p. 14). 2 pproximation (see Bohr and 
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7.2. The Decay of 24Cm 


The decay of *48Cm has been investigated by Asaro, Thomson and 
Perlman (1953) and Asaro and Perlman (1954). They found alpha- 
particle groups having energies of 5-985, 5-777, 5-732 and 5-677 Mev and 
photons with energies of 104 key, 226 kev and 278 kev, the latter two 
being in coincidence with the 5-777 Mev alpha group. A number of 
authors (Graham and Bell 1951, Engelkemeier and Magnusson 1955, 
Fulbright 1949, Slatis 1948) have studied the beta decay of 23°N p which 
leads to the same residual nucleus 78°Pu as is reached by the alpha decay 
of Cm. The Coulomb excitation of 28°Pu has been studied by Newton 
(1955). A level scheme derived from all these measurements is shown 
in fig. 10. 

It can be seen that the 277, 228 and 210 key radiations all arise from the 
metastable state labelled 285 kev. If the decay scheme of fig. 10, is 
correct, the predominant part of the plutonium K—X-radiation which we. 
observe must arise from the internal conversion of the 277, 228 and 210. 
kev transitions. 


Table 4. Comparison of experimental and theoretical energy ratios for 
the excited states of 7°4U and 238Pu 


Nucleus E,/ Ey E,/E, 
Bee (yi 3:295 +-010 6-816 + -022 
Pe Pu 3312-005 6-891 +-011 
theoretical 10/3 7-00 


The lifetimes for the individual gamma rays are approximately 10-8 sec. 
whilst the Weisskopf single particle formula predicts lifetimes of about 
10-2 sec and 10-! sec for #3 and M3, 10-® sec and 10-® sec for H2 
and M2 and 10-14 sec and 10-! sec for #1 and M1 radiations. Transi- 
tions which are faster than the single particle predictions by as much as a 
factor of 10? are rarely observed except in the H2 case where collective 
effects are important. It therefore seems safe to conclude that these 
transitions are dipole or quadrupole. . : 

From our figure for the intensity of the K—X-radiation relative to that 
of the three gamma rays and with the assumption that the fluorescent: 
~ yield of the K-shell in plutonium is 0-94 Burhop (1952) we obtain an average 
conversion coefficient of 2:48-0-25 for the three transitions. According 
to Slatis (1948) and Fulbright (1949) the K-conversion lines of the three 
transitions have roughly equal intensities so that the A-conver- 
sion coefficients must all be of the order of 3. The # ! and #2 
K-conversion coefficients for gamma rays of this energy are of the order 
of 0-1 (Rose et al. 1951) hence it follows that the observed transitions must 
be predominantly magnetic in character. If the transitions are all pure 


Seay) 
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M1 the average conversion coefficient should be 3-2=-0-2 whilst if they are 
all pure M2 it should be 7-2--0-5. The errors on these quantities are 
difficult to estimate because the conversion coefficients for the 210 and 
228 kev transitions have to be derived by extrapolation from the calculated 
results of Rose et al, (1951). 

Fig. 10 
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In order to obtain agreement with the experimental average conversion 
coefficient, we have to assume that on the average the transitions are 
either (a) 0-75+0-25 H2 or (b) 0-33 M2+0-67 Fl. Of these (a) is more 
likely since (b) would require M2 transitions about 30 times faster than 
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the single particle estimates. We may obtain further evidence on this. 
by estimating K/L ratios for the three transitions from the curve given by 
Slatis (1948) and comparing them with the theoretical K /L ratios. This 
comparison is given in table 5 ; the L-conversion coefficients were obtained 
from calculated data privately circulated by M. E. Rose. The M1-H2: 


Table 5. Comparison of experimental (SlAtis, 1948) and calculated 
K/L ratios. The errors are our estimates 


Ey kev 209 226 277 | 

___ K/L (expt) 4:440-7 4-7 0-9 6541-5 | 
K/L 0-75 M1 £0-25 £2 5-0 5-4 5-4 
K/L 0-67 E1 £0-:33 M2 2-4 2-4 2.9 


mixture is clearly favoured by this evidence. More detailed measure- 
ments on the conversion spectrum are required to obtain the individual 
mixture ratios for each transition. 

It follows from the above result that the average M1 lifetime is about 
10-° sec and that the average H2 life is about 410-8 sec. The M1 
transitions are therefore highly retarded by a factor of 2104, whilst 
the #2 transitions are retarded by a factor of 20, relative to the single: 
particle estimates. 

It seems likely from the results reported by Asaro and Perlman (1954) 
that the states at 285, 331 and 387 kev in *°°Pu (see fig. 10) belong to a. 
rotational band having K=5/2 or 7/2; the four low-lying states probably 
belong to a rotational band having K=4 (Newton 1955). Magnetic: 
dipole radiation between these bands will therefore be A-forbidden 
(Alaga et al. 1955). This fact may possibly account for the large retarda- 
tion factor of the M1 radiation. The #2 radiation will be allowed if the 
upper band has K=5/2 and K-forbidden if it has K=7/2. The smallness. 
of the retardation factor and the scarcity of data on K-forbidden transi- 
tions makes it unsafe to draw any conclusion from this result at present. 
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ABSTRACT 


The neutron elastic scattering cross sections for lead, cadmium, copper, 
aluminium, oxygen, nitrogen, carbon, beryllium and lithium have been 
measured for an effective neutron energy of 136 Mev. A large liquid 
scintillator was used, and the cross sections were obtained on an absolute 
scale. 

The integrated total elastic cross sections are consistent with previous 
measurements of total and inelastic cross sections. 


§ 1. INTRODUCTION 


HiGuH energy neutrons are expected to show diffraction patterns in small 
angle scattering off nuclei. This has been studied at 290 Mev by Ball 
(1952), 84 Mev by Bratenahl e¢ al. (1950), and at several lower energies. 
The width of the central diffraction maximum depends upon the nuclear 
radius. Consequently a study of these maxima enables determinations of 
nuclear radii to be made. The earlier measurements have not been made 
at angles small enough to give a reliable extrapolation to 0°. On the 
assumption of an optical model, the value of the scattering cross section 
at very small angles can be used to yield a value for the volume integral 
of the nuclear potential which is largely independent of other parameters, 
as discussed in an accompanying paper (Wilson 1956). 


§ 2. METHOD 


The same method was employed as in the n—p scattering experiment 
(Thresher et al. 1955 a) . The absolute differential scattering cross section 
was obtained by measuring the ratio of the scattered intensity to the 
direct beam intensity. 

With the counter at an angle to the direct beam the maximum available 
neutron intensity was employed, while the beam was reduced by a factor 
of about 300 in order to obtain direct beam normalization. When the 
full beam was employed, the monitor counting rate no longer varied 
linearly with neutron beam intensity. This effect, together with a possible 
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change in proton energy spectrum with intensity, were not separated and a 
correction for both was determined experimentally. This correction was 
16+1%, for the two effects combined, and compares with 14+1% found 
in the n-p experiment. We allow 2°/, uncertainty in this. 


§ 3. SCATTERERS 


Lead, cadmium, copper, aluminium, oxygen, nitrogen, carbon, beryllium 
and lithium scatterers were used. They were chosen for the following 
reasons : (1) the values of A are fairly evenly spaced through the periodic 
table, (2) a comparison could be made directly with previously determined 
neutron cross sections (Taylor and Wood 1953, Culler 1955), elastic 
scattering cross sections (Bratenahl et al. 1950, Ball 1952), and inelastic 
scattering cross sections (Voss and Wilson 1956), (3) the elements were 
readily available in adequate purity. 

The area of the neutron beam was 7-6 cm x 2-54 cm, so that the areas 
of the scatterers were chosen conveniently larger than this. The thick- 
nesses of the scatterers were chosen to strike a compromise between a 
conveniently high counting rate and a reasonably small multiple scattering 
correction. The lead, cadmium, copper, aluminium and lithium were cast 
into rectangular forms, while beryllium metal powder was obtained 
already pressed into this form. Liquid nitrogen and oxygen were held in 
the liquid target container as in the n—p scattering experiment (Thresher 
et al. 1955 a). The lead, cadmium, copper, and aluminium targets were 
chemically pure to better than 99-99%. These samples were also 
without flaws. Stereoscopic x-ray photographs were taken of the 
beryllium and lithium to test for internal flaws and impurities. The 
beryllium had several flaws. Iron filing impurities were seen with the 
x-rays and their amount was estimated as 4% by volume. A further 3% 
oxygen impurity was also estimated. The density of the lithium was 
3°/, high compared with the expected value, while x-ray analysis indicated 
about 5°, impurity, probably lithium oxide from oxidization in casting by 
volume. 

The nitrogen and oxygen scatterers were moved into and out of the 
beam manually, while the remaining scatterers were clamped into position 
on a remote controlled target changer. These scatterers, together with a 
blank, could be interchanged in any order. This enabled sets of angular 
runs and the corresponding backgrounds to be determined with a minimum 
time lapse between runs. 

The counter and electronics were the same as in the n—p scattering 
experiment. Some considerable time had elapsed since the n—p scattering 
measurements were made (‘Thresher et al. 1955 a) and, moreover, the beryl- 
lium target for producing the neutron beam for the present experiment 
was placed at a slightly smaller radius in the cyclotron than before. 
Consequently the effective energy of the neutron beam, the effective area of 
the counter, the linearity factor of the counter and monitor, and the pile- 
up correction were all redetermined. 


CARBON 


ALUMINIUM 


CADMIUM 


GOREER 


Differential scattering cross sections for 


C, Al, Cu, Cd and Pb. 
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§ 4, EFFECTIVE AREA OF THE COUNTER 


The geometrical cross section of the counter was 14 x 7-6 cm. Owing to 
edge effects, the effective area of the counter was less than the geometrical 
cross section, and this had to be determined for the absolute cross section 
measurements. 

A uniform neutron flux was made available at a distance from the counter 
The neutrons which passed through a precision 7-6 x 2-54 cm steel colli- 
mator in a steel well were counted by the counter and then neutrons 
counted without the collimator. The ratio of counting rates gave the 
ratio of the effective area of the counter to the 7-6 x 2-54cm area. The 
effective area in the present experiment was measured as 72:313-0 em? 
as compared with 73-0-+1-3 cm? at a slightly higher effective energy in the 
n-—p scattering experiment (Thresher et al. 1955 a). 

The effective area measurement only applied when the area of the 
neutron beam striking the counter at the 0° normalization position was 
exactly 7-6 x 2:54 cm?. This was very nearly the case for the 6 ft position. 
At 12 ft however, the area of the beam was 25°% larger, and the counts 
in the direct beam were 8°, lower. Accordingly the 12 ft points were 
normalized to the 6 ft points. Further, the edges of the beam were not 
sharply defined at the 6 ft position, so that there was an uncertainty in the 
application of the effective area at this position. A correction and error 
of (8--2)°% was made on the direct beam normalization to allow for this 


§ 5. Prrz-Up 
Pile-up phenomena have been fully discussed in an earlier paper 
(Thresher e¢ al. 1955 b). Pile-up corrections of counts were only significant 
for normalization runs and for scattering off the heavier elements at smaller 


angles. In no single count did this correction exceed 3°/ in the present 
experiment. 


§ 6. EFFECTIVE ENERGY 


The spread of neutron energies accepted in this experiment was about 
20% full width at half maximum. Had the biases been set at higher 
values in order to reduce the energy spread, the counting rate would have 
been reduced and become extremely dependent on electronic stability. 
The * effective energy ’ was determined by measuring the total attenuation 
cross section of neutrons by carbon in good geometry (‘Thresher et al. 
1955 a, b) since this is well known as a function of neutron energy (Taylor 
and Wood 1953, Culler and Waniek 1955). A value of 136-+-4 Mev was 
obtained by this method, assuming that the cross section varies little 
with energy. 

At angles near the diffraction minimum the cross section can vary as 
fast as 1/H* and an effective energy as low as 129 Mev was possible. For 


this and other reasons little accuracy was attached to the points near the 
diffraction minimum. 
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The error quoted in the effective energy is +4 Mev. Of this, +3 mev 
was due to the uncertainty in Taylor and Wood’s energy scale, while the 
rest was due to statistical errors in the effective energy determination. 


§ 7. BACKGROUND CORRECTION 


The neutron background for the counter set at an angle to the neutron 
beam was determined in a separate run for each angle. At small angles 
the background decreased rapidly with increasing angle, but at larger 
angles was nearly constant. 


Fig. 2 
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Differential scattering cross sections extrapolated to 0° plotted as a 
function of A. The lines are do/dQ proportional to A*/? and A’. 


At small angles the background was appreciably affected by the presence 
of the scatterer. This has been discussed in the n—p scattering paper 
(Thresher et al. 1955 a). The background consisted of two parts. The 
first was due to neutrons which had filtered through the shielding direct 
from the cyclotron target. This was approximately constant irrespective 
of counter position. The second part was due to small angle scattering of 
the direct beam, partly by the air and partly by the collimators. This 
was appreciable only at small angles. The first part was not affected by 
the presence of a scatterer, whereas the second part, due to the arrangement 
of the shielding, had to pass through the scatterer and consequently was 
attenuated. This attenuation was calculated from the total cross sections 
of Taylor and Wood (1953) for eachelement. The value of the background 
at the 10° position was assumed to be the constant background, and hence 
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the corrected background could be calculated. In no case did the 
attenuation of the background modify the scattering cross section by 
more than 3%. 

The background due to small angle air scattering was reduced by 
inserting a long polythene bag containing helium in the beam path, in 
order to displace the air. 

Since at large angles the presence of the scatterers did not affect the 
background appreciably, a straightforward subtraction of the background 
gave the necessary correction. 

The uncertainty in this background subtraction procedure was small, 
and was included in the errors. 


§ 8. ANGULAR RESOLUTION 


In order to determine the angular resolution of the counter the following 
procedure was adopted. The counter was moved across the beam, keeping 
the axis parallel to the beam, and the counting rate was plotted as a 
function of position. The resulting curve is a combination of the angular 
resolution of the counter and the size of the beam. This was both to 
find the angular resolution and the effective centre of the counter. The 
full resolution width at half maximum was found to be 2° and 1° when the 
counter was 6 ft and 12 ft respectively from the counters. 

No correction for the angular resolution of the counter was necessary 
where the rate of change of elastic scattering cross section with angle was 
constant. However, where this varied rapidly with angle, corrections 
were necessary. The correction became greater than 3°, only for the 
scattering cross section for lead for angles larger than 7°. 


§ 9. MULTIPLE SCATTERING 


It was necessary to apply a correction to the data to account for multiple 
scattering of the neutrons in the scatterers. The method used to evaluate 
this correction has been described by Hildebrand (1951). A neutron which 
is not absorbed makes on an average ¢/A, collisions in passing through a 
slab of thickness t, where \, is the mean free path for neutron elastic 
scattering for the element. The probability of just » elastic collisions is 
given by the Poisson distribution 


fap We t 
aE exp(—t . 


If this is multiplied by the probability that the neutron is not absorbed in 
the thickness ¢ then the probability of n elastic collisions, and no other 
collisions, in passing through the slab is 


1 ita\G t t Le /e\t t 
Teg) i ta)l-ale) o(—) 


where \, is the mean free path for inelastic collisions, and dp is the mean 
free path for a scattering event which is either elastic or inelastic. The 


above expression assumes that only small angle scattering is involved, so 
that ¢ is unique. 
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The central maximum of the scattering pattern may be approximated 
by a Gaussian distribution. Since transverse momenta comprise a 
‘random walk ’ problem, the customary assumption can be made that the 
angular spread after scatterings is \/n times the spread from a single 
scattering. This implies that the forward scattering intensity is reduced 
by a factor 1/n. 

Thus the relative number of neutrons scattered into unit solid angle at 
an angle @ to the direct beam by a slab of thickness ¢ is given by 


= S Alas) o(~3) Ilse, t0(za)} 


In calculating this correction it was found that only double scattering 
contributed appreciably. In the region of central diffraction maxima this 
correction was about 4%. 


§ 10. Swati ANGLE Macnetic Moment—Covutoms Frieup SCATTERING 


The elastic scattering cross section at small angles includes, in addition 
to the nuclear scattering, a contribution due to the interaction of the 
magnetic moment of the neutron with the Coulomb field of the nucleus 
(Schwinger 1948, Voss and Wilson 1956b). A calculated value was 
subtracted from the measured cross section to yield the nuclear scattering. 


§ 11. Toran Exastic ScaTrEeRtneG Cross SECTION 


The total elastic scattering cross section o, was obtained by graphical 
integration of the differential cross section for scattering of neutrons into a 
forward cone extending to the largest angle at which measurements were 
made, namely 173°. 

The intensity of neutrons scattered into wide angles (greater than 173°) 
by the heavy elements is so low that measurements become very inaccurate. 
However, because of the sine factor which weights the differential cross 
section, the contribution due to these neutrons may not be negligible. 
Consequently, the wide angle scattering was estimated using the differential 
cross section measurements made with protons (Dickson ef al. 1955). 
The effect of the wide angle scattering is more serious for the lighter 
elements, and hence the measured total elastic cross sections for these 
elements are less accurate. The measured values of c, are given in table 2, 
together with the contribution due to wide angle scattering, and an 
estimated error. 

The differential scattering cross sections were expected to vary as 1/H 
and the effective energy of the total cross sections became 134+-4 Mev. 


§ 12. CONTRIBUTIONS FROM INELASTICALLY SCATTERED NEUTRONS 


The energy spread allowed in this experiment permitted an appreciable 
number of inelastically scattered neutrons to be counted. Accordingly 
it was necessary to determine the possible inelastic contribution to the 
elastic scattering cross sections. 
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An indication of the magnitude of this contribution for lead was obtained 
by a comparison of the upper and lower bias determinations of the cross 
sections. A diffraction minimum as exhibited by lead, would be obscured 
by inelastic events if the latter were present in appreciable amounts. 
Hence, assuming that the upper bias recorded elastic events only, a 
comparison of the values of the cross sections determined by the two biases 
gave the contribution due to inelastic events. By this method, the 
inelastic contribution was less than 20% of the effect in the region of the 
apparent lead minimum. Assuming the inelastic contribution to be 
isotropic, this formed less than 0-4°, of the effect at small angles. As the 
upper bias might have been recording a small number of inelastic events an 
error of 1°, was allowed for inelastic scattering. 

For carbon a larger error might be present. There are two ways of 
arriving at an estimate of the error involved. Experimentally the ratio 
of upper to lower bias determinations could be compared for lead, where 
the effect is small, and for carbon. This shows that 5+-5°% of the counts 
at small angles from carbon came from inelastically scattered neutrons. 
Theoretically, the data on inelastic proton scattering for carbon (Strauch 
and Titus 1955) can be extrapolated to zero and assumed to apply to 
neutron scattering. The relative efficiency of counting the inelastic 
neutrons can be estimated from the bias curve. From this calculation 
3-+2°, of the small angle counts are due to inelastically scattered neutrons 
and about 10-+5° at 15°. No correction has been made for this error, 
which is expected to vary as 1/A. 


Note added in proof.—Recent data from Upssala suggest that at smail 
angles the contribution of inelastic scattering is sensibly zero. 


Table 1 
Elements G Al Cu Cd Pb 
co, up to 17 1/2° : 123 St 691 1176 1911 
o, beyond 17 1/2° 36+10 37-210 52+14 69-+17 117 +30 
o, 159+10 | 348+10 | 743+14 | 1247+17 | 2028+30 | 
1 —Oinel 144+11 | 334+17 | 788+27 | 1336+30 | 2050-+50 | 
Discrepancy (%) 10+11 4+6 —613:5} —7+3 —1+-3 


Weighted average discrepancy —3+2%. 
Units of 10-2? em?. 

It can be seen from table 1 that for light elements there is a discrepancy 
between the integrated elastic cross sections and the value of PES 
This discrepancy is within the quoted error, but might disappear if the 
correction is made for inelastically scattered neutrons. 


§ 13. ERRors 
‘Two types of errors may be distinguished. 
Hs - : ; 
Firstly, relative errors on each point. In every case, except lithium 
and beryllium, the relative errors quoted on all points were statistical. 
Due to impurities, an additional error of 4°, and 8% was attached to the 
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small angle lithium and beryllium points respectively. In addition there 
is the error due to inelastically scattered neutrons discussed in § 12. 
Secondly, errors on the absolute scale. These absolute errors were 
determined by two independent methods : 
(i) the direct beam normalization procedure, in which the error was 
6% and included the uncertainty in the effective area, linearity factor of 
the counter and monitor and contributions from inelastic scattering; 


(ii) a comparison of the total elastic scattering cross sections with the 
value of o7-c;,, from the measurements of Taylor (1953) and Voss (1956a). 
This comparison is shown in table 2 and yields a value accurate to 4%. A 
discrepancy of 3% exists between the total cross sections determined by 
the two methods. 

The direct average of methods (i) and (ii) has been taken and an error 
of 4° is assigned to the absolute scale. 

Only relative errors have been drawn on the cross sections for the angular 
distributions (fig. 1 (a to e)). The angular uncertainties of these distri- 
butions were due to the finite thickness of each scatterer. The relative 
and absolute errors were combined, however, to give the quoted errors on 
the total cross sections and forward scattering cross sections. Theoretical 


Table 2. 


Differential Cross Sections in units of 10-78 em?2; an additional 
4°, absolute error should be assigned 


6 Cc Al Cu Cd Pb 
0° 1:190+0-06 | 4:15+0-2 |14:60+0-5 |33-00+1-4 | 73-10+3-2 
aes 1:180+0:05 | 4:08+0-2 |14:30+0-4 |31-40+1-2 |71:10+2-8 
2.6 1-:140+0-04 | 3:87+0-1 |13-50+0-4 |29-50+0-5 |65-40+1-3 
Se 1-:150+0-04 | 3:50+0-08 | 11-30+0-3 | 27-60+0-5 | 55-90+1-2 
4° 1:010+0-03 | 3-47+0-07 | 10-90+0-2 | 23-40+0-4 | 45-30+0-9 
5° 0:970+0-03 | 3:19+0-06 | 9-40+0-2 |18-30+0-4 | 30-20+0-6 
6° 0-790-+0-02 | 2-66+0-08 | 7-38+0-25]14-50+0-4 | 18-70+0-5 
i 0-780+0-02 | 2-31+0-08 | 5:97+0-25}) 9-00+0-4 | 11-30+0-4 
8° 0-680 +0-02 | 1:86-+0-06 | 4:26+0-20|] 5-70+0-22} 5-70+0:3 
10° 0:550-+0-02 | 1:37+0-04 | 2-08+0-08| 2:11+0-06| 1-20+0-06 
iN ierg Bays 0-70 +0-13 
12° 30’ | 0-390+0-02 | 0-71+40-02 | 0-70+0-03} 0-52+0-04] 0-97 +0-06 
13° 45’ 0:93 +0-06 
15° 0-226 -+0-007 | 0:39-+0-02 | 0:27+0-03} 0-33+0-03| 0-99+0-07 
17° 39’ | 0-138+0-009 | 0-17+0-02 | 0-19-+0-03| 0-40+0-03} 0:53 +0-06 


8 Li Be N, oF 
: 0520-40-03 | 0-760-0-05 | 1-72040-08 | 2-190+0-09 
2°12 | 0.51020-02 | 0-7304.0-05 | 1-630-.005 | 2070-40-06 
4° 0-500-0-02 | 0:700-40-04 | 1-480+0-05 | 1-8004.0-05 
15° 0-140-.0-010 | 0-194-40-010 | 0-270+0-02_ | 0-288 +-0-009 
17°39’ | 0-095-£0-010 | 0-11340-007 | 0-138£0-007 | 0-172-£0-009 


me ger 
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curves for a square potential well are fitted to the data with arbitrary 
normalization on both axes. Three curves fit well, but the carbon and 
cadmium curves if fitted at zero degrees, fall off the experimental points. 
This may indicate a systematic error in these elements. 


§ 14. ExTRAPOLATION TO ZERO. 


The points were extrapolated to zero by using the 1° and 2° points and 
assuming that the cross section varies as 1-«6 in this region (except for 
cadmium where a fit to the whole curve was used). 

The nitrogen, oxygen, beryllium and lithium data were extrapolated to 
zero assuming corrections intermediate to those for carbon and oxygen. 
The results are included in table 2. 

A transition from the transparent nucleus formula (A?) to the opaque 
nucleus (44/3) is seen on proceeding from light to heavy nuclei. 


§ 15. N—-P DIrrERENTIAL CROSS SECTION 


In the work of Thresher e¢ al. (1955 a) preliminary results of this experi- 
ment threw doubt upon the normalization and a large error was assigned 
to the n-p scattering. The discrepancies were found in the effective area 
measurement and are now removed, and we believe that the differential 
cross section scale for the 137 Mev n—p scattering is accurate to 5% if the 
points are raised by 2%. It is possible that the similar, but smaller 
discrepancy at the 105 Mev points is due to the same cause. We cannot 
be sure about this, so we believe that the quoted uncertainty of 6°, should 
be retained. 
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ABSTRACT 


Neutron diffraction scattering at high energies is discussed and it is 
shown that two parameters can be derived from the data; the nuclear 
radius and the nuclear potential, both with little effect of other para- 
meters. These parameters are derived for each nucleus individually. 
The radius, so derived, is larger than the charge distribution radius. The 
potential agrees with that derived from the analysis of total cross sections 
of all elements simultaneously when a nuclear potential with a reasonable 
taper is assumed. 


§ 1. INTRODUCTION 


Neutron diffraction scattering measurements have been made at 84 Mev, 
136 Mev and at 290 Mev and at several lower energies (Van Zyl et al. 
1956 preceding paper, and references therein). The purpose of this paper 
is to discuss their interpretation on the optical model of the nucleus, using 
the Huyghen’s principle approximation for computation (Fernbach et al. 
1949). In particular two experimental features will be discussed, the 
value of the differential cross section extrapolated to zero angle, and the 
width of the diffraction pattern. 

In this paper a departure from the usual course will be made. Hach 
nucleus will be considered independently upon the optical model, and 
only later will the parameters of two nuclei be compared. This is possible 
for this paper, where more than one measurement is considered for each 
nucleus. It was impossible for previous optical model comparisons (e.g. 
Taylor 1954) which considered only total cross sections. In this respect, 
the work may be regarded as a continuation of that of Voss and Wilson 
(1956). 

Fernbach’s original equations can be generalized to include effects of a 
diffuse nuclear surface, and a spin-orbit coupling. This generalization 
has been discussed before (Malenka 1954, Wilson 1955, Voss and Wilson 
1956). The notation used here will be the same as that used in the last 
nae ee een ees ee ee ee 
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two references. The scattering amplitude for scattering from a potential 
(Ltée) Vin)+8 : LU(r) may be written in the form 
f(@) +e . ng(@) 


where f(@)=ik |” pap {1—exp [—2i84(9)] eosthp®, (a) olkP®) (1) 


(@)=ik |” pdp {exp [—2id4(p)] sin(kydr()}Ikp®) - (2) 


where Bu(p)=— ou (1446) le Vine. > 
FO ed ge 
and pias p?- 28, <0) a ee eee cue 3 Ol) 


k is the wave number of a neutron of energy #. 

In the analysis which follows the contribution to U(r) from the inter- 
action of the magnetic moment of the neutron and a coulomb field 
will be ignored, and the appropriate correction made to the experimental 
results, so that, when interpreted, they will yield the nuclear potential 
without the magnetic moment interaction. A correction of about 1% at 
the small angles has already been subtracted from the experimental points — 
in the work of Van Zyl. 

U(r) may now be assumed to have a reasonable form which falls to zero 
outside the nucleus. (U(r)=0 for r>R). Then g(6) varies as @ for small 
@ and if the experimental neutron scattering cross sections are extrapolated 
to zero angle, the value determined, which will be called J (0), is deter- 
mined theoretically by f(0) alone. The angles over which the measure- 
ments are made are small enough that 2 sin }@=@ within experimental 
error, so that Huyghen’s principle should be valid. 


§ 2. VALIDITY 


The approximation is only expected to be valid, even for small 0, if V/H 
is sufficiently small. Shapiro (1955) has compared the angular distribution 
of neutron diffraction sattering using both the small angle approximation 
and an exact phase shift analysis. He finds that J(0) is larger for the 
exact calculations, and the value of J(@)/I(0) is reduced by about 10°, in © 
the region of angles of interest here. This alters the value of R deter- 
mined by the width of the diffraction pattern by 14%. For Shapiro’s 
calculations V/H was of the order of 0:25 whereas V/H here is approxi- 
mately 0-12. Moreover the important calculations of this paper assume 
a tapered potential whereas Shapiro used a square well potential so that 
neglect of reflection at the nuclear boundary is justified. Accordingly 
the correction due to the use of the Huyghen’s principle approximation 
is less than 1% and is neglected. It is of interest to note that the exact 
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calculations quoted give a value for o7 about 10°% greater than the optical 
model. Ifo p were used to derive a value of R, there would be a 2% error 
instead. of the $°% here. 

In the determination of V from small angle scattering, the experimental 
errors are also sufficiently large to mask any errors. 

The neutron diffraction cross sections are determined experimentally 
in the frame of reference where the nucleus is at rest ; and so are the total 
neutron cross sections. Thus Re f(0) is determined in this frame. It 
would be possible to transform Re f(0) into the centre of mass system. 
The energies would now be different for the different nuclei, and would 
fall from 136 Mev to 120 Mev for carbon. The question naturally arises, 
in what frame of reference are the optical model parameters expected to 
be most nearly constant from nucleus to nucleus? So long as attention 
is confined to the parameters R, k, and K, where k,+i3K is the change 
in wave number on entering the nucleus, then there is no difference so 
long as Huyghen’s principle is valid. The situation is different when it is 
desired to derive a value for the potential V. This question can be 
answered by considering the derivation of the nuclear potential from 
individual nucleon-nucleon interactions (Jastrow 1951, Mandl and Skyrme 
1952.) The same problem arises in the interpretation of nucleon—deuteron 
scattering (Chamberlain and Stern 1954). The nuclear potential may be 
derived from the formula 


V = So (ALB mp0) +4 Fnp(0)]-+(A—Z)[Efnn"(0) +4 fnne(0)] - (6) 


which is derived by expressing the forward scattering amplitude for 
neutron—nucleus scattering as the sum of singlet and triplet neutron— 
proton ( fnps|(0)) and fnp(0) and neutron—neutron ( fnns(0) and fnn'(0)) 
scattering amplitudes (with the spin-flip terms removed). This form is 
valid in the impulse approximation and from its structure clearly applies 
to the frame of reference in which the nucleus is at rest. Thus the 
appropriate frame of reference to analyse the zero angle cross section is 
expected to be that in which the target nucleus is stationary. This is the 
procedure that has been adopted in the past. This argument must be 
used with caution, however. The same approximation would show that 
the nuclear potential has the same radius as the matter distribution 
which will be shown to be incorrect. 

Tt remains to be shown that the use of a potential to describe the inter- 
action of the neutron with a nucleus is sufficiently accurate. No exact 
analysis is attempted here, but the problem is related to that of X-ray 
scattering (e.g. Compton and Allison 1935). Here Born approximation is 
valid for the scattering from an atom, and the form factor calculated is 
the same as would be calculated upon the optical model using a potential. 
But in this case, the form factor must be multiplied by the cross section. 
for scattering from an individual electron, which is assumed isotropic in 
the optical model. 
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By analogy it would appear that the calculations using the optical 

model must be multiplied by the square of the scattering amplitude 

Z [3fnp'(0)+4fnps(6)] +(A—Z) [Bfnn'(6)+ afnns(6)] .  - (7) 
The fnn(6) are probably nearly constant with 6 near 0°, but the fnp(@) vary 
with @ (Thresher et al. 1955). For light elements we are concerned with 
fnp to 15°, where Inp (8)=0-7 Inp(0). For lead we are only concerned 
with fnp to 6° where Inp (@)=0-95 Inp(0). 

The internal momentum of the nucleus in the nucleus will reduce the 
effect somewhat since it spreads out the angles involved. Thus an internal 
momentum of 20 Mev will spread the angles over 21°. Experimentally, 
Chamberlain and Segre (1952) has found an angular spread of 45° (at half 
maximum) for quasi p—p collisions in Lithium nuclei, using 340 Mev 
bombarding energy. At 135 Mev this spread should be increased to 70° 
which reduces the effect to insignificance. 

If this is regarded as a correction to the calculations of this paper and 
the momentum distribution is ignored, then the radius determined for 
lead must be multiplied by 0-985 and that for Lithium by 0-86. Due to 
the uncertainties involved, no correction has been assigned to the data, 
but this consideration should be considered as giving the order of accuracy 
to be expected from the calculations. 


§ 3. WipTH OF DIFFRACTION PATTERN 


For a nucleus which is nearly opaque to the incident neutrons, the 
diffraction pattern has the angular distribution 


J,(kR8) 


kk 


Transparency of the nucleus broadens the diffraction pattern slightly 
and reduces the diffraction minima associated with eqn. (8). But in the 
limit of a completely transparent nucleus (Born approximation) the 
value of 6 for 1(0)/I(6)=4 is only increased by 10% from that for an 
opaque nucleus. The diffuseness of the nuclear boundary also has an 
important effect upon the diffraction pattern, filling out the minimum. 

Here, we choose the angle at which I(0)/7(#)=4 and endeavour to extract 
from this value information about the nuclear radius. This choice is 
influenced by the fact that at such an angle the plot of (@) versus @ is 
steep, so that he accuracy of determination is quite large. Also, it is the 
hope that at this angle the experimental errors arising from the detection 
of inelastically scattered neutrons (discussed at some length by Van Zyl 
1956) are small. It has only been thought worthwhile to carry through 
this analysis for 136 Mev where the best data are available. The other 
data scale to the 136 Mev data by the relation 6?/E=constant but the 
errors are considerable. 

When the effect of the diffuseness of the boundary is taken into con- 
sideration, it is found using Born approximation that the value deduced of 


(72)= | Virydr / | V(r)dr oN ee a 


2 


1(0)= eR, a ee 
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is not appreciably affected by the degree of taper used. This result has 
already been found for electron scattering. 

For reasonable optical model parameters the difference is not large and 
can be readily computed. Therefore the most accurate radius parameter 
deduced from our data is the mean square radius. The detailed calcu- 
lations support this conclusion. The angular distributions of Van Zyl 
are in agreement with the Born Approximation shape, as may be seen 
from the figures of the previous paper, but addition of a taper slightly 
improves the agreement. 


§ 4. SPrIN—-ORBIT CoUPLING 


A correction is applied for the presence of spin-orbit coupling in the 
nucleon—nucleus interaction which broadens the angular distribution. 
This correction is semi-empirical and is derived as follows. 

The unpolarized cross section is given by 

es aa et a (10) 
and the polarization in the scattering is given by 
__ 2{Lmf(0)Lg(6)-+ Ref(@)Reg(6)} 


PO) = (11) 
For light nuclei within the central diffraction maximum, 
Imf(0)/Re f(0)=Imf(9)/Re f(0)=A ee 62) 


where A is close to unity. g(@) is mainly imaginary but its real component 
is of such a sign as to reduce the polarization (e.g. Wilson 1955). Then 
an approximate minimum value of g(@) becomes 
ee et ee EAS el yes 3. (13) 
The values of P(@) for the various elements are taken from the data of 
Dickson e¢ al. (1955). These data are for protons and it is assumed that 
the plot of P(@) is the same for neutrons and protons except for the 
coulomb interference region. Here the data are extrapolated assuming 
P(@) is proportional to 6 (Wilson 1955). Thus for carbon at 14°, which is 
the relevant angle, we find |g/f|=0-14 so that the correction to [(@) is 2%. 
If appreciable phase shift occurs in g(@) (which will only occur if a spin— 
orbit potential is assumed which is not the derivative of the nuclear 
potential) then the correction might be up to 5%. 33%+1)% has been 
taken as the correction. For heavy nuclei there is certainly some phase 
shift occurring ; the correction on the above basis at 6° becomes 0:02% ; 
if phase shift occurs the correction could be 1% which is still negligible. 
There remains the effect that the spin- orbit coupling has on f(#). The 
effect here is to reduce the contribution that large p make to Re f(@), and 
to increase the contribution of large p to Imf(6), making the nucleus more 
nearly opaque.t The effect on a determination of the nuclear radius is, 


therefore, small. 


eae ee a Nee SIA I PTT ee” Re Wane ese rar eEn 
+ A discussion of the effect on /mf(0) and therefore op in Wilson (1955) 
contains an error which did not, however, alter the arguments of that paper. 


The discussion here is correct. 
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In order to make an exact estimate of the effect of the spin-orbit 
coupling a magnitude and share of the potential U(r) must be assumed. 
There is by now considerable evidence that U(r) must be nearly pro- 
portional to 1/r{dV(r)/dr} (Wilson 1955, Sternheimer 1955), although an 
imaginary component of U(r) has been assumed by Heckrotte (1956). 
Accordingly a value of U(r) which gives a good fit to the data on polari- 
zation at small angles was assumed. It is not possible to do this con- 
sistently for a square well potential V(r) since with such a potential it is 
not possible to fit the polarization. So that the same U(r) was used for 
the square well as fits the data using a tapered potential. 


§ 5. TAPERED POTENTIALS 


Since a square potential well is unrealistic, the calculations were 
performed for tapered potential wells. The value of R determined by 
the data has been calculated using tapered potential distributions with 
parameters for the potential and absorption coefficient that fit the 
absorption cross section data, the small angle cross sections, total cross 
sections, and the polarization at small angles. 

For lead the potential is assumed constant to a value 0° 8hy,, tapering 
linearly to 1-2Ry. For copper potential assumed is constant to 0-71Ry, 
tapering to 1:29Ry. The total cross sections have been calculated for these 
potentials for a range of parameters (Voss and Wilson 1956), though 
it should be noted that the notation here is slightly different. For carbon 
where the details of the taper become important, a Woods and Saxon 
(1954) potential is assumed with a shape 


= Vo 
"= Thexp (Rana 


(14) 
where a=0-39 fermis. 

All these potentials have the common feature that, for the appropriate 
nucleus, the distance for the potential to fall from 90% to 10% of its 
central value is 2-5 fermis in agreement with the taper of the charge 
density distribution (Hofstadter 1956). 

The taper for lead is less than that assumed by Williams (1955) at high 
energies, but seems the most that can be allowed by total cross section 
measurements (Voss and Wilson 1956). 

For each case the value of (7?) has been calculated and the radius 
Req of a square potential well that will yield the same <r?). This Req 
should be nearly as accurate as the experimental data. 

The value Ry is dependent upon the taper. A variation of the taper 
distance by 0-1 fermi gives less that 1/2°/, change in R,, for heavy nuclei 
and about 1% for light nuclei. No significant error is expected from the 
choice of the other parameters which come from experiment, except the 
polarization which can give 1°, uncertainty for light elements. An over- 


all uncertainty of 5% is therefore applied to Ry for the light elements 
and about 3% to the heavy elements. 
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it was not necessary to carry through detailed calculations for all nuclei 
since the corrections are small. The corrections were computed for C, 
Cu and Pb and applied to the other nuclei by graphical interpolation. 


§ 6. Dara AND COMPARISON 
Table 1 shows the experimental data of Van Zyl and the calculated 
nuclear radii. The four columns give the values according to the following 
assumptions. (1) Born approximation (shown as R,,), (2) Opaque nucleus 
(Roy) (3) A square well potential with Taylor’s parameters for V and K 
and the polarization fitted by a tapering nucleus (shown as Rgy), (4) and 
(5) the * best ’ fit being that calculated from the nuclear potentials just 


Table 1. Angle for which /(0)//(0)=4 for 136 Mev and values of 
R(in fermis) (k= 2-654 x 1018 cm—!) 


Li Cc 


Angle 15°5 14-9 13-9 
(degrees) +0-25} +0:2 | 40-2 


R(B.A.)| 3:48] 3-62| 3-88 
R(O.N.)| 3-07} 3-20] 3-43 
R(S.W.)| 3-41| 3-55] 3-76 


Eu : 2-65 2-80 3-06 
40-2 | +0-2 | +0-15 


REQ 3°37 3°48 3°72 
40-1 | 40-1 | 40-07 


An additional 14% error is applicable to all the values because 
of the uncertainty in the energy. 


described yielding values Reg (the radius of a square well which would 
yield the same r.m.s. radius) and Ry (the radius of the point at which 
the potential falls to one half of its maximum value). 

Figures 1 and 2 show plots of R versus A1/3 for Reg and for Ry. The 
line is the radius of the electron density distribution as determined by 
Hofstadter. 


§ 7. ANALYSIS OF ZERO ZNGLE CROSS SECTIONS 


I(0) may be written |Ref(0)?+[Jm/(0)P. This suggests at once the 
use of the cross section theorem to determine I[m/(0). This gives 
Imf(0)=kog/47. og is known accurately from a variety of experiments 
over the energy range of interest (Taylor 1952, Culler and Waniek 1955, 
Nedzel 1954), so that |Zmf(0)]? may be determined to about 2%. The 
accuracy of the energy measured for op and for (0) does not enter into 


fermis 


The mean radius as deduced in this paper. The line is R=1-35 AY®. These 
values depend upon the taper assumed as discussed in the text. 


Fig. 2 


Req, 


fermis 


The ‘ effective’ cath , 
e peer radius Req of a square well that would produce the same 
ae, ennes The line is 1-2+1-18 413, All points have an additional 

© uncertainty due to the energy determination. 
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this determination to any significant extent for the data of Van Zyl, 
since for that data, the effective energy was determined by a total cross 
section measurement. A slight additional uncertainty arises, however, 
in interpreting the data of Bratenahl and Ball. 

The use of this cross section theorem is more general than the optical 
model or any particular computational approximation, since it may be 
derived entirely from a conservation law. Consideration must, however, 
be given to the effect of exchange interactions of neutrons with the 
neutrons in the nuclei, which is neglected upon the optical model, and is 
tentatively neglected in this paper. 

The Huyghen’s principle method of calculation is not accurate for very 
low energies. By this method of subtracting |Jm/(0)|? from 1(0) to obtain 
[Ref(0)|? we obtain Ref(0) with no uncertainty, and any analysis of 
fef(0) becomes more accurate. This is seen particularly in an analysis of 
the neutron scattering from lead at 84 Mev where Ref (0) is small compared 
with /mf(0) and could not be determined at all within the uncertainty of 
an optical model fit to /(0). Here, however, it is possible to determine 
Ref(9) even in this unfavourable case. 

From (2) we can write 


Ref (0)=k | pdp sin [25,(p)] cos [kpS,(p)]. . . . (15) 
0 
This can be expressed in the form 
bens 
Ref (0) = ay (FV, #, R, U)}, ae LG) 


where F is a function of V, H, and U, which is less than unity but reaches 
unity in the limit of the Born approximation. The term cos [kpd,(p)] 
acts, for a derivative type spin-orbit coupling, to reduce the integrand 
for large p where 5)(p) already vanishes. Thus even though the spin— 
orbit coupling and the imaginary parts of the potential are not well known, 
this measurement can still be used to determine V R?. 


§ 8. MEASUREMENTS OF ZERO ANGLE CROSS SECTIONS 


In table 2 are the experimental results extrapolated to zero degrees. 
Also tabulated are the derived values of |Ref(0)|?. No cross section 
determinations were available at the exact energies required, so that 
Imf(0) was determined by interpolation between different energies. 
Cross section determinations were available for the elements Li, Be, N 
and O at very few energies, so that there is a small error in the interpolation 
at these energies. This has been estimated and included. The author 
has arbitrarily assigned errors to the data of Bratenahl and of Ball. we 
errors at 84 Mev in |Ref(0)|? range from 30%, for Be, C, Al, to about 108 hs 
for Cu, Ag and Pb. At 136 Mev errors range from 15% for Li to 30% io 
Cd and 50% for Pb. At 290 Mev errors are 40% for C and Al and 100% 
for Cuand Pb. For a given value of Req, a realistic taper gives a different 
value for Ref(0) than the square well for two reasons. 
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Table 2. Experimental and derived quantities at zero degrees 


1(0) (10-*4 cm?/ster) 


Element Li Be C N Ovi VAL Cu | Ag | Cd Pb 


84 Mev 1-2 bey] 5-1 | 15-5 | 29 63 
136 Mev 0-52| 0-76} 1-:19]1-72 | 2-19 | 4-15} 14-6 33:0 | 73-1 
290 Mev 1-25 DD 130 


Element i C N Cu 


84 Mev 0-508 | 0-878 ; 13-7 
136 Mev yy 0-336 | 0-585 | 0-768 | 0-96 | 2- 10-22 
290 Mev 0-827 “2 14-50 


Klement 


84 Mev 
136 Mev 
290 Mev 


Firstly J Vr?dr is smaller ; secondly the edges of the nucleus are more 
transparent and consequently the average value of F is more nearly equal 
to unity. The two effects tend to cancel. 

Accordingly Ref(0) has been calculated for the tapered potentials used 
before using the values of Rk previously determined and the experimental 
value has been used to assign a value to V. This is tabulated in table 3. 
Also tabulated is the value of V that is calculated in Born approximation 
(showing therefore the effect of nuclear absorption). The results are 
in satisfactory agreement with a reasonably low value of V as required by 
Taylor's fit to total cross sections. In table 3 the errors assigned are those 
of the experiment only ; in addition are errors in R and in the other para- 
meters. These yield an overall error of about 3 Mev in the values of V. 


§ 9. Discussion 


The values of the nuclear radius and of the potential derived in this 
paper are derived for each nucleus individually. By contrast, the values 
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of the radius or potential derived by Taylor or by Shapiro are derived 
assuming that the potential is the same for all nuclei and that R varies as 
a+bA™. It is gratifying that such good agreement is obtained, once an 
appropriate tapered potential has been chosen. This paper has therefore 
verified that the optical model parameters are nearly constant for all 
nuclei. The most accurate radius parameter determined is Req. This 
has been plotted in fig. 1 versus 41/3, R,, and Reg for the charge density 
distribution have been accurately determined by Hofstadter (1956) and 
are Ry=1-:07<x Al? fermis and Reg.=1-:18x A? fermis for lead and 
1-35 x A fermis for light elements. The difference in the ratio between 


Hofstadter’s results and those assumed here are due to slight differences 
in the taper assumed. 


Table 3. Nuclear Potential (mev) deduced in this paper 


Vy 35 


Born Approx. V iss Vo90 Veg | Vis from op 

Li 14-2 18-441-5 

Be 14-2 19-9+1-3 28+3 

C 13-8 18-141-5 |14+4 | 2643 

N 14-3 19-:9+1-3 

O 14:3 20-1+1-4 

Al 10-1 15-4+1-7 }15+3 | 20+4 

Cu 6-6 12:-3+1-5 | 6+6 13-:2+1-0 
Cd 4-9 11-242-3 

Pb 3°6 14 45 |234+15 14-0+1-0 

Taylor 13-5 10 20 


This plot of the nuclear radius shows clearly that the radius of the 
potential is larger than the radius of the charge distribution. The amount 
of this difference varies from 0-5 fermis for low atomic numbers to 1-25 
fermis for high atomic numbers. The difference for low atomic numbers 
must be attributed to the effect of nuclear forces (Drell 1955). However 
if the radius for lithium were multiplied by 0-86, to allow for a variation 
of V with angle (§2), then the 0-5 fermi difference between the potential 
well and the charge density radii would vanish. The effect of nuclear 
forces should be about the same for high atomic numbers; a difference 
between the radius of the charge distribution and the radius of a matter 

‘distribution, if real, would give an additional effect which can be as large 
as 0-7 fermis. 

The plot of nuclear radius versus A'* seems quite smooth in the region 
of carbon. The high binding energy of carbon seems to have less than 2% 
effect on the nuclear radius though there is an appreciable effect in the 
apparent radius calculated from deuteron stripping cross sections (Holt 
and Marsham 1953). 
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$10. COMPARISON WITH ABSORPTION Cross SECTIONS 


The radii obtained at 140 mev from absorption cross sections are some- 
what dependent on the value of «, the value assumed for the effect of the 
exclusion principle upon nucleon-nucleon scattering inside nuclei. 
Nevertheless, it is interesting to compare the results of Voss and Wilson 
(1956) for the nuclear radius, for the heavy elements. 

Table 4 compares the nuclear radii Reg deduced from the angular 
distributions in this paper, with the values Ry and R, deduced from 
neutron and proton absorption cross sections and also Renae the 


Table 4. Comparison of Radius determinations units of L0~8 cm 


(fermis) 
Pb Cd Cu 
Req 8-5 +0-2 | 6:92+0-2 6-03 40-15 
Ry 7:65+0:] 6-16+0-1 5:35 +0-1 
Bp 8:3 +0-4 | 7:-23+0-45 | 5-61+0-4 
core 6-99 5-69 4-68 


charge distribution radius. The absorption cross sections only give 
reasonable values for heavy elements. The same taper has been assumed 
as in this paper. 

It is seen that the values of R deduced fron neutron absorption 
experiments lie between the value deduced from neutron scattering and 
electron scattering. This also is consistent with the idea that the neutron 
radius is larger than the proton radius. Alternatively it might suggest 
that the real part of the nuclear potential has a larger radius than the 
imaginary part, the latter being fairly close to the charge radius. 

Note added in proof.— Recent measurements of the ratio of absorption 
of 700 Mev positive and negative 7-mesons by Cool show that the neutron 
and proton radii are essentially the same. 

It should be noted that although these radii are in agreement with the 
radii deduced from 84 and 290 Mev neutron data (see § 3), radii deduced 
from proton scattering depend upon the energy, and at higher energies 
are smaller than those deduced here. This fact, noted before by Shapiro 
(1955), may indicate a dependence of radius upon impact parameter. 


§ 11. Discussion OF THE POTENTIAL 
The value of the nuclear potential derived from this experiment is high 
compared with that obtained from total cross section measurements. It 
is however, more definite in the present analysis ; attempts have been 
made to fit 260 Mev and 400 Mev neutron total cross sections by assuming 
that V=0 at these energies (Nedzel 1954). By this means a linear plot of 
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f versus A’3 was obtained. However if V=-0 then Ref(0)=0 and it is 
then impossible to fit the neutron scattering data of Ball. 

The discrepancy between the values of V found from this paper and 
those of Taylor may lie in two of Taylor’s assumptions. Firstly the total 
cross sections for light nuclei are not very sensitive to the nuclear potential, 
and the value of V was fitted to the heavy nuclei. It is possible that V 
is slightly different for light and heavy nuclei. Secondly, part of the 
discrepancy is removed by considering the taper of the nucleus. 

There remains the effect of the exchange of particles upon the cross 
section theorem. In the examples of the breakdown of the theorem (in 
its elementary form) given by Mapleton (1954) the forward scattering 
amplitude is larger than that giver by kop/47. This would seem to be a 
general result. 

An order of magnitude estimate of its failure can be given by the 
following indirect argument. The production of neutrons by proton 
bombardment of nuclei at high energies has been interpreted as a quasi 
pn scattering with a neutron inside the nucleus. The energy distribution. 
of the outgoing neutrons is in qualitative agreement with a reasonable 
momentum distribution of the bound neutrons (Mandl and Skyrme 1952). 
A similar effect should occur for quasi n—n scattering inside the nucleus, 
and the residual nucleus should be in one of the lower lying levels according 
to the momentum of the struck nucleon. This inelastic scattering of 
protons has been observed. Some of the neutrons will be elastic, and the 
extent to which these are exchange neutrons will govern the failure of the 
cross section theorem. 

By a plausible extrapolation of Strauch’s data, we find for carbon that 
inelastically scattered neutrons leaving the nucleus in the 4:4 Mev excited 
state have an intensity about 4°, of the elastic neutrons. We might 
tentatively estimate the contribution of exchange neutrons to the forward 
elastic scattering intensity to be about half this or 2%. The effect on 
the scattering amplitude will then be about 10%, with a corresponding 
effect upon the cross section theorem. ; 

Accordingly, it is possible that the values of V deduced in this paper 
are high by about 10°—which could eliminate the discrepancy between 
the results of this paper and those of Taylor. 

This conclusion is supported by the calculations of Chew for neutron— 
deuteron scattering (Chew 1948). Our neglect of exchange corresponds 
to his most elementary treatment. Inclusion of exchange gives a larger 
cross section at small angles. The exchange effects in the deuteron also 
give a narrower diffraction pattern, so that the value of the radius, as 
determined in this paper, may have to be modified. — 


§ 12. COMPARISON wiITH ToTaL Cross SECTIONS 
It is not certain that the radii and potential shapes of the real and 
imaginary parts of the potential are the same, and indeed the comparison 
in table 4 suggests that they are not. In view of this, an extensive 
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comparison with experiment of total cross sections calculated from the 
parameters here deduced seems unjustified. For light nuclei also, the 
shape of the potential is in considerable doubt. 

Assuming the inelastic cross sections measured by Voss and Wilson 
(1956) and the radii here presented and using also the taper assumed, the 
real part of the nuclear potential is deduced for copper and lead. The 
results are included in table 3. For carbon the error in this procedure is 
too great to be meaningful. The total cross section has been calculated 
from ail the parameters in this paper and agrees with experiment to 5% 
within an error of 15%. These comparisons show that there is nothing 
drastically inconsistent in the procedure. 
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SUMMARY 

Angular distributions and energy spectra of the neutrons from the 
reaction 1B (d,n)!1 C have been measured at 1-08 Mev bombarding 
energy using the photographic plate technique. Energy levels were 
found in C at 2-02+0-04, 4-24+-0-03 and 4-73+-0-03 Mev. A large part 
of the reaction cross section must be due to compound nucleus formation. 
The experimental angular distributions can be adequately represented if a 
pure stripping curve is added to an isotropic distribution from compound 
nucleus formation. The orbital angular momentum of the captured 
protons has been determined. 


§ 1. INTRODUCTION 


PREVIOUS measurements have been made in this laboratory (Graue 1955) 
on the angular distribution of the ground state neutron group from 
0B (d,n) 14C. The distribution obtained at 0-85 Mev bombarding energy 
was found to be almost isotropic, thus indicating a large contribution 
from compound nucleus formation. Paris and Hndt (1954), however, 
investigated the same reaction at 0-58 mev deuteron energy, and found a 
pronounced stripping contribution to the reaction cross section even at 
this low energy. 

The purpose of the present investigation is to determine the relative 
importance of stripping and compound nucleus formation at 1:08 Mev 
bombarding energy and to determine the angular momentum of the 
protons captured by the stripping process. 


§ 2, EXPERIMENTAL CONDITIONS 


The deuterons were accelerated in the Van de Graaff generator at the 
University of Bergen. A 1B target supplied from the Atomic Energy 
Research Establishment, Harwell, England was used. The target 
thickness was approximately 80 kev for 1-08 Mev deuterons. ford C2 
Nuclear Research Plates with emulsion thickness 200 were used. The 
neutrons entered the plates at grazing incidence. They were arranged in 
such a way that they covered the range of the angle 0°<0@< 160", 6 being 
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the angle between deuteron and neutron directions. The distance between 
the target and the scanned areas was 120-130 mm, and independent of 6. 
The total deuteron charge on the target was 35 A-hours, which gave a 
suitable density of recoil proton tracks in the emulsions. An improved 
photographic method described by Gibson (1953) was used, thus 
simplifying the measurements and calculations. 


§ 3. NEUTRON SPECTRA 


Neutron energy spectra have been measured at the following angles : 
0°, 15°, 30°, 53°, 75°, 90°, 105° and 135°. Figure 1 shows the 30°-spectrum 
after correction for the variation with energy of the neutron-proton 


Fig. 1 
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scattering cross section and for the escape of protons from the emulsion. 
The errors indicated are statistical errors only. The four neutron groups 
correspond to levels in 14C at 0, 2-02+-0-04, 4-24+.0-03 and 4:-73-40-03 Mey. 
The first excited level seems to be a little higher than 1-90 mev, as pre- 
viously reported (Ajzenberg and Lauritsen 1955) 


§ 4. ANGULAR DISTRIBUTION OF THE NEUTRONS 


Burke et al. (1954) studied the angular distribution of neutrons from the 
0B (d,n) 'C reaction at 1-06 Mev bombarding energy, but could not 
differentiate between the various neutron groups. It was, therefore, 
desirable to measure the angular distribution anew, using neutron 
detectors with energy discrimination. 

The angular distributions obtained are shown in fig. 2. The distribu- 
tion of the ground state neutrons in fig. 2 (a) is almost isotropic, thereby 
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indicating a large contribution from compound nucleus formation. This 
distribution is similar to that obtained at 0-85 mev bombarding energy 
(Graue 1955). Figure 2 (b) shows the angular distribution of the neutrons 
corresponding to the 2-02 mev level in "C. In spite of the poor statistics 


Fig. 2 
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it can be said that the compound nucleus formation is responsible for a 
substantial part of the reaction yield. Figure 2 (c) and (d) shows the 
distributions of neutrons corresponding to levels in "C at 4:24 and 
4-73 Mev, respectively. The experimental points can be adequately 
represented if a pure stripping curve is added to an isotropic distribution 
from compound nucleus formation. The theoretical stripping curves 
(stippled) were calculated by means of a formula from the non-coulomb 
stripping theory (Friedman and Tobocman 1953). The formula has been 
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rearranged in a form convenient for numerical calculations (Enge and 
Graue 1955). The nuclear radius of the target nucleus used in the 
calculation was R=0-43x10-12cm, which was computed from the 
expression R—=(0-17-+0-12241%)+10-% em. The J-values found by 
comparison between the experimental angular distributions and the theore- 
tical curves were 0 and 2 for the 4:24 Mev and 4-73 Mev level, respectively. 
These values are in agreement with previous measurements on the mirror- 
reaction 1°B (d, p) 14B (Paris et al. 1954). 

The large contribution from compound nucleus formation in all the 
angular distributions may be explained by the fact that our exposure is 
made at a deuteron energy where there is some indication of a very broad | 
resonance (Ajzenberg and Lauritsen 1955). 


ACKNOWLEDGMENTS 


The author is indebted to Professor B. Trumpy for his interest in this 
investigation, to Mr. A. Otfedal and Mr. 8. Hauge for much of the micro- 
scope work and to Det Videnskapelige Forskningsfond ay 1919 for financial 


support. 


REFERENCES 


AJZENBERG, F., and Lauritsen, T., 1955, Rev. Mod. Phys., 27, 77. 

Burke, W. H., Risser, J. R., and Paruies, G. C., 1954, Phys. Rev., 93, 188. 
Ener, H. A., and Grave, A., 1955, Universitetet i Bergen, Arbok 1955, nr. 13. 
Frrepman, F. L., and Topocman, W., 1953, Phys. Rev., 92, 93. 

Gipson, W. M., 1953, Phil. Mag., 44, 297. - 

GRAUE, A., 1955, Universitetet i Bergen, Arbok 1955, nr. 1. 

Paris, C. H., Vauckx, F. P. G., and Enpt, P. M., 1954, Physica, 20, 573. 
Paris, C. H., and Enpt, P. M., 1954, Physica, 20, 585. 


fe lOs tem 


CVI. Isotopic Spin Selection Rules 
VIII: Charge Independence and the Comparison of Isobaric Triplets 


By D. H. Winxryson 


Cavendish Laboratory, Cambridget+ 
[Received December 20, 1955] 


ABSTRACT 


The experimentally-determined positions of the first 7’=! states of self- 
conjugate nuclei up to “°Ca are compared with their positions calculated 
from the known ground state mass differences of the 7’,=-+-1 and 7,=0 
isobars. The Coulomb correction is made in a completely automatic 
manner that does not involve an assumed nuclear radius but uses the 
experimental energy difference in the corresponding mirror pair. Close 
correspondence between calculated and experimental positions is found, 
the former being on the average higher by 35 kev. If this difference is 
taken to reflect a departure from equality of the n-n and n—p singlet 
forces the latter is stronger by about 1°/ (square well). This is to be 
compared with the difference of 3°, (square well) in the same sense found 
between p—p and n—p and n—p singlet forces from the zero energy scattering 
lengths. 


$1. InrRopUCTION 
In the series of which this paper is a part we have been chiefly concerned 
with testing the rule that if electric dipole radiation is to be freely emitted 
by a self-conjugate nucleus then the isotopic spin must change by one 
unit. The operating of the isotopic spin selection rules has in fact only 
been tested in such experiments on self-conjugate nuclei and in experi- 
ments in which 7'=0 particles provoke reactions and are produced as a 
result of reactions in self-conjugate nuclei. As is well known the isotopic 
spin selection rules which demand full charge independence (n—-n=n-p= 
p-p) for their formulation, can for all practical purposes be replaced, 
when we deal with self-conjugate nuclei, by the charge parity selection 
rules (Kroll and Foldy 1952) which are based on the much less restrictive 
requirement of charge symmetry (n-n=p-—p). So although we continue 
to talk of isotopic spin selection rules because we should like charge 
independence to be true we must not lose sight of the fact that all experi- 
ments so far performed, with one or two poorly-founded exceptions, are 
strictly relevant only for establishing charge parity selection rules and. 
charge symmetry in which we are much more willing to believe a priort 
anyway. Our understanding of the situation in self-conjugate nuclei 
is now fairly good both experimentally and theoretically and, for example, 
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paper VII of this series illustrates the general way in which the operation 
of the selection rules changes over a wide range of energy of excitation in 
a single nucleus (160). Charge symmetry, is then well established 
experimentally on the basis of the selection rules and on the basis of the 
well-known comparison of the level schemes of mirror nuclei. It is also 
very strongly to be expected on rather general theoretical grounds 
connected with the supposed origin of nuclear, forces and so it is now 
appropriate to ask how we may extend our investigations to full charge 
independence which is comparatively poorly established. experimentally 
and which cannot be expected theoretically with the confidence that we 
feel for charge symmetry. 

There are three chief ways in which we may approach this problem of 
establishing charge independence from a study of light nuclei : 

(i) We may study what we believe to be corresponding isobaric states 
in even nuclei of different 7’, and see to what degree their dynamical 
properties are related in the way required by the isotopic spin being a 
good quantum number. Alternatively we may study the competing 
modes of decay of an excited system when the branching ratios for the 
emission of different particles to different members of an isobaric multiplet 
are related by, among other things, vector coupling coefficients in the 
isotopic spin. Both these methods suffer from the disadvantage that 
before comparison of the interesting quantities (reduced widths) may be 
made we must eliminate extraneous factors such as transmission 
coefficients. This is usually a very uncertain process, often depending 
rather sharply on what we take as the radius of the system. And often we 
donot know how to weight different orbital angular momenta of the product 
particles in the second type of investigation. Because of such considera- 
tions and because of the fact that the systems usually pass through rather 
highly excited states where the Coulomb perturbation may effectively 
relax the isotopic spin it seems unlikely that such measurements, although 
they are very interesting, will ever give quantitative information about 
the degree of charge independence. It is more probable that they will 
be useful as a means for studying the interplay of the Coulomb perturba- 
tion and the finite lifetime of the intermediate system in determing the 
‘ effective ’ purity of isotopic spin states. 

(ii) We may study the properties of light nuclei of odd mass that we 
may make by bombarding a 7',=+1 target nucleus with protons, so 
forming both 7’=} and 7'=}$ states. Full isotopic spin rather than charge ; 
parity selection rules are then required to discuss the decay of these two 
classes of states. The two systems that suggest themselves are 4C--p 
and *O--p ; already there is some suggestion that the first 7 =? state of 
1°N, characterized by a small reduced width for neutron emission to the 
ground state of N, may have been found in the first of these systems 
(Bartholomew, Brown, Gove, Litherland and Paul 1955). This type of 
experiment must certainly be pressed hard for it offers rather good hope 
of seeing isotopic spin selection rules proper at work but, like the experi- 
ments in class (i), it suffers from passing through regions of rather high 
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excitation where, even if specifically nuclear charge independence were 
complete, the Coulomb forces might blur the definition of the isotopic 
spin beyond meaningfulness. The preliminary work on the 180 +p 
system for example (unpublished work at the University of Minnesota and 
the Rice Institute) contains no sign of any isotopic spin selection rules at 
work ; here the Coulomb perturbation is expected to be more severe than 
in the former system. 

The suggestion of Peaslee and Telegdi (1953) that the resonances of 
(y.n) and (y,t) processes be compared is a similar one but great technical 
difficulty is involved in the experiments. 

(iii) The third possibility is to study rather more closely than has been 
done so far the evidence contained in the relative energies of the isobaric 
multiplets themselves. But for the Coulomb energy and the neutron— 
proton mass difference the members of, say, a 7’=1 triplet should have the 
same mass if charge independence holds and we neglect small effects. 
This would not be so if we had merely charge symmetry. It has been 
obvious for some time (e.g. Sherr, Muether and White 1949, Lauritsen 
1952, Staehlin 1953, Inglis 1953, Ajzenberg and Lauritsen 1954, Burcham 
1952, 1955) that there is at any rate approximate equivalence of mass 
in some cases but the very great uncertainties involved in extracting the 
Coulomb energy have hindered a serious examination of charge indepen- 
dence from this point of view. We shall now see how, if we are prepared 
to begin by accepting charge symmetry as complete, we may make quite 
a strong statement about charge independence based on the relative 
energies of the isobaric multiplets in nuclei of even mass. 


§2. THE CouLtomB ENERGY DIFFERENCES 


There are several reasons why it is difficult actually to compute the 
Coulomb energy : | 

(i) The expression for the Coulomb energy contains the ‘radius’ of 
the nucleus in the denominator. It is not very clear what ‘radius ’ 
means and how radii determined in different contexts are related ; they 
differ by 50% or more for the same nucleus. 

(ii) The usual evaluation of the Coulomb energy neglects exchange 
effects which is certainly incorrect. It is not clear how they should most 


properly be allowed for (see e.g. Cooper and Henley 1953, Carlson and 


Talmi 1954, Peaslee 1954) but they can change the deduced radius by 
20% or more. 

(iii) The usual assumption that approximates the nucleus by a uniform 
sphere of charge is incorrect both from the point of view of the shell 
model and of direct experience with the scattering of fast electrons which 
shows that light nuclei are largely ‘ edge ’ with only a small region of 
constant charge density if any at all (e.g. Fregeau and Hofstadter 1955). 
Owing to the mutual repulsion of the protons the form of the distribution 
will depend on Z for a given A. Sieacen 

(iv) Even if we knew the effective radial charge distribution we should 
still not know sufficient details of the wave functions to calculate the 
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Coulomb energy accurately. Nuclear forces favour the pairing of spins 
and so nuclei with even numbers of protons tend to have relatively greater 
Coulomb energies than those with odd Z. This was remarked on by 
Feenberg and Goerzel (1946) and shown up very clearly by Carlson and 
Talmi (1954). The effect is particularly striking in the p-shell. 

(v) There is a complicated interaction between Coulomb and nuclear » 
forces in which the former tend to inercase the radius of the nucleus 
thereby changing the binding energy due to the latter. This effect will 
differ between members of an isobaric triplet because of the differing 
nuclear charges and so there will be a ‘nuclear’ contribution to the 
‘Coulomb’ energy differences which is not taken into account in the ~ 
conventional calculation which uses a constant radius parameter * ry’. 

Because of these various effects and uncertainties it is clear that we 
cannot accurately calculate the required differences between the Coulomb 
energies of the members of the isobaric triplets. We may find an empirical 
recipe in terms of parameters for the conventional * uniform sphere * 
calculation but even apart from the fundamentally unsatisfactory nature 
of this approach we could not look to it for better accuracy than several 
hundred kilovolts because of the finer effects we have listed. 

If we are now prepared to assume the charge symmetry of nuclear 
forces we can do much better than this. Suppose we wish to discuss the 
lowest T=1 triplet T,=+1, 0, —1 (Z—1, Z, Z+1) of mass A (=2Z). 
Z—1and Z of mass A —1 are now mirror nuclei and their mass difference, 
corrected for the n—p mass difference, is due to the difference in “ Coulomb * 
energies of the two nuclei where in ‘ Coulomb ’ we include effects such as 
(v) above. We now multiply this energy by (1—1/A)"* to allow for the 
change in radius on adding another particle and suppose that it repre- 
sents the ‘Coulomb ’ energy difference between the 7',=+1 and 7',=0 
members of the triplet. This procedure eliminates most of the difficulties 
mentioned above and has the advantage of being completely automatic, 
with no parameters to be manipulated to secure ‘ best fit ’. The question 
of the radius and its definition are never raised. A particular advantage 
is that we may hope that we are allowing for the pairing effects in Z 
since the even—-odd character of the charges is the same in the even A 
nuclei in which we are interested and in the mirror nuclei that we use for 
comparison. We are also automatically allowing for the various effects 
listed above which amount to the statement that ‘7 )’ depends on Z 
for a fixed A. The most serious defect of the method is probably that it 
makes no allowance for the difference in symmetry properties of the 
protons in the A and A—1 systems. We must hope that the influence of 
these differences is small. 

Further lack of confidence may come from a comparison of the level 
structures of mirror nuclei where levels get out of step in the two members 
by as much as a few hundred kilovolts at excitations of a few Mev. To 
some degree this may be seen to be associated with the Thomas shift 
(see later) but it is also doubtless due in part to the Coulomb energy 
depending on the level we deal with. When the states concerned differ 
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considerably in structure from the ground states we must anticipate that 
they have different Coulomb energies from the ground states. In this 
connection our present procedure is not so bad as it looks since the 
lowest 7'=1 state of a T,=0 nucleus is closely related to the ground state 
even though the 7=0 levels at the same excitation as it may not be. 
Nor are the levels at the same exictation in the mirror A—1 nuclei neces- 
sarily closely related to their own ground states so the relative displace- 
ments in the A—1 nuclei that we have commented on as discouraging are 
not necessarily relevant. 
A similar method to the present one has been used by Inglis (1953). 


§3. EXPERIMENTAL RESULTS 


Our procedure is now to extract from experiment the energy difference 
between the mirror A—1 nuclei and use it, corrected by the small factor 
(1—1/A)%, together with further experimental knowledge about the mass 
difference between the ground state of the 7,—+1 member of the A 
triplet and that of the T7,—0 member to predict the energy above the 
ground state of the 7’,—0 nucleus at which we expect to find the first 7=1 
state if we have full charge independence. Call this energy Ep_, (calc.). 
From experiment we now find the actual excitation energy Ep_, (exp.) 
of this state and hence the difference 4=H,_, (calc.)—Hp_, (exp.) 
between calculated and experimental positions. If no other factors 
entered into consideration 4 would now be a measure of the departure 
of nuclear forces from charge independence. 

The two tables present the material available as far as *°Ca. 

Table 1 lists the energy differences between members of the mass A—1 
mirror pair (allowing for the n—p mass difference). We also list the energy 
difference between the ground states of the 7’,=-+1 and 7',—0 members 
of the mass A triplet (allowing for the n—p mass difference); here a 
positive difference means that the 7',=-+ 1 member is the heavier. For the 
cases of A =26 and 38 we do not need this latter information and although 
it may be known we do not list it ; the measurement of energy difference 
in the mass 26 and 38 nuclei refers directly to the 7'=1 states of 7°Al 
and 38K as is signified by the subscripts. For masses 10 and 14 it is also 
useful to know the 7,——1 to T,=0 mass differeuce so we list this as 
well. 

Table 2 lists as H,_, (calc.) the calculated excitation energy of the first 
7'=1 state in 7',=0 using the data of table 1 and the automatic method 
described above. If the energy of the 7,——1 member of the triplet 
were accurately known we could calculate the expected position of the 
7'—1 state in 7',=0 without recourse to the A—1 murror pair, basing 
ourselves on the assumption of charge symmetry only. It is in fact 
known in several cases but usually only very inaccurately. In the two 
cases, namely A=10 and 14, the energy differences are well enough known 
from the 7',=—1 side (°C and QO) to use this more satisfactory pro- 
cedure. In just these two cases Hp_, (calc.) refers to the energy computed 
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Table 1. Energy differences, allowing for n—p mass difference, between 


mirror nuclei of mass A=1 and between the ground states of the 
T.=-+1 (or —1) and 7,=0 nuclei of mass A (with the exception 
of 2°Al and *8K—see text). 


Mirror 
pair 


Energy 
difference 
(Mev) 


Ref. 


Even 
mass 
nuclei 


Energy 
difference 
(Mev) 


5He—*Li 
"Li—7Be 


9Be—°B 
11pB 110 
13(1__ 13) 


J SN ame: 50 
170 a7) 17 
19} — 19Ne 
21Ne—21Na 
23Na—°Mg 
25Mo —25A] 
27 A] 27S} 
29G} __29p 
31p 819 
339) — 830] 
351] 35 


0-89 +0-20 
1-645 0-001 


1-852 +.0-002 
2-762 +.0-003 
3-003 +.0-002 


3°487 +0-005 
3-549 +.0-006 
4-027 +0-008 
4:30 +0-08 

4-841 +0-010 
5-084 +.0-025 
5-584 +0-010 
5-749 +0-010 


‘He—*Li 
SLi—®Be 
1 Be pa} oR 
100) a3 10B 
2B — 12¢3 
1a(Se 14N 
se T@) pO) 
16N ce 169 
180) ==18h) 


208 —20Ne 
2Ne—2Na 
24N a, —4Mo 
26Mg =e SAlr_y 


28 A] —28Qj 
309) —30p 
32P 326 
34S __34Q)] 

86] = 386A 


2:77 +0-05 
15-20 +0-04 
—0-226 +0-003 

4-42 +0-06 
12-588 +0-025 
—0-627 +0-001 

5-952 +0-008 

9-62 +0-03 
—2-453 +0-009 

6-261 +0-012 
—3-623 +0-006 
4-730 + 0-005 
—5-006 +0-010 
3-868 +0-011 
—5-08 +0-04 
0-925 +0-004 
—6-30 +0-03 


—0-068 + 0-005 


— 
7 


Sel 6 Come 


387A 37K 


39K —39Ca, 


—6-37 +0-13 
0-543 40-015 


8A 8K y_, 
40K A059, 


(a) This figure is the mean of two. One derives from the observed *He—*Li mass 
difference—see Ajzenberg and Lauritsen (ref. (j) below). The other comes from 
assuming that the even—odd alternations in 7, observed by Carlson and Talmi (1954) 
between A=7 and 9, and A=11 and 13 hold also between A=8 and 5 which enables 
us to compute the energy difference according to their recipe. 

(b) C.W. Li, W. Wuatine, W. A. Fowter, and C. C. Lauritsen, 1951, Phys. Rev., 83, 512. 

(c) J.D. Kiyeron, J. K. Barr, H. O. Conn, and H. B. Wizxarp, 1955, Phys. Rev., 99, 1393. 

(d) C. W. Li, 1952, Phys. Rev., 88, 1038. 

(e) H. Roprriox, and C. Wone, 1953, Phys. Rev., 92, 204. 

) F.J. Boney, and D, J. Zarrarano, 1951, Phys. Rev., 84, 1059 who measured *18(f*)31P, 
also R. N. H. Hastam, R. G. Summers-Gir1, and E. H. Crossy, 1952, Canad. J. Phys., 
30, 257 who measured the threshold for **S(y,n)?48 from which, together with the 
828 and *1P mass defects, the same datum may be deduced. 

(g) No reliable data exist here. This energy difference is computed from the semi-empirical 
formula for Peaslee (1954) for mirror nuclei and the probable error taken from his 
estimate of its reliability. 

(h) The old concordant but unreliable cloud chamber data which have been shown to be 
systematically too low have been discarded and this estimate follows (g). 

(i) KR. G. Summers-Girt, R. N. H. Hasiam, and L. Karz, 1953, Canad. J. Phys., 31, 70 
measured the threshold for 4°Ca(y, n)?®Ca which combines with the mass defects of 
“Ca and *°K to give the required datum. 

(j) F. Aszenserc, and T. Laurrrsen, 1955, Rev. Mod. Phys., 27, 77. 

(k) D. E. AtBurcER, 1952, Phys. Rev., 88, 1257. 

C. Wona, 1954, Phys. Rev., 95, 765. 
(1) K. Stmenann, 1946, Phys. Rev., 70, 127. 
A. Heperan, and D. Linp., 1952, Ark. Fys., 5, 177. 
(m) ae en pe PATTER, hae W. BurcuHNER, 1952, Phys. Rev., 87, 51. 
poe NDEVILLE, ©. P. Swan, 8. D. CHarrerRJEn d D. M. 
Phys. Rev., 85, 193. is ol 
(n) D. Green, and J. R. Ricuarpson, 1954, Phys. Rev., 96, 858. 
(o) L. Fe~pmay, and C. 8, Wu, 1952, Phys. Rev., 87, 1091. 
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Table 2. Calculated and Experimental Excitations of the First 7=1 


state of Self-conjugate Nuclei and the Differences AH 2y)( Caley) 
—Hp_, (exp.). 


Nucleus Er_, (cale.) Ep_, (exp.) A Ref 
(Mev) (Mev) (Mev) Sag 
3°6] +0-2 3°57 +0-0] 0-04 +0-2 (q) 
16:77 +0-04 16-72 +0-01 0:05 +0-04 (r) 
1-775 +0-027 1-739 +0-003 0-036 +0-027 (q) 
15-27 +0-03 15-10 +0-02 0-17 40-04 (s) 
2-336 +0-004 2-313 0-004 0-023 +0-006 (q) 
13-08 +0-03 12:95 +0-01 0-13 0-03 (t) 
1-029 + 0-009 1-075 +0-010 —(0-046 +0-014 (u) 
10-221 +0-015 unknown — = 
0-61 +0-03 near near 
0-592 +0-003 0:02 +0-03 (v) 
9-504 +0-012 9-47 +0-10 0-03 +-0-10 (w) 
5:018 +0-025—H| 5-006-+-0-010—Z 0-012 40-027 (x) 
9-384 +0-017 9-27 == 0-11 Q-11 +0-11 (w), (y) 
0-61 +0-04 0-688 +.0-007 —0:08 +0-04 (z) 
6-86 0-10 7:02 +0-05 —0-16 +0-12 (w) 
—0:02 +0-12 0 —0:02 +0-12 (a’) 
6:06 +0-06 unknown — — 
6°72) 220: 12=8 6:37 +0:13—H 0-35 +0-18 (b’) 
fo ee Or unknown — = 


(p) We use the measurement of Boley and Zaffarano (f) above which, in view of the 


(q) 
(r) 


(8) 


(t) 


results of P. Srarniin, 1953, Helv. Phys. Acta, 26, 691 we re-interpret as referring 
to the 7 =1 state of **K. This entry is therefore subject to some doubt, particularly 
since *“K to which it was originally supposed the measurements referred may well 
have a lifetime and positron end point close to those of *K7=1. 


First T=1 state in 7'z=0 well established—see reference (j). 


The evidence for the first 7’=1 state being around 17 Mev was given by M. GrLL-MANN 
and V. L. TetEapi, 1953, Phys. Rev., 91, 169, on the basis of its formation by alpha- 
particle emission from 12C following E1 absorption. F. K.Gowarp, and J. J. WILKINS, 
1955, Proc. Roy, Soc. A, 228, 376, propose the same assignment and find more than 
one 7’=1 level in this region, the lowest well-established one being at 16-8-+-0-2 Mev, 
our quoted figure comes from T. W. Bonner, and C. F. Coox, 1954, Phys. Rev., 96, 
122. The assignment 2+ for this level from the !°C(y, «)*Be work is consistent with 
the known properties of *Li (j). 

That this first 7 —=1 state is around 16 mev is suggested by C. A. Hsrao, and V. L. 
TrLEepi, 1953, Phys. Rev., 90, 494, on the basis of its frequent production following 
El absorption in 1%O(y, «)12C. Many workers (see (j)) have observed that a level at 
15-10 mev (this energy comes from C. F. Cook, J. Marion and T. W. Bonner (j)) emits 
y-rays to the ground state in successful competition with alpha-particles. This 
strongly suggests 7’=1 and the emission to the ground state is consistent with the 
1+ assignment usually given to B. 

This assignment is not sure. 1°N is 2— so we seek a T=1 2— state in *°O around 
13 Mev. Two 2— states are known, at 12-51 and 12-95 mev. Both have very large 
reduced proton widths for d-wave emission to the ground state of 1°N as we should 
expect if they were simple states such as the ground state of '°N is presumed to be. 
This suggests that they may be the 7’=0 and 7'=1 linkages of a dj, proton to the 
ground state of 15N just as we suppose the ground state of 16N to be the linkage of 
a d,/, neutron (neglecting configuration mixing in both cases) to the ground state of 
15N. The lower 1°O state has a reduced alpha-particle width for emission to the 
4:43 mev state of 120 of 80% of the single-particle value, the upper state has a 
corresponding width of 4%. So if we are to choose our 7'=1 state from these two 
we choose that at 12-95 mev. This choice is supported by the associated 4 value 
of 130 kev which would become much bigger than any other example if we chose 
the lower state as 7=1. We may note as a comment on the large alpha-particle 
width of our supposed first Z’=1 2— state that this 77=1 1— state also shows an 
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(w) 


(x) 


(y) 


(z) 


(a’) 
(b’) 


é - icle width of 5% of the single particle value. It seems that the high 
Pay : Mn region is Selecta the definition of isotopic spin rather 


density of 7=0 states in L « 
aides note in particular that if we are correct in supposing that the 12-51 


12-95 mev 2— levels are closely-related 7’=0 and 7’'=1 states we might expect 

ee to contaminate each other hee effectively—see IDR iEb, WILKINSON, and 
A. B. Ciece, 1956, Phil. Mag., 1, 91, part VI of this series. We must admit that 
this identification is as yet unsure; it is clear from the density o levels near the 
ground state of 16N that all 71 states in this neighbourhood in 18O have not yet 
turned up. 
Only this state is known within 800 kev of the expected energy. Its properties are 
consistent with the 0+ required for its being 7’=1 (P. C. Pricer, 1955, Proc. Phys. 
Soc. A, 68, 553). It is formed in *°Ne(d, «)!*F and this has always stood ageie” its 
being recognized as of 7=1. However the total excitation m the intermediate Na 
system in these latter experiments is only 18-5 Mev at which energy we may still 
expect some resonance structure to persist. Under these conditions we must be 
prepared to find severe breakdown of the isotopic spin selection rules (A. M. Lang, and 
R. G. Tuomas, 1956, Rev. Mod. Phys. (in press); D. H. Wirxrnson, Phil. Mag. 
(paper VII of this series) ; C. P. BRowNE (private communication)) and until another 
and more acceptable state is found in this immediate energy region it seems that we 
may tentatively look on the present one as of 7’=1. 


N. P. HeypEenpoure, and G. M. TemmMer, 1954, Phys. Rev., 94, 1252, have observed 
a gamma-ray of 592+3 kev which they supposed to come from the 0+7=1 state. 
However they also observed that its half-life against presumed M3 radiation to the 
3+ ground state was less than 10-*sec. about ten times the single-particle speed 
which makes the identification doubtful. S. VEGors, and P. AxEL (private communi- 
cation) have sought an isomeric state in **Na following **Na(y,n)**Na and are able 
with good probability to place an upper limit of 10sec on the lifetime. It therefore 
appears probable that the 592 kev state is not 0+ but perhaps 1+ with the true 
0+7'=1 state lying closely above it, being fed but with its M1 radiation being of 
low enough energy to have passed unnoticed in Heydenburg and Temmer’s work. 
Again C. P. Browne, and W. C. Coss, 1955, Bull. Amer. Phys. Soc., No. 3, p. 44, 
find the 592 kev state strongly excited in **Mg(d, «)??Na which may tell against its 
being of 7’=1. Here the arguments of (u) about the production of 7=1 states in 
such reactions again apply although we should only use them to excuse a 7'=1 state’s 
appearance rather than to lead ourselves to expect it (¢e.f. the similar reaction 
*8Si(d, «)?®Al where the 7’'=1 state is not produced—C. P. Browne, 1954, Phys. Rev., 
95, 860, see also R. W. Kavanacu, W. R. Mints, and R. Suerr, 1955, Phys. Rev., 
97, 248). We therefore consider it likely that the 7=1 state is at a slightly higher 
excitation than the 592 kev state but close enough for the cascade radiation to have 
been missed, 

N. W. Guass, and J. R. Ricuarpson, 1955, Phys. Rev., 98, 1251, have studied the 
positron emitters *4Al, ?8P and **Cl and find in each of them evidence for allowed 
favoured transitions such as are expected from them to the analogue 7’=1 level in 
the 7z=0 nucleus. The quoted level positions then come from the appropriate 
compounding of measured gamma-ray energies (averaged with another datum for 
28Si—-see (y)). 

The situation in ?°Al is now quite well understood (see M. J. Lausrrz, 1955, Proc. Phys. 
Soc. A, 68, 1033 for a summary). The experimental datum here of greatest accuracy 
is a measurement of the Q-value for ?*Mg(p, n)2*Al7=1 direct to the first 7=1 state 


of **Al (ref. (c) above) so we do not need to bother with determining the position of 


this level relative to the ground state of ?°Al. In the entry of table 2 E is the energy 
difference, corrected for the n—p mass difference, between the ground state of °6Mg 
and that of 2°AI. 


J. M. Catverr, A. A. Jarre, A. E. Liraernanp, and EF. E. Mastin, 1955, Proc. Phys. 
Soc. A, 68, 1008, have measured the stripping properties of a level in 288i at 
9-16--0-17 Mev seen in the reaction 27Al(d, n)28Si and find a close correspondence 
with those of the ground state of *8Al seen in 27Al(d, p)*sAl. This leads them to 
suggest this level as the first 77=1 state of *8Si. The energy is in good agreement 
with the 9-37--0-15 Mev deriving from the decay of *8P—see (w)— and we have 
combined these values. 


Again the only justification for regarding this level as of 7’—1 is that it is the only 
known level within 600 kev of the expected energy. P. M. Enpt, J. C. Kiuynr, 
and C. vAN DER Leun, 1954, Phys. Rev., 95, 580—see also reference (w). 

D. C, PEASLEE, 1953, Nuovo Cimento, 10, 1349. 


STAEHLIN (1953)—main references. The calculated entry, hence 4, is subject to 
some doubt—see reference (p). In this entry of table 2 # is the energy difference, 
corrected for the n—p mass difference, between the ground state of 88A and that of 3K. 
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using the data from 7,=-++1 and without reference to the data from the 
mirror pair of table 1. In making this computation we have followed the 
method of Carlson and Talmi (1954) and estimated the (relative) Coulomb 
corrections between 7’,=+1 and 7',=0 and between 7',=—1 and 1220 
using LS wave-functions ; these are the more valid approximation in 
light nuclei. It is now interesting to check these results for B and 14N 
against what we should have obtained using the method based on the 
mirror transition which we have used in all other cases. In the case of 
MN the agreement is good and we get 2-302-+0-003 mev for Ep_, (calc.) 
from the mirror pair as against 2-336-+0-004 mev from the better PT =e 
method, This is good confirmation of the validity of our general method 
of using the mirror pair. In the case of !°B we must expect a bad estimate 
from using the mirror pair because the proton-rich member, °B, is actually 
unstable against proton emission. This will cause its mass to be depressed 
by the Thomas (1952) shift representing the fact that the loose proton can 
spend much of its time at the edge of the nucleus where its Coulomb 
potential energy is low. So we expect the °B—*Be entry of table 1 to be 
too low if we use the table 1 data. In fact this is so and #7_, (calc.) 
by the mirror pair method would have been 1-56 Mev as against 1:78 Mev 
of the better method and the 1-74 Mey of experiment. It is fortunate 
that in this, the case where we would have least confidence in the 
mirror pair method (with the exception of 4—1=5) we have the better 
J =—=- 1 method to rely on: 

The experimental 7’=1 assignments are discussed in the table 2 
references. 


§ 4. Discussion 


The entry A of table 2 shows that when we make the Coulomb correction 
in this automatic fashion we do indeed find extremely close agreement 
between the calculated and experimental positions of the first 7’=1 
states of 7',—0. In most cases the value and probable error in 4 fit 
together so well that it is difficult to see that there is any significant 
discrepancy at all. Only in C and 160 do we seem to find substantial 
disagreement. These cases may be somewhat special in that the excita- 
tion energy involved in 7',=0 is considerable and we are in a region of 
quite closely spaced 7’'=0 states (not so for the otherwise-comparable 
8Be) whose density is increasing with increasing excitation. Under 
these circumstances we may expect considerable isotopic spin mixing with 
these T’=0 states, and so displacement of the 7’=1 level by the usual 
mechanism of configuration interaction. This will tend to be towards 
lower energies because of the greater abundance of mixing (and repelling) 
states above than below. This is to some degree confirmed by the observed 
high impurity of the 7=1 states in this region of *O—see reference (t) 
of the table 2 References—although we should comment that in contrast 
the purity in the corresponding case of #C appears to be high enough to 
allow gamma-ray emission to compete with alpha-particle decay— 
reference (s) of the table 2 References. In O a special factor may be 
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seen at work to produce a displacement in the observed direction. We 
take the Coulomb shift from 1=N—1O where in the latter body we have 
6 lp protons. However in the 7'=1 state of 16Q) we have one proton 
removed from the Ip to the Id shell. This means that it sits somewhat 
nearer the edge of the nucleus and so we have a smaller Coulomb shift 
between the 7’=1 state of 16O and the ground state of *N than we calcu- 
late from 4—15. A rough estimate of the quantitative importance of 
this effect, using shell-model wave functions for the last proton and a 
uniform sphere for the rest of the nucleus, gives 40 kev. However no 
such argument can be made for ?C. Another factor tending to produce 
a positive 4 for these highly-excited 7’=1 states is the Thomas shift. 
These states are either unstable against proton emission (1®O) or only 
just stable (*Be, 2C) and so we must anticipate a downward shift as we 
noticed in our discussion of 9B in §2. This is in the right sense but it 
appears to be substantially too small to account for the observed 4 in 
these cases. 

A significant tendency is for positive 4 values at any rate in the p-shell. 
There seems no general reason for this unless it be that the proton binding 
energy tends to be less for the 7’=1 state of 7,—0 than in the mirror 
nucleus with an attendant relative Thomas shift. This would have to be 
surprisingly large to account for the observed 4 values. 

There are two small effects which go in the wrong sense—tend to raise 
the 7=1 state of T,—0. As we noted in §1 part of the ‘Coulomb’ 
energy expresses the distortion of the nucleus by the Coulomb forces 
with attendant change in the binding by the specifically nuclear forces. 
This is allowed for by our use of data from the mirror pair. However in 
the even A nuclei the specifically-nuclear forces tend to be relatively 
stronger because of the extra neutron so we should expect the even A 
nuclei to be relatively smaller (smaller effective ‘7)’) and so have a bigger 
Coulomb shift. The second effect is closely related to this and is again 
that the even A nuclei will tend to be relatively smaller—because of the 
pairing effect of an even number of particles. An attendant effect to this 
last one goes sometimes one way sometimes the other and its failure to 
show up as an alternation in 4 may indicate that these fine effects are 
indeed as small as we should expect them to be. It is that when 7,=0 is 
even-even, 7’.=-++1 is odd—odd and vice versa. In the first case we 
should expect 7',—0 to be relatively smaller and so A to be decreased ; 
in the second case 4 should be increased by this effect. As we see from 
table 2, if anything even-even 7',=0 have the larger 4 values. 

Beyond the p-shell the probable errors are larger and no significant 
trend of 4 can be discerned. 

We may sum up this discussion by saying that, although there seems to 
be an unexplained tendency towards positive values, 4 is on the whole 
remarkably small. If we award a weight of 1 to 4 values with probable 
errors of 50 kev or less and divide the weight by 4 for each range of 
doubling of the probable error we find a mean A value for all nuclei 
listed in table 2 of +35 kev. It seems then that we can state that the 
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average discrepancy between calculated and expected positions of the 
7’=1 states in less than 50 kev. It is significantly positive. 


§5. CHarGE INDEPENDENCE 


Our 4 values are very small and so imply that the n—p bond of fi aW) 
is closely equivalent to the n-n bond of 7',=+-1. These nuclei differ 
only in the single bond so that if the n—p force is not the same as the n-n 
there will be a gross shift in both nuclear masses due to all the n—n (or 
p-p) and n—p bonds differing but the relative shift with which we are 
concerned here is due to the single bond alone. If we regard the departure 
from complete charge independence as a perturbation, presuming charge 
symmetry, we expect a relative shift of the order 


A~ | b* ASibdr 


where AS is the difference in strength between the n-n and n-p (chiefly 
318) interactions or 
AS 


where 4/8 is the fractional departure from charge independence, V is the 
depth of the 348 interactions (say 10 Mev for a square-well range of a 
=2-810~%cem) while P is the fraction of the time for which the two 
nucleons are within range of each other’s attraction (~(a/R)* for a big 
nucleus). If we insert our value of 35 kev for A for light nuclei with say 
P~4—4 we find that in the square-well approximation the n—p interaction 
is stronger than the n—n by 0-7—1°%. Allowing for the very rough nature 
of this estimate it still seems that these data on the positions of the first 
T=1 states cannot admit of more than a very few per cent departure 
from charge independence. 

It is interesting to compare this tentative conclusion that n—p forces 
may be stronger than n—n by about 1°% with the results deriving from 
n-p and p-p singlet scattering. The two singlet scattering lengths are 
very significantly different but, since both are very large, reflecting the 
very small ‘energy of the singlet deuteron ’, the difference in n—p and 
p—p well strengths is small. Again in the square well approximation the 
n-p strength is greater than the p-p by 3-3% after allowance for the 
Coulomb effect (Jackson and Blatt 1950) to be compared with our 1% 
for the excess of n—p force over n-n. The effects are seen to have the 
same sign and to be of the same order. Schwinger (1950) has pointed out 
that the difference between the n—p and p-—p forces may be semi-quantita- 
tively accounted for as due to the different magnetic interactions of the 
two systems, an explanation that would leave the specifically-nuclear 
charge independence complete. We should then expect just about the 
same difference between the n—p and n—n potentials as between the n—p 
and p-p. If we omit the relatively small orbit—orbit contribution to the 
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p-p magnetic potential the three magnetic contributions stand in the 
ratio : 


Vass ves es 
= bn” : ry Ta) : Pa(Hp—2) 
3°66 : 5-51 : —4-39 


The p’s are the magnetic moments in nuclear magnetons. It is not at 
all clear that this magnetic explanation for the difference between n—p 
and p-p scattering lengths is adequate since Schwinger found that a 
square well gave unsatisfactory results and a Yukawa-type potential 
was needed. In view of the current favouring of a repulsive core to the 
nucleon-nucleon potential it is unlikely that the satisfactory agreement 
of the Yukawa potential can be maintained (Salpeter 1953). But in 
any case our present result on the n—n versus n—p forces is consistent with 
the information from scattering lengths and the basic assumption of 
charge symmetry. It is also consistent with the much less accurate 
information relating to comparison of the n—p and p-—p effective ranges in 
the singlet state (see e.g. Storrs and Frisch 1954). 

If we are to invoke magnetic forces in our comparison of the n—n and 
n-p interactions from the isobar structure we must also consider the 
magnetic contribution to the mass of the isobars due for example to the 
coupling between the magnetic moments of the nucleons and those of the 
mirror nuclei to which they are added or, more properly, to the moments 
of the various parents in A—1. These effects are of comparable impor- 
tance with the nucleon-nucleon magnetic energies but vary in sign from 
case to case and so tend to cancel out in our evaluation of the mean J. 
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ABSTRACT 


The basic ideas of Brueckner’s self-consistent nuclear model are applied 
in a simplified, approximate form to the case of a finite nucleus. It is 
shown that it is possible to reconcile the observed values of surface energy 
and surface thickness, to explain the greater extension of the nuclear 
potential compared with the charge distribution, and also to fit the well- 
depth of the optical model for nucleon scattering. 


§ 1. INTRODUCTION 


A DIFFICULTY has recently been pointed out by Swiatecki (1955 b) in 
reconciling the value of the nuclear surface energy with the measurements 
of surface thickness obtained from electron scattering experiments. This 
problem is reinvestigated here in the light of the probable existence of a 
velocity-dependent mean potential, as suggseted by Brueckner ef al. 
(1954) and by Johnson and Teller (1955). 

In the self-consistent model of Brueckner, whose ideas are followed here, 
nuclear saturation is achieved by the reduction in the attractive force 
between pairs of nucleons as their relative velocity increases. In the 
idealized case of indefinitely extended nuclear matter, the effective 
potential between pairs of nucleons is completely described in terms of 
their relative momenta; in the case of a finite system the effective 
potential must depend upon both momenta and coordinates, and as has 
been pointed out by Bethe (1956), the mean central potential should 
then properly be taken as a potential matrix, diagonal neither in co- 
ordinate nor momentum space. 

To provide a simple model for the present discussion it has here been 
assumed that the effective interaction between pairs of nucleons may be 
described by delta-functions of position and differential operators (to 
express momentum dependence) ; the resultant mean potential is then a 
sum of products of functions of position with momentum operators, With 
these assumptions a relatively simple set of self-consistent equations may 
be written down for the individual particle states of the assumed. model 


wave-function. 


+ Communicated by Dr. B. H. Flowers. — ; 
{ The principal reference for this paper is Brueckner (1954 b). 
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The problem of determining self-consistent scattering parameters from 
an assumed inter-nucleon potential is not discussed here ; a simple form 
for the scattering amplitude is taken with parameters adjusted to fit 
some of the experimental data. 

An estimate of the nature of the solution of these self-consistent equa- 
tions is made by using the Thomas—Fermi approximation for the 
momentum distribution. This leads to certain definite results for the 
surface energy and for the shapes of the density distribution and of the 


mean potential well. 


§2. Tor Mopret Equations (INFINITE SYSTEM) 


In the case of infinitely extended nuclear matter the determinantal, 
model, wave function introduced by Brueckner is built out of individual 
plane-wave states, and Brueckner has shown how in good approximation 
the energy of the system may then be calculated as the expectation value 
of kinetic energy plus an effective potential proportional to the forward 
scattering amplitude between the plane wave states. In a complete 
calculation, self-consistent ab initio, this scattermg has to be calculated 
from the assumed two-body potential, and evaluated for particles inside 
nuclear matter, subject therefore to the mean potential field and to the 
exclusion principle. This full problem is not considered in the present 
paper ; only the form of the theory is used, slightly generalized to apply 
to a finite system, and with the constants in the scattering amplitude 
adjusted to give the correct volume saturation effects. 

In this model saturation is achieved through the velocity dependence of 
the scattering amplitude. For a determinantal wave function built out 
of plane-wave states, with momenta k; and normalized in a volume JV, 
the effective interaction energy between a pair of nucleons with momenta 
k; and k; is a function only of the relative momenta, 

(1/V)i(kj—ky). a eh ee 
In Born approximation, and neglecting the effects of exchange, ¢ is equal 
to the volume integral of the assumed two-body potential. In general the 
velocity-dependence of ¢ arises in part from the contributions of the 
exchange integrals, and in part from the energy dependence of the 
scattering amplitudes, when calculated exactly. 

On the average, for the momenta occurring within the nucleus, ¢ is 
negative, expressing the binding of the nucleons ; saturation will obtain 
ift becomes more positive, sufficiently rapidly, with increasing momentum. 
A simple form that shows this behaviour is 


t(k)=ty+t,hk? *.5 pated 1:4 ce 
with negative f, and positive t,. This might be regarded either as the 
first two terms in the complete expansion of ¢ in powers of k, or perhaps 
more realistically as an approximation describing the behaviour of ¢ 
over the important momentum range; in either case the parameters 
tg, t; will be functions of the density in the self consistent calculation. 
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Tt will be convenient, for the subsequent consideration of the finite 
system, to introduce symbols p and 7 to denote the particle and square of 
momentum densities; for plane waves in a volume V these are 


ped) eis ve (LIV) 2,h?: reall (3) 
Then the interaction energy density (of the infinite system) is 
(UV) 2c, WK —k;)=dtop?+ typ. ch st emcee C9) 


Here the cross terms in k;k; are assumed to vanish on summation.t 
Addition of the kinetic energy to this gives for the total energy density 
A g=(1/(2M)+t1p)7+htp? . . . . . . (5) 
where the coefficient of t may be interpreted in terms of an effective mass 
M* defined by 
LOMO A At pie Se ys * (8) 
In this infinite system where the states are plane waves with momenta 
fillmg up a sphere (in momentum space), the kinetic energy is a known 
function of the density and with four nucleons to each state. 
T= (s'O)p oo, a holt eS Bee eee 8) 
Then the energy per particle H=H/p is a known function of the density, 
and the equilibrium density and binding energy are obtained by 
- minimizing H(p). . 


§ 3. THe Mopet Equations (FINITE SYSTEM) 


The appropriate generalizations of these ideas to the case of finite 
system have been considered by Bethe (1956) and by Eden (1956). The 
‘suitability of any particular prescription for describing the model wave 
function depends ultimately upon an examination and evaluation of the 
correction terms that subsequently arise. As this analysis involves the 
relation of the scattering matrix to the potential it will not be discussed 
further here ; the form of the equation to be used will be the simplest 
possible generalization of those of the infinite system. 

The scattering described by t is a local phenomenon ; roughly speaking 
i(k) may be taken as the matrix element of a ‘ pseudopotential ’ which is 
its Fourier transform, and if the approximation (2) is used this potential 
will be described by delta-functions and differential operators; the 
- coefficient of k? is a measure of the range of the potential. To analyes 
_ this more clearly it is desirable to consider the complete scattering matrix 
for the two nucleon system, regarding the nucleons as distinguishable ; 
suppose then that the scattering from relative momentum k to k’ has the 
amplitude 
(AStETD ky emer ee tes he ns 1 (8) 
where 7’ is also an operator in spin and isotopic spin. Then the amplitude 
used in the description of the infinite system is 
t(k)=(k|T|k)—(—k|TP,P,|k) . . - - «+ + (9) 


+ It is not necessary to exclude the terms with i—j, because these will give no 
contribution when t properly includes the effects of exchange. 
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. . . “e ‘ ™ 

where P., P. are the spin and isotopic spin exchange operators, and 7 
oe MS . . . . = 

denotes the average of 7’ over spins and isotopic spms. t(k) may be 

regarded as the diagonal element of 


— 


(k’ [él (k' [7 ]k)—(—k' [TPP |k) . 
and it will be supposed that this can be represented approximately by 
(k’ |élk)=to+tyk?ttyk’k 2. . . . . = (LY) 


where f,, and t,, are then the parts of t, associated with scattering in even 
and odd states respectively. This amplitude (11) is equal to the direct 
matrix element of the ‘ pseudopotential ’ 
iy (x—x')— Ht, [S(x—x’) (V—V'P-+(V—V'}? 8(a—2)] 
=ty(V—V’) . o(x—x') (VV. 2) 

The finite system will be discussed by using (12) to evaluate the effective 
potential energy of a determinantal model wave function. The coeffi- 
cients tf), etc. may be functions of position; it will here be assumed that 
they are the same functions of the local density as for an infinite system, 
It will be shown below that it is not necessary for consistency to suppose 
that these coefficients vary greatly with density, so that the assumption is 
not obviously unreasonable ; further discussion would involve the (real) 
potential and the correction terms of the theory. 

The model wave function is assumed to be a determinant of orthogonal 
(real) single-particle states ¢, (7). The local densities p and 7 are defined 
by the natural generalizations of the definitions (3), 


p=Z,1b(@)P; t=S(V@)P. . . . . OB 
The interaction energy of the system is the expectation value of 7’, eqn. (8), 
for the determinantal wave function, and this is equal to the expectation 


value of (12) for the simple product wave function of ¢,. The interaction 
energy density is easily evaluated in terms of p and 7, and is 


atop? tpt +t [2 V (p)*—3p°Vp]—FtaolVp) - - - (14) 
If it is supposed that ¢,, is a function only of p the term in V2p can be 
transformed into a multiple of (Vp), and the total energy density written 
as H=[1/(2M*)]r+4op?+4B(Vp)®. . 2 2. . (15) 
of which the first two terms are identical with the energy density of the 
infinite system. In particular if it is supposed that t, is constant, and 
the scattering only in even states, then 
B=(3/2)t, and 1/(2M*)=1/(2M)+3pB. . . . (16) 
It is with the effects of the last term in the energy that this paper is 
chiefly concerned ; it has manifestly the form of a surface energy, and 
appears in the present treatment as a consequence of the higher momentum 
components needed in the wave function to describe the surface. 
The statement of the problem in the form (15), leads immediately, by 
variation of the total energy, to a set of self consistent equations, of the 
Hartree type, for the best set of single particle wave functions > (x) 
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The single particle states are the eigenstates of the single particle 
Hamiltonian 
h=—V [1/(2M*)].V+U(@). . . . . . (14) 


with the mean central potential defined self consistently by 
d afl 


M™ is defined by eqn. (6) and p and 7 are functions of x given by (13). 


§ 4. FURTHER APPROXIMATIONS 
The solution of the self consistent field eqns. (17) and (18) would require 
much computation ; to get a simple qualitative picture of their implica- 
tions it is necessary to make further approximations. The obvious 
approximation in the present context is that of Thomas and Fermi, 
relating the local kinetic energy to the density ; this will be used here in 
the form of eqn. (7) 7=(3°6)p58 eR eel 


and with this substitution eqn. (15) leads to an ordinary differential 
equation for the density distribution. 

To keep the algebra as simple as possible, a further approximation has 
been introduced to represent the volume energy terms. These are 

Ff = pH (p)={1/(2M)+ pty}(3-6)p®3-+ dtyp? - . (19) 
and the minimum of H(p) should give the density and volume binding 
energy per particle; the experimental values of these quantities are 
CES po=(4r7,2/3)-1, with ry>=1-2x10-Bem;. . . (20) 

€y)=15 Mev. 

As mentioned in § 2 the coefficients must in general be functions of the 
density, but it is not in fact necessary to assume a strong dependence to 
achieve the saturation values given by (20). Again for simplicity, it will 
be assumed that ¢, is a constant, and the value that will be chosen later 
(from consideration of the surface effects) is such that 

Pe OM CML) le Pee. ge (21) 
Then the values (20) will be obtained at the minimum of H(p) if 
—Ey=E (po) = {1-5/(2)}(3-6 po") + 3foPo 
0 = pol" (pq) ={1-5](2M)}(3-6p92)-+ Hopo t+ Ho'po? - (22) 
where i,’ is the value of df/dp at p=p ; and these equations then 
determine Deas Lee inet Hy 8 9 oe. (23) 
These values suggest that the variation of ty with density is not excessive. 

As it is not the purpose of this paper to discuss the relation of tf) to 
the internucleon potentials, a simple form has been chosen for H(p) that 
satisfies the conditions (20) and also that H(0)=0; such is 


F(p)=2y {1 plea es eee (24) 


442 
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The form of H(p) for low densities is not important in the sequel, and the 
most critical assumption hidden in the form (24) is the value of the 
compressibility coefficient ; if this is defined, as by Feenberg (1947), 
as 9p2d2#/dp2, its value is 270 Mev ; in terms of fy, as in eqn. (22), this value 
would imply a second derivative given by 


40,2)’ = — 0-47 By: « ur 
The Hamiltonian density of the system is therefore taken as 
H={1/(2M*)} (r—3-6p°)-+pE(p)+4B(Vp)? . . . (26) 


and the first term is dropped when the Thomas—Fermi approximation is 
made. Ifthe last term were absent, this approximation would lead toa 
sharp edge to the density distribution, and would certainly involve con 
siderable error near the surface. In the present treatment however the 
presence of the last, potential surface energy, term already ensures 
finite density and potential gradients, so that the error involved in the 
approximation is smaller. Indeed the calculations made by Swiatecki 
(1951) indicate that the contributions of kinetic energy to the total surface 
energy are given quite well by the Thomas—Fermi approximation for 
surface thicknesses of 210-1! cm or more (see fig. 2). The validity of 
the approximation will be discussed further in § 7. 


§ 5. SURFACE THICKNESS AND ENERGY 


With the above assumptions the energy density is given by 


H=pH,{(\—p/po)"—l} +a BV ph) 2 ee 

It is now convenient to introduce the dimensionless variable 
Y= plpg 9. Oh 

and the length b defined by 
b°=$Bp [Ey :. > 2 oon 
so that . 

H=p EE ,{y®—2y?+b7(Vy)*}. . . . . . . (80) 
The total energy must be minimized subject to the normalizing condition 
[yd (47/3)rQ2 A) AT, RS a 


(for a nucleus of a large number A of nucleons), and this leads to the 
differential equation for the density distribution 


2b2V%y=d/dy(y®—Iye+rAy) . . . . . . (32) 


where A is the Lagrange multiplier to be determined by the condition (31). 

It is evident that in a large nucleus the left side of (32) must be very 
small at the centre, where also y1, so that \@1 ; and the error made 
by putting A=1 is in fact of the relative order A-2/, Again for considera- 
tion of surface effects V2 may be replaced by d?/dx?, where 2 is the distance 
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measured inwards from the-surface ; these assumptions together are 
equivalent to the discussion of a semi-infinite slab as made by Swiatecki. 
This gives b2(dy/dx)?=y(1—y)? 

y=(tanh 2/2b)*=F(x), . . 2. . . (83) 


Then it will be approximately true for a large spherical nucleus of outer 
radius R that 


y(r) =F (R—?). ee ee ee 4) 
The total energy of the system is then 
[Hr =—AE, +2) 3B(V pdr Aimee Self P35) 
of which the second term represents the surface energy, equal to 
2H, | b2(dy/da)? (4nr,2A23)da—=H,A? (8b/5ry). . . . (36) 
Agreement with the experimental value requires therefore that 
De (S/o)re=010 «10722 Cran a, Be) (87) 


and with this vaiue of b it happens that the surface thickness also agrees 
with experiment. A reasonable measure of the surface thickness is 
provided by the value of fy(1—y) dx; and if comparison is made on 
this basis between the distribution (33) and the often used form 
y= [1+ exp (—2/a)]~1, the corresponding value of a should be 26/3, 
i.e. 0-5 10-4 cm, in good agreement with the evidence from electron 
scattering. 

It should be emphasized that this excellent agreement is rather 
fortuitous and linked with the assumption of the simple form (24) for 
E(p); reasonable changes in the latter might alter the values of surface 
energy and thickness by some 20% ; it would then however be difficult 
to choose a value of b to obtain agreement for both, and the above result 
therefore suggests that the value of the compressibility coefficient implied 
by (24) is near to the truth. 

Finally the outer radius R determined by (31) and (34) is found to be 


R=r, A484 2b=(1-2A¥84-1-5)x10-B em. . . . (38) 


The r.m.s. radius has the pure A? dependence. 


§ 6. THe PoTeNTIAL 
It is interesting now to examine the mean potential U, given by eqn. (18), 
of which an estimate has been made by the Thomas—Fermi approximation. 
When the energy density is written in the form (26), with constant B, 
eqn. (18) should be replaced by 


d 6-0 
Ggre Pl BV ie oat . Saar) . : . (39) 
and for the actual density distribution satisfying eqn. (32) this becomes 
U=—H,—{6-0/(2M*)}p®? . . . . . . (40) 


inside the nucleus, and of course U=0 outside. A noteworthy feature of 
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this potential is the sharp edge at r—R, where the potential well is im- 
mediately 15 Mev deep. This result is independent of the detailed 
assumptions about the form of H(p). Taken at its face value this would 
remove a difficulty, the apparent extension of the potential field signifi- 
cantly beyond the density distribution, that has been discussed by 
Drell (1955). This conclusion however may be rather sensitive to the 
use of the Thomas—Fermi approximation, and also to the representation 
of the pseudopotential in the form (12). 

If the relation (16) is also employed, involving the additional assump- 
tion that scattering in even states is dominant, then the choice (37) implies. 


MjM*=110-5y |. (ee 


and the potential (40) is completely determined ; the central depth of the 
well is V)=—67 Mev. This potential and the density distribution are 
shown in fig. .1. 

Consider now the effect of this potential on an incident nucleon. 
Leaving aside for a moment the finite surface thickness, the particle will 
satisfy the wave equation 


{—(1/2M)V?—E}P=0 9 
outside the nucleus, and inside 
{—(1/2M*)V?— V,—E} ¥=0 hte welt en ea 
which may also be written as 
{—(1/2M)V?—(M*/M)V.+H(1—M*/M)-- EE} ¥=0 . (44) 


so that the effective potential acting on the nucleon, when its inertia is 
supposed constant, is 


Verr=—(M*/M)V.+H(1—M*/M) . . . . . (45) 
and with the corresponding numerical value of M*/M=0-67. 
V ore= —(45—0-33H#) Mev ..... . (46) 


which is in good agreement with the experimental data at low and 
moderate energies (Adair 1954, Feshbach et al. 1954, Melkanoff et al. 
1956). 

In speaking of this effective potential it must also be remembered that 
the boundary conditions at the nuclear surface are changed. The 
proper conditions are that Y is continuous and that 


(1/M) (dP /dr) y= (1/M*) (dP /dr) ine 5 
this discontinuity in gradient has the effect of strengthening the wave 
reflected from the surface. When the finite thickness of the surface 
is taken into account this means that the edge of the potential well 
appears a little steeper than it actually is, as may also be seen from the 


following analysis. In the neighbourhood of the surface the wave equation 
may be approximated in the one-dimensional form 


{—d/da (1/2M*) djdx+U(a)—E}W=0 . . . . (47) 
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and this may be transformed into a Schrédinger equation with constant 
mass by the transformations 


de=fidey=f 
with pea ( MIMS IAEO Gy ee ,  , (48) 


The resultant wave equation for y as a function of Z is 


Belt: +.4 U(x) : ya eapdf is ( 
The potential Uege occurring in this equation is shown in fig. 2 as a 
function of z in comparison with U(x) ; the main effect is a foreshortening 
one due to the change from x to z. 

This phenomenon would be observed as an apparent difference between 
measurements of the surface thickness obtained from electron scattering 
and from the naive interpretation of nucleon scattering. The effect 
expected is of the order of 10°% ; experimentally, the thicknesses deduced 
from the interpretation of proton scattering experiments appear slightly 
smaller than those deduced from electron scattering (Melkanoff et al. 
1956). 


§ 7. CONSISTENCY 


The potential U(x) has been estimated by the Thomas—Fermi approxi- 
mation, and it should now be asked whether the actual density 
distribution corresponding to the use of this potential would reproduce the 
density shown in fig. 1. 

It can be seen from fig. 2 that the variation of M* has a rather small 
effect on the shape of the potential. The topmost occupied level has 
an energy—F,, so that the depth of the potential, measured from this 
level, is equivalent, very nearly, to the potential (M*/M){£,+U(x)} 
for a particle of fixed mass M ; this potential is also shown in fig. 2. 

According to eqn. (40), 


—2M*(E,+U)=(6-0)p2 Renae teen (50) 


which is another expression of the Thomas—Fermi approximation that 
has been made. Swiatecki (1955 a) has suggested an improved form of 
this relation between density and potential depth, derived from considera- 
tion of the density in a potential of uniform slope. For the present 
problem, with the effective potential given by eqn. (50), this relation 
(eqn. (6) ) takes the form 


plpp=(2-v/3)-? (ro/A°P{1-03 (Ura)'y23} 2... (51) 


where ¢(é) is a function given numerically by Swiatecki, and / is a 
measure of the inverse slope of the potential ; in this case 


(ro/2)8=1-08 ry (dy*8/dx)=1-10 y¥6 (1—y) . . . . (82) 


from eqns. (33) and (37). 
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For large 1, eqn. (51) reproduces p/po=y- The density calculated with 
/ given by eqn. (52) is shown in fig. 3 and compared with the original 
distribution y(x). The difference: is quite small, and the jump in the 
potential U at x=0 can be expected further to reduce the discrepancy. 
The third curve in fig. 3 shows the density distribution corresponding to a 
square well of infinite depth with its edge at x=0, as calculated by 
Swiatecki (1951). 


Fig. 3 
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§ 8. SUMMARY 
The following assumptions have been made in this simplified model : 
(i) The general assumptions of Brueckner’s self-consistent field model. 


(ii) The description of the two-particle scatterings in the finite system. 
by a point interaction derived from consideration of the infinite medium, — 


(iii) The retention only of quadratic terms in the velocity dependence 
of the potential, and that the coefficient ¢, is constant. 

(iv) The replacement of the volume energy term by the approximation 
(24). 

(v) The use of the Thomas—Fermi approximation. 

Assumptions (i) and (ii) are the essential hypotheses upon which the 
present model is based, and no attempt has been made to relate the 
‘ pseudopotential’ used here to any actual potential. Assumptions 
(iii) and (iv) are made primarily to simplify the algebra. Assumption (v) 
is made to obtain a simple picture of the form of the solutions of the self 
consistent field eqns. (17) and (18); the arguments § 7 suggest that it is 
not seriously in error. 

There are three remaining parameters in the ‘ pseudopotential ’ that is 
employed ; two of these Hy and py have been chosen to accord with the 
experimental values of the binding energy per particle and with the 
observed nuclear radius. The remaining parameter B and 6 can then be 
chosen to give both the correct surface energy and surface thickness. 
(In the discussion of the surface thickness the distinction between the 
model and true wave functions has been ignored.) 
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With this simple form it is also found that the potential well has a 
depth of 15 Mev at a distance 1-5 10-1 cm beyond the mean radius ; 
this result may be sensitive to assumption (ii). The effective central 
well-depth for nucleons of incident energy # is (45—0-33 E) Mev, in fair 
agreement with observation. The model also suggests that the potential 
should appear to have a slightly thinner surface for incident nucleons 
than would be expected from the measurements of density distribution by 
electron scattering. 
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ABSTRACT 
Analysis of the nucleon-nucleon scattering around 100 Mev has 
determined the spin-orbit coupling part of the two-body scattering 
matrix at that energy, and a reasonable extrapolation to lower energies 
is possible. This scattering amplitude has been used, in the spirit of 
Brueckner’s nuclear model, to estimate the resultant single-body 
spin—orbit coupling for a single nucleon interacting with a large nucleus. 
This resultant potential has a radial dependence approximately propor- 
tional to r~ dp/dr, and with a magnitude in good agreement with that 
required to explain the doublet splittings in nuclei and the polarization 

of nucleons scattered elastically off nuclei. 


§ 1. INTRODUCTION 

THE coupling between the spin of a single nucleon and its orbit within 
a nucleus is manifested by the large splitting between the levels of a 
doublet in the nuclear excitation spectrum, interpreted in terms of a shell 
model, and by the strong polarization detected in the scattering of nucleons 
by nuclei. It is well known (Sternheimer 1955 and references given 
therein) that both these effects can be described approximately by an 
attractive potential (c ./)U(r) between a single nucleon and the whole 
nucleus, and that the strength of this potential is some twenty times 
greater than the relativistic Thomas term associated with the mean 
central potential. 

Attempts have been made to explain this large value of the interaction 
as an amplification due to mesonic effects; in terms of radiative 
corrections (Chisholm and Touschek 1953) or in terms of a reduced 
inertial mass of the nucleon (Johnson and Teller 1955, Skyrme 1955). 
Our present knowledge of the mesonic fields within a nucleus is so un- 
certain, however, that an analysis along these lines remains speculative. 
A simpler, and not necessarily inconsistent, approach lies in an analysis 
in terms of two-body interactions. The only type of interaction which 
gives a spin-orbit coupling in the first approximation is the two-body 
spin-orbit force (Hughes and Le Couteur 1950, Elliott and Lane 1954, 
Blin-Stoyle 1955). The introduction of this force is however an additional 
assumption, and attempts have been made to explain the spin-orbit 
coupling in terms of higher order effects of the two-body tensor force, 
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which latter must also be present ; the calculations of Keilson (1951) 
do not suggest however that this is a good explanation, but the evaluation 
of higher order effects is a very complicated matter. 

In the model of nuclear matter that has been discussed by Brueckner 
(1954) and co-workers in a series of papers, the two-body interactions 
are introduced through the two-body scattering matrix. From their 
original viewpoint the scattering matrix was to be determined by the. 
experimental information on nucleon-nucleon scattering ; subsequently 
it was seen, that it was necessary to recalculate the scattering matrix 
for two particles interacting in the presence of the mean field of the other 
nucleons and subject to the exclusion principle. Apart from the 
complication introduced by this extra calculation there is the additional 
uncertainty due to ignorance of the potentials which give rise to the 
scattering and of how they may themselves be modified within nuclear 
matter. In this paper we have adopted the original point of view with 
regard to the relevant part of the scattering matrix. 

An approximate determination of the scattering matrix has recently 
become possible as a result of experimental and theoretical developments. 
The angular distribution of the n—p polarization at 95 Mev has been 
measured by Hillman and Stafford (1956), and this has made possible 
a phase-shift analysis at that energy by Phillips (private communication), 
these phase-shifts are in close agreement with those given by Feshbach 
and Lomon (1956), whose results give also an estimate of the energy 
variation of the phases. The energies important in the present analysis 
lie below 100 Mev, and approximate analytical formulae have been 
determined to represent the spin-orbit part of the scattering amplitude 
over this energy range ; the uncertainties in this could be reduced by the 
extension to lower energies of the measurements of polarization. 

This scattering amplitude has been interpreted for the finite nuclear 
system in a manner similar to that employed by Skyrme (1956) in his 
discussion of the nuclear surface ; as in that work the present discussion 
assumes a standard nucleus with equal numbers of protons and neutrons 
and does not allow for a possible difference between their distributions. 
Our analysis is closely related to that for a two-body spin-orbit potential, 
and may in some respects be regarded as an extension of the work of 
Blin-Stoyle (1955) ; the relation between these two approaches is discussed 
below. 


§ 2. Tob ScaTTERING AMPLITUDE 


The amplitude / that describes the scattering of two free nucleons is 
a scalar function of k, k’, the initial and final relative momenta of the 
nucleons, and of their spins 6, and o,. It is unnecessary to introduce 
the isotopic spin, provided that the initial and final states are properly 
antisymmetrized, because the p—p amplitude is part of the n—p one. 
The possible forms of spin-dependence of f have often been given, e.g. by 
Stapp (1955); those parts independent of spin are obviously irrelevant 
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to our purpose ; terms that are bilinear in the two spins also give no 
contribution to the resultant single-particle spin—orbit coupling, because 
the direct integrals vanish after averaging over the spins of one nucleon 
and the exchange integrals only give spin independent contributions after 
averaging. The only term that remains is that one linear in the spins, 
of the form 

A=1(6,+0,) . (k’X k)F(k?, k’2, k . k’) ee es 
where F is a scalar function of its arguments; on the energy shell 
Kem Kki2: 

For the comparison with the analysis of the experimental information 
the amplitude (1) must be related to the phase shifts. In terms of the 
scattering angle #, and with the notation of Stapp (loc. cit., eqn. (10) and 
table B) 

K? sin OF (k?, k?, k® cos 0)=—4»/2(Myy—My,). . . . (2) 


The matrix elements M are expressible as series of harmonics with 
coefficients depending upon the phase shifts; if only those terms are 
retained that appear to be important over the energy range considered, 


My,=(1/40k)(sin 6/4/2)[3R,,1—3R,2—+/6R? 
+ cos 6(3.Ro1+5Ro)?—8Ro3+34/2Rk!—44/3R)] 

M y= (1/4tk)(sin O/4/2)[—2R1)°+2R1 > — /6R? = 63) 
+ cos 6(—6R,,!+ 6R#,,°4+34/2R1— 44/33) 


where the R coefficients are functions of the phase shifts 6 and admixture 
parameters « (arising from the tensor interactions) given by Stapp 
(eqn. (39)). In the present problem it is not however really the scattering 
amplitude that is wanted, but the ‘ reactance amplitude ’ corresponding 
to a scattered wave orthogonal to the incident wave; this is obtained 
from the scattering amplitude by replacing [exp (276)—1] by [27 tan 4]. 
Then, for L=/J, 

Me Apa gE ASN, OU AS OE wees ua hh aire 4) 


and for the mixed triplet states with L=J-+1, 
| RY, y= 26(cos? <7 tan 84.,,+sin® <7 tan 55,1) 
R’ =i sin? (tan 67 ,,—tan 84_,). Ped con cane 80 J (5) 


(The quantities R7 do not appear in the expression (2).) 

We want now to represent A as a simple analytic function of the 
momenta k and k’. Comparison of (1), (2) and (3) suggests that we 
should take F as a linear function of k . k’ at a given energy, 


F=A (ke, k2)+0(k2, k2)k.K. ee. (8) 


In principle the values of A and C on the energy shell, when k?=k’”?, 
can be found from the experimental phase shifts (by using formulae (2) 
to (5) above), but at present information sufficient to determine the phase 
shifts is available at only one energy, around 95 Mev, in the important 
range, so that this must be supplemented by semi-empirical estimates of 
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the energy-dependence of the phase shifts ; for this we have used the 
work of Feshbach and Lomon (1956). 

The phase shifts that occur in R can be expanded in a series of odd 
powers of k, so in principle A and C in eqn. (6) could be expanded, on 
the energy shell, in powers of k?. As a rough approximation A can be 
represented by a linear expression, and this will naturally be generalized 
to points off the energy shell, using the condition of hermiticity, by 
assuming 


A(K, 9) =Ay—HAVPE+ AH). 2. 2 (7) 


The energy variation of C cannot be represented properly in this way, 
owing to the behaviour of the 3S phase-shift 6)'; this phase goes to 7 at 
zero energy and must pass through 47 somewhere around 20 Mey, at 
which point C will have a pole. Such low energies do not however make 
a significant contribution to the spin-orbit coupling so we look for a rough 
approximation to C over the range 40-100 Mev. In our subsequent 
analysis we shall first make the crude approximation 


O(i*, k'*)=constant- <2 ee 


which is very poor, but sufficient as a guide to the expected effects of the 
C term. 

It will appear that in fact the most important part of F is the 
‘ forward-scattering ’ part, and it is convenient to include in this part 
the contributions that will arise from exchange. Then, in a ‘ standard ’ 
nucleus with four nucleons in each orbital state, the amplitude concerned 


is 4 
Fy =F (0=0)-+4F @=7) 
= (3/2)A(k?, k®)4+-(1/2)0(k%, W)K2 . . (9) 


where the factor } is the expectation value of the charge exchange 


operator P,. Itis possible to handle somewhat more complicated approxi- 
mations to fy; in particular a quadratic polynomial in k? seems to be 
a fair enough approximation over the important energy range, and the 
inadequacy of experimental or theoretical knowledge of the correct 
value of F does not justify more precise calculations at present. Such 
an approximation is used for our numerical estimates in § 4. 


§ 3. NucLEAR CALCULATION 


The method of calculation is similar to that employed by Skyrme 
(1956) ; as in eqn. (12) of that reference a pseudopotential is introduced 
defined by 

Vi2=(—4h?/M)A(k, k’)8(x,—x,) . . . . (10) 


where A is the reactance amplitude, and the normalization is so chosen 
as to reproduce the scattering amplitude in Born approximation. In 
a matrix element of V 12 k and k’ are to be interpreted as differential 
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operators acting upon the states to the right and to the left respectively, 


k=}(k,—k,)=—4i(V1— Vo), on the right 
k’=3(k,'—k,)= H(y,—¥,), on the left \ aL) 
In a matrix element integrated over all space, integration by parts leads 
to a relation between these operators, analogous to conservation of 
momentum, 
k,’+k,’=k,+kg. SA emer Ate 
In the spirit of Brueckner’s model of the nucleus, the contribution of 
the spin-orbit interaction to the total energy of the system is found from 
the expectation value of > JV,, for a determinantal wave-function of 
tw 
single-particle states ¢,(x). The terms that involve a particular ¢,(x) 
are 
Di [J Po* ag *(% 2) V ral bo(%1)h (Xz) —bo(%a)hi(X1)] dx, dx, . (13) 
and the summation over states 7 may include the state ¢, without error 
because of the antisymmetry. 
When the amplitude A is expressed as a polynomial in k and k’, the 
direct integrals in (13) will lead to sums such as > | 6,(x) |?, ¥4;*(«) 74,(x), 
etc.; those that occur in our subsequent analysis can be expressed in 


terms of 
p= 2 | ¢,(x) P 
Tuv= > (db ;*/da,)(d¢d da ,), T=Thp ree Mae 
f= > (vb*)(V76,) 
The exchange integrals involve similar sums that are outer products in 
spin and isotopic spin space, instead of the inner products that occur in 
(14). With the assumption of a standard nucleus with four particles 
in each state, this outer product will be a unit operator J in the spin 
spaces, so that 
Yd AIH DAW )09. 2... (5) 
Finally in this perturbation calculation all these sums may be evaluated 
for the uncoupled states, and provided that the spin states are equally 
occupied these may be taken to be real in space, so that the conjugation 
sign may be dropped in the above definitions. . 
As a first step the energy (13) has been evaluated with the assumptions 
(7) and (8); we define the constants a, 6 and ¢ so that 
(—4nh?/M)A =i(o,+6,) . (k’ x k)[a—}(bk?+-6*k”)+ck .k’]. . (16) 
The calculation is straightforward ; the use of the relation (15) enables 
the exchange integrals to be expressed as multiples of the direct ones : 
for the a and 6 terms this multiple is }, for the c term it is —}. The 
result may be expressed in the form 
Peet OOV (X)o, eX. eee ey el! 7) 
where then V(a) is the resultant single-particle potential. Provided 
that the sums (14) have spherical symmetry, 1.e. 
(x)=p(r) 
ral =G)eUrB part TH(7) (ey 19 of3) es) 
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V(x) can be expressed in the form (o . 1)U(7). The expression for U oan 

be simplified further if it is supposed that the radial part of do is real ; 

this restriction essentially means that ¢9 describes a stationary state. 
Then we obtain from (16), 


U(r) =(Za)r-tdp/dr 
— (3/16)(b—e/3)|(5/3)r—dz/dr — (2/3)rdz,/dr 
— (10/3); —r-4dp/dr(-72)] 
+-(3/32)br-1d( 2p)/dr 
++ (3/32)(b+ 2¢/3)[(rd/dr+-5)(r-1d/drr-tdp/dr)]  . . (19) 


where } has been written for the real part of b, since only that occurs. 
The last term can be simplified if only terms of order 1/r at the surface 
are retained ; then the last two terms together give 


(3/16)(b+¢/3)r—ld( V2p)/dr. . . . «. . (20) 


These last terms involve second or higher derivatives of the density, 
so that their average through the surface vanishes (for sufficiently large r) ; 
and as might be expected their contribution to the spin—orbit coupling 
is rather small; it is also rather uncertain, depending sensitively upon 
the density distribution and upon the shape of the wave function ¢5. 
If these are neglected the remaining terms of U can be obtained in a simple 
way that can conveniently be extended to more complicated expressions 
for the amplitude A. This method depends upon the observations that 
the vector product k’ xk is equivalent with (k,—k,’) xk, on account of 
the relation (12) that may be used inside the integral (13), and that the 
factor (k,—k,’) is equivalent to an operation of differentiation upon the 
product ¢;*(v)¢,(~). It is this factor which produces the first derivatives 
occurring in the principal part of U ; the later terms containing higher 
derivatives arise from the difference between k, and k,’ in the terms of F. 
Thus the retention only of the leading terms of U is equivalent with 
putting k,—k,’, and so k=k’, in F, after the inclusion of the exchange 
contributions. 

This last prescription means that then only the ‘ forward’ part Fy of 
F is required as in eqn. (9). Then the evaluation of (13) involves first 
the averaging of kFo(k*) over ky at the position x ignoring the gradients 
of density ete., and secondly a differentiation with respect to position 
to take care of the factor (k,—k,’); the result can be expressed in terms 
of gradients of the quantities defined in eqn. (14). It can easily be 
verified that the assumption F')=(3/2)(a—bk®)-+1ck2 indeed reproduces 
the first two terms in the result (19) that we found above. 


As was mentioned in the previous section a fair approximation can 
be made to F,, in the form 


(—40h?7/M)Fy»=a—Bh?tykt . . 1... (21) 
and then this method can easily be applied to find U(r). If the 
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non-diagonal terms such as 7, are omitted the result is 
O(r)=ter—dpjdr 
—(1/8)8[(5/3)r—4d7/dr —r—ldp/dr( 7 2)] 
+(1/32)y1(7/3)r-Mdé/dir— (14/3)rAdr/dr( 2) 
a CPST ATA, Seanad hae aen ee eet (22) 
As written this expression is not properly hermitian but this does not 


cause any difficulty in the following applications where V? is equivalent 
to a real operator. 


The evaluation of these formulae requires a knowledge not only of the 
density distribution, but also of the distribution of + etc. For this 
' we have employed the Thomas—Fermi approximation, as in Skyrme 
(1956), according to which 

7=(3/5)(6-0)p9>;  E=(3/7)(36-0)p7.. (28) 
Furthermore with the same approximation for a particle at the top of 
the Fermi distribution, energy — EZ, 
See) pee Ses FA) (24) 
so that for a particle at energy level H 
— y2=(6-0)p2/8 42 
k®?—=(2M/h2(E+H,). . . . . . «© (25) 

With the help of (23) and (25), U(r) defined by eqn. (22) can be written 

in the form 
U(r)=3r-*(dp/dr) {x —B(4-0p?/9+ $k?) 
+y[18-0p!/?-+ (10-0p?/2)(fk2) +-(fh2)2]}. 2. (26) 

Compared with eqn. (19) this neglects among other terms, those in 7, 
(defined in eqn. (18)). An estimate of the error made can be obtained 
by considering the momentum distribution in a Fermi gas bounded by 
an infinite plane wall, a model suggested by Swiatecki (1951). The 
anisotropy of the momentum distribution is most marked actually at 
the wall, where the density goes to zero; at this point the effect would 
be to replace the factor (5/3) that multiplies r—'dz/dr in the second term 
of (22) by 1:4, and correspondingly in eqn. (26) the term 4-0p”/? that 
multiplies 8 would become 3-6p?*. Since this is only a 10% effect at 
its most important it is a small matter compared with other errors that 
can arise from our approximations. 


§ 4, Toe DouBLET SPLITTING 


We consider in this section the doublet splitting for a nucleon at the 
top of the Fermi distribution, for which, in eqn. (25), k2=0. Then the 
formula (26) gives 

U(r) =4r-Mdp[ar[a— (1-1y?/3)B-+(1-3y#)y] . . + (27) 


where we have put 
p=Ypo: po=(4a7,3/8)-2=07188 5. + (28) 


with r,=1-2, and the unit of length equal to 10=)* com. 
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If the form (21) is assumed to be a good approximation the result (27) 
shows by comparison that the important range of k° is from 0 to 1-1, 
corresponding to laboratory energies of 0 to 90 Mev, with the greatest 
contribution near y=4, around 55 mev. A more physical argument is 
that the important values of k? are those around the average of 4(k,—k,)’, 
where k, is the momentum of the particle at the top of the Fermi distribu- 
tion and k, of another with which it interacts ; this average has the value 
(2/5)kp?=0-63, with rp>=1-2, giving an energy of 52 Mev. 

Unfortunately polarization experiments have not yet been done at 
this low energy, so an extrapolation downwards has been needed of the 
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phase-shift analysis available at 95 mev (Phillips, private communication). 
This has been made following the estimates given by Feshbach and Lomon 
(1956), which are in agreement with the phases determined at 95 mev 
and higher energies. The function F4(k?), defined by eqn. (9), determined 
empirically in this way is shown in fig. 1. An approximation to this 
curve for the lower energies has been taken in the form (21) with the 
coefficients 

a=105,  Pp=y=20 Mey =.) J ae 
with 10-™ cm as the unit of length ; this approximation is also shown in 
the figure. 

With the values given by (29), U,(7) becomes 


U(r) =r-dp/dr[50-24-18(y?3—0-42)2], 2 2. (30) 
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The second term is rather small, reflecting the small dependence of ee 
on k? in this range ; a simple average over y; from 0 to 1, gives the value 
1-3. Since the estimation of its expectation value is rather uncertain 
anyway, we shall adopt the simple formula 
Oe j—olbrdpjdr se BT) 
With the perturbing potential o .1U(r), the splitting between levels 
with j=/++ is 
Abr (el) (Ur) ye. ee, a ov (32) 
For a light nucleus it may not be unreasonable to assume a parabolic 
density distribution 
Pion lal), eat haihin eh 4a aa 88) 
for which, independently of the state of the nucleon considered, (31) 
and (32) give 
AE=(21-+1)(51-5)(2p,/R2). . . . . . (34) 
If we write the standard radius 7,4 1/3 as Ry, the proper normalization of 
the density (33) gives the condition 


P1/Po=2'5(R/,R)5 ihe ee Saleh ond be ALO) 
so that p,; may be eliminated from (34) to give 
AE =24-8(R,/R)3(21+-1)A-28, . . . 2. (36) 


If now & be chosen so that the mean square radius is equal to that of 
a sphere of radius Ry, (R/R,)=1:18 and 4H=10-7 (214+-1)A-?8. This 
gives for 1°N and 1’O splittings equal to 5-3 and 8-1 Mev respectively, 
to be compared with the experimental values of 6-3 and 5-1 Mev 
respectively. 

Another estimate of the density distribution in light nuclei may be 
found from a shell model with oscillator wave functions: the calculation 
of the expectation value in (32) is then equivalent to the calculation 
made by Lane and Elliott (1954) in the limit of zero range, as will be 
discussed below. The calculated splittings depend upon the scale 
constant 7, of the wave functions (supposed to contain the exponential 
factor exp (—4r?/r,?)), and are for 

ISN; 23°9(79/73)° 5 TCE S: Oats) eee eee eva, (ST) 
For the same value of 7, the ratio of these is close to the ratio found for 
a parabolic distribution. To reproduce the observed splittings it would 
be necessary to assume 
LON 7 = 1-56; WO) verge lhe, a 1S a (88) 
As has been pointed out by Lane and Elliott (1954) such a jump in the 
value of 7, going across the 160 closed shell appears to be necessary to 
explain the Coulomb energy differences of mirror nuclei ; their unpublished 
calculation gives the values 1-67 and 1-90 for the two nuclei, with a ratio 
similar to that required in (38). The equation of the mean square radius 
with (3/5)R,? for 1°O gives r;=1-56. 
4B2 
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Since the spin-orbit splitting is such a sensitive function of the density 
distribution, exemplified by the fifth power of 7; that occurs in (37), 
and since it is a rather different measure of it from the Coulomb energy 
difference or the mean square radius, it is hard to say whether or no the 
formula (31) leads to results in accord with experiment. A tentative 
conclusion from the analysis given would be that the constant in that 
formula may be too small, possibly by as much as 30%. 

In the region of heavy nuclei the experimental information near 
208Ph has been examined by Blin-Stoyle (1955), who employed the form 
(31) with an adjustable numerical constant. He found that the data 
could be well fitted by taking for the constant (K, in his notation) the 
value 57 (assuming that r),=1-2), which is 10° greater than our calculated 
value. 

In noting these approximate agreements we should not overlook the 
assumptions made in deriving the simplified formula (26); these were 
first the ‘ forward scattering’ approximation used to obtain (22) and 
secondly the use of the Thomas—Fermi approximation to express all 
quantities in terms of the density distribution. These approximations 
are very similar and both depend for their validity on the smallness of 
the density gradient. This condition is hardly satisfied at the nuclear 
surface, but once the forward-scattering approximation has been made 
the remaining error will be small because of the flatness of fF’) as a function 
of k? ; on the other hand the separate parts A and Ck? of Fo, as in eqn. (9), 
vary more rapidly with energy and the error in the forward scattering 
approximation may be appreciable. It will be shown in the § 5 how an 
alternative evaluation of (16) may be made which indicates, as might 
be expected, that the approximations have tended to underestimate 
the magnitude of the important momenta. 


§ 5. RELATION witH SPrN—ORBIT POTENTIAL 


The method that we have employed is essentially the same as that 
generally used in shell-model calculations, namely the calculation of the 
expectation value of a perturbing interaction in an unperturbed individual 
particle state ; the difference is only that we have used a singular point- 
type interaction instead of a potential. The singularity arises just 
because we have chosen to use a finite polynomial expression for the 
scattering amplitudes ; in so far as these polynomials can be regarded 
as approximations to the matrix elements of some potential, equivalent 
results should be obtained by working with that potential. 

In our problem the amplitude A, eqn. (1), resembles closely the matrix 
element of the two-body spin-orbit potential, and we shall now enquire 


how far this correspondence can be carried. We shall assume a potential 
of the type 


—(1—-%+2P,)(6,+6,) . (r3—Frs) X(Pi—P2)V(r7y2) « - « (39) 


such as has been considered by Lane and Elliott (1954) and other authors. 
In this expression P, is the charge exchange operator and x measures the 
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fraction of exchange potential. The matrix element of this potential 
between states of momentum k and k’ for the n—p system is easily found 


to b : 
8 ie, +-ey) «(kx W)[(1—2) W(| kK |) Fae WER]. (40) 
where W is related to Vir) by 


W()—=(Sn a8) | Vi(r)(sin gr—qr cosqr)rdr. . . . (41) 


In particular for the form of V used by Lane and Elliott, the Yukawa 
form V=V,(A/r) exp (—r/A), we obtain 
WAG == Oman l= eg?) 3222 ole 8 WP er (42) 

As compared with the general form of amplitude given by (1), the 
expression (40) derived from a static potential in Born approximation 
is of a special form, and we cannot expect to be able to express the 
observed scattering amplitude in this way. As we mentioned in the 
discussion following eqn. (6) the coefficient C in that equation has a 
singular behaviour at low energies, and this certainly cannot be reproduced 
from a potential. The physical meaning behind this is that ‘ final state 
interactions ’ play an important role in the spin-orbit coupling between 
pairs of nucleons, and the use of a spin-orbit potential in Born approxi- 
mation can only be regarded as a mathematical artifice. At higher 
energies where the phase shifts become small a direct physical 
interpretation as a potential might be more reasonable, but then the 
whole concept of a potential seems to be unsatisfactory for analysing 
nucleon—nucleon scattering. 

It is possible however to use this idea of a potential to find the effect 
of using the form (19) for U(r). This was derived from the form (16), 
which can be reproduced by expanding W to order g? with suitable values 
for the parameters V5, A, x. Thus to this order 

(1—a)W(| k—k’ |)+ W(| k--k’ ]) 
—167V,A3[1—2A%(k2-+k’2) 4+ 4\2(1—2a)k kJ. 6 we (48) 
The results of Lane and Elliott may also be expanded to relative order A?, 
giving then the splitting due to the amplitude F given by (48); their 
results are for 
15N : AH =(108/4/27)(VoA7/r4°)[1 —[(70+ 40x)/9](A7/7")], 
179: AE=(150/4/27)(VoA5/r4°)[1— (104 3:22) (A?/ry?)]. . (44) 

In terms of the constants a, b and c of eqn. (16), determined by 
comparison with (43), the final factors in these expressions for the 
splittings are proportional respectively to 


a—(b—0:22c)(2-5/r4") ; Gan Dae 140) (2:0/04") eee tn (40) 
which are to be compared with the estimate 
€—(b220-33¢)(4:067) 1. ee one (46) 


derived from eqn. (26). The comparison suggests that the coefficient of b, 
and therefore the nuclear momenta, have been somewhat underestimated. 

As we mentioned before the form (43), equivalent to (16), is not a 
realistic representation of the amplitude so the analysis in terms of a 
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potential is limited in application. Because the amplitude 5 is a slowly 
varying function of energy the analysis in § 4 leads to substantially the 
same result as the use of a potential of zero range (when the exchange 
character is irrelevant). Such a potential implies that all the interaction 
between nucleons takes place in the 8P state, but this is not really so and 
there is a considerable contribution from the 3S—*D coupled state which 
gives the C term in eqn. (6). 


§ 6. PorENTIAL FOR SCATTERING 


Another possible comparison of our results with experiment comes 
from the analysis of the polarization of nucleons elastically scattered 
from target nuclei. Our analysis has been based upon the possibility 
of using the form (6) for the amplitude (although this could easily be 
generalized) ; accordingly comparison should be limited to those energies 
for which it may be a fair approximation. In this section we shall refer 
to the calculations of Sternheimer (1955) analysing the scattering of 
130 Mev protons from light nuclei; this energy is already rather high, 
so that we should be prepared for some inaccuracy in our hypotheses. 

The values (29) do not give a good representation of Fy over the larger 
energy range that is relevant for incident nucleons of 130 Mev ; a better 
overall fit is obtained with the choice 

a=105, Pp=y=8 kw 2 pe 
which is illustrated in fig. 1. 

In this case we must include the k? terms in the formula (26) for U(r) ; 
k? was defined by eqn. (25), and in that equation the mass M should 
properly be the effective mass M*. We adopt the values used by 
Skyrme (1956) for M* and Ey, which give 

Ak2—=1-75(140-5y)-2 . 2 www. (48) 
where as before y is the density in units of py; we also use the same 
density distribution, so that in the ‘ large nucleus ’ approximation 


—4r-ldp/dr=0-0767A-1/3y1/2(1 —y) 


y=plpo=ltanh (R—r/(1-5))P. 5. . (49) 
The formula (26) thus gives 
—AlWU(r)=—K(1—y)y* ee ee ee 


where the coefficient K varies from 8-85 to 10-25 as y goes from 0 to 1. 
For comparison we note that eqn. (31) gives K=7-9, so that there is 
a 20% increase in the strength of the potentials compared with that felt 
by a particle at the top of the nuclear well. 

Before we compare our potential with that used by Sternheimer, we 
should make allowance for the mass variation. In the approximate 
expression for the phase shifts used by Sternheimer the perturbing 
potential is multiplied by the factor (4k/T'), where 7 is the kinetic energy 
of the nucleon and k its wave number. This factor is equal to 
(2M*/(H—V))*/2, where E is the energy of the incident nucleon and V 
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is the nuclear potential, and when the effective mass changes appreciably 
through the surface the potential V is increased so that the variation of 
this factor is not negligible. An effective potential has accordingly been 
estimated by multiplying (50) by this factor, expressed as a function of y 
with the same hypotheses. The resulting modified coefficient K’ then 
varies from 8-4 to 6-8; this potential is shown in fig. 2 as a function of 
distance from the outer edge of the nucleus, and it is not significantly 
different from the potential (31) 

The maximum value of the potential is about 34~-1/3 Mev and its 
width is very similar to that assumed by Sternheimer, whose potential 
is also shown in the figure, adjusted to have a coincident maximum ; 


Fig. 2 
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the maximum value is close to that found necessary by Sternheimer to 
explain the experimental results. The main difference is that the 
maximum of our potential occurs 0-5 x 10~'8 cm outside the mean radius 
r=r,A¥? instead of at that point (always speaking of the potential in 
a large nucleus, ignoring the r~1 factor outside). It is not possible to say 
without detailed calculation whether this and other small differences 
would help to remove the difficulties that Sternheimer found in fitting 
properly the angular distributions. 


SUMMARY 


The general shape of the phase-shift analysis of the nucleon-nucleon 
scattering for energies up to 100 Mev now seems well established and it 
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has been possible to fix the spin-orbit part of the scattering amplitude 
with reasonable certainty. 

This empirically determined amplitude has been applied in the spirit 
of Brueckner’s model of the nucleus to find the effective single-particle 
spin-orbit coupling. The essential assumptions are that 

(i) The energy of the nuclear system may be determined by using 
a ‘model’ wave function which is approximately that of a set of 
independent particles moving in a central scalar potential. 

(ii) The energy change due to a perturbing interaction is approximately 
the sum of changes due to two-body interactions and can be evaluated 
by using the reactance amplitudes and model wave function. 

(iii) The scattering associated with these two-body interactions is 
approximately the same as that between free nucleons, so far as the 
spin-orbit amplitude is concerned. 

The numerical result that we obtain from our analysis appears to be in 
agreement with experimental evidence, within the limitations arising 
from ignorance of the nuclear wave functions. However this agreement 
really involves only one quantity, the strength of the coupling, so that it 
would be premature to claim that the method of calculation was thereby 
justified. 
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Stacking Faults in Cold-Worked Cobalt—Nickel Alloys 


By J. W. CurisTian and J. SpREADBOROUGH 
The Inorganic Chemistry Laboratory, Oxford + 


[Received March 22, 1956 and in revised form May 30, 1956] 


WARREN AND WAREKOIS (1955) have shown that heavily cold-worked 
alpha-brasses contain deformation stacking faults in sufficient amounts 
to give measurable x-ray diffraction effects. The stacking fault para- 
meter (Paterson 1952) increases with zinc content from «0-005 at 
10% zinc to «=0-039 at 35% zinc. In this note, we briefly record 
similar results we have obtained for cobalt—nickel alloys, and we discuss 
a possible correlation between the production of faults and the abnormally 
large increases in electrical resistivity on cold-working, and large work- 
hardening capacities, shown by certain alloys of critical composition. 
Cobalt-nickel alloys were chosen because the free energy difference 
between the bulk face-centred cubic and hexagonal close packed forms 
changes sign in a certain composition range, thus indicating that the 
specific energy of a fault should be low. Electrical resistance measure- 
ments by Broom and Barrett (1953) have already suggested that stacking 
faults probably form in these alloys. Growth stacking faults ((14) 
faults in Frank and Nicholas’ 1953 notation) are found in considerable 
amounts in hexagonal cobalt naturally transformed from the cubic 
phase (Edwards and Lipson 1942, Wilson 1942), but recent work in this 
laboratory has shown that when the martensitic transformation is stress- 
aided, the faults are predominantly of the hexagonal deformation (24) 
type. These (24) faults are the hexagonal equivalent of the cubic 
deformation faults studied in the present work. Alloys containing from 
50 to 69 weight-°/, cobalt were filed, passed through 250 mesh sieves and 
formed into diffractometer specimens in specially designed holders. 
Accurate line profile measurements were made by stepwise counting on a 
Norelco diffractometer, using cobalt radiation. The specimens were given 
various intermediate annealing treatments before final recrystallization 
at 800°, and line profiles were measured at each stage. Values of the 
stacking fault parameter were obtained from the shifts in the positions of 
the Ka, peaks, which were located by the method described by Anan- 
tharaman and Christian (1953). Consistent values could be obtained from 
the movements of individual lines, as well as from the changes in separa- 
tion, and the method is considered to be more reliable than that of 
measuring the centre of gravity of a line used by Warren and Warekois. 


+ Communicated by the Authors. 
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The stacking fault parameters for cobalt—nickel alloys are ue 
tig. 1; the probable error is believed to be less than 10°4 when a~0-02. 
The result for pure nickel was obtained by T. R. Anantharaman in pre- 
liminary measurements, and is probably just significant. The general 
shape of the curve is very similar to that found by Warren and Warekois 
for copper-zine alloys, and there is a steep rise in the concentration of 
stacking fault near the critical composition where the free energies of the 
two phases become equal at room temperature. A similar high value of 


Fig. 1 
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Stacking Faults and Changes in Electrical Resistance in Cobalt—Nickel Alloys. 

x Measured Stacking Fault Parameter. 

4 Percentage Resistance Change after 35°, Reduction in Area (Broom and 
Barrett 1953). 


a(0-024) was obtained from filings of a cobalt—iron alloy containing 7-9°% 
iron, which is near the critical composition in this system. As a check on 
the methods, we measured « in a copper-30% zine alloy, obtaining a 
result (0-023) in very good: agreement with that of the American workers. 
Very careful measurements were made on heavily cold-worked pure 
copper filings, but the shifts in peak positions were so small that only the 
order of magntiude of « could be determined. Our value of 0-004 is thus 
compatible with that of 0-007 obtained by Greenough and Smith (1955). 
It should be noted that the diffraction theory is based on the assumption 
that faults form on one set of planes only, whereas in the present experi- 
ments they are formed equally on all sets of {111} planes. To a first 
approximation, however, the concentration of faults on each set of planes 
may be taken to be «/4. 

Broom and Barrett (1953) found that the increase in electrical resistivity 
caused by cold-working cobalt-nickel alloys had a maximum in the 
region where heavy faulting is expected (fig. 1). This result had been 
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predicted by Broom (1952 a) who first suggested that the large increases 
in resistance found on cold-working certain face-centred cubic alloys 
might be due to the formation of stacking faults. Such faults on one 
set of planes constitute a formidable obstacle to slip on intersecting 
planes, and in fact the alloys concerned are those which show markedly 
unequal hardening of latent and active slip planes. Recent work on 
alpha-brass (Piercy et al. 1955) has shown that there is a hardening of the 
latent planes, as required by this explanation, rather than a softening of 
the active planes. The connection between percentage increases in 
resistivity and hardness at room temperature and the formation of 
stacking faults is also evident for the copper—zinc alloys, as shown in fig. 2. 


Fig. 2 


0:04 
AH 
DH % 
250 
(DQ) 
0-03 
230 
Ox 
210 (0, g 
0:02 
190 
1704 O-Ol 
ISO . 
wt.% Zinc 
“1 I 
Gu 10 2S 30 


Stacking Faults and Changes in Electrical Resistance and Hardness in Copper— 
Zine Alloys. 

© Measured Stacking Fault Parameter (Warren and Warekois 1955). 

x Measured Stacking Fault Parameter (Present Work). 

A Percentage Resistance Increase for Log Strain 2-3 (Boas and Nicholas 1953). 

vy Percentage Resistance Increase for 84° Reduction in Area (Crampton, et al. 
1941). ; 

CO Pants Increase in Hardness after 90°% Reduction by Rolling (Copper 
Development Association). 


The measurements are, of course, made at temperatures where there are 
still some point defects (probably single vacancies) contributing to the 
resistance, but it seems likely that stacking faults are responsible for the 
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rapid variation with composition. A survey of the literature suggests 
that whenever there is a tendency to have low fault energy (as shown by 
transformations to hexagonal phases, or by prolific formation of annealing 
twins), the increase of resistance on cold-working and the rate of work- 
hardening are unusually high. 

Since 4a (110) dislocations are dissociated into narrow ribbons of 
stacking fault on {111} planes, it is tempting to equate the measured m 
values with the fractional areas of the close-packed planes contained 
within these extended dislocations. Any such interpretation is, however, 
highly speculative, since the diffraction theory is based on the assumption 
that faults occupy whole planes, and the diffraction effects of narrow 
ribbons of stacking fault cannot be predicted in detail at present. If 
the separation of the partial dislocations is nevertheless calculated on this 
basis, the results obtained are in good agreement with theoretical estimates 
of the separations (Seeger and Schoeck 1953, Seeger 1955 a). A density 
of 10! dislocations cm~? corresponds to separations of ~8 atoms for 
¢=0-004 (the lower limit observable) and ~80 atoms for «=0-040 (the 
highest value yet observed). The actual variation in the width of an 
extended dislocation must be less than indicated by the « values, however, 
since the dislocation density characteristic of the heavily cold-worked 
state also varies with composition (Smallman, private communication). 
The present position thus appears to be that the diffraction effects are in 
accordance with Paterson’s theory, but the physical significance of the 
faulting parameter is uncertain. An alternative possibility is that large 
faults are produced irreversibly during the deformation by processes such 
as the breaking of Cottrell-—Lomer locks in the way suggested by Friedel 
(1955). 

Our results for the cobalt—nickel alloys, and those of the Americans 
for copper-zine alloys, show that the measurable stacking faults are 
removed gradually as the annealing temperature is raised. Typical 
results from one of our alloys were that about 20°% of the measurable faults 
were removed by a short anneal at 270°c, and 60% at 400°c ; longer 
anneals at these temperatures produced no further changes. The 
faulting parameter became too small to be detected in all alloys after a 
half-hour anneal at 500°c, although appreciable line broadening remained 
to much higher temperatures. These results seem to be more readily 
explained if the faults are extended dislocations. Reduction in the 
numbers of faults would then occur both by the mutual elimination of 
opposite dislocations and (possibly) by the aligning of dislocations into 
sub-boundaries, where their widths would be reduced. There is also some 
evidence (‘Tamman 1936, Maddigan and Blank 1940) that in alpha—brass, 
a major part of the extra electrical resistance anneals out in the tempera- 
ture range where the faults diminish. 

_ Seeger (1955 a, b) has emphasized that materials with very narrow 
dislocations are expected to have plastic properties quite different from 
those with wide extended dislocations. The known differences in the 
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I Al log strain 2-3, D 0-5 Boas and Nicholas 
: é (1953). 
II Cu, OFHC 84%, D 3 Crampton, et al. 
; (1941). 
Cu, spec. log strain 2-3, D 3 Boas and Nicholas 
pure: (1953). 
Cu log strain 2:3, D 3 Boas and Nicholas 
1953. 
Cu 40-80°%, D 2 Takahasi (1930). 
Ni OOo DD 8 Geiss and van Liempt 
1 (1927). 
Ni log strain 2:3, D 6 Boas and Nicholas 
(1953). 
Ni 30% 4 Broom (1952 a). 
Ag 98%, D 5 Tamman and Dreyer 
(1933). 
Ag OUTe, v 3 Takahasi (1930). 
Pt J9%,,.0 6 Geiss and van Liempt 
(1927). 
Pd 96%, D 4. Tamman and Dreyer 
(1933). 
Cu-Ni 
(various 84%, D 4 Crampton et al. 
compn.) (1941). 
Ag—25% Au 97°95. BR 1 Tamman and Dreyer 
(1933). 
Ag-35% Au 90%, D 1 Broom (1952 a). 
III Ag-30°%Cd 97%, R 28 Tamman and Dreyer 
(1933). 
Ag-23%Zn 97%, B 19 Tamman and Dreyer 
(1933). 
Cu-20%Zn | log strain 2-3, D 20 Boas and Nicholas 
(1953). 
Cu-25%Zn | log strain 2-3, D 25 Boas and Nicholas 
(1953). 
Cu-30%Zn 84%, D 26 Crampton et al. 
(1941). 
Cu-5-7%Al | log strain 2-3, D 28 Boas and Nicholas 
(1953). 
Cu-6% Al log strain 1-96, D 27 Broom and Clothier 
(1952). 
Cu-7:5%Al 84%, D 33 Crampton et al. 
(1941). 
Co-31%Ni 30%, 8 Broom and Barrett 
(1953). 


ee en ee eee 
In the deformation column, the percentages give the reduction in area caused 


Group I. 
Group II. 


Group III. 


by drawing (D) or in thickness by rolling (R). 
Fault energy >100 ergs cm~’, dislocation widths ~1 atoms. 
Fault energy ~40 ergs cm~’, dislocation widths ~10 atoms. Faults 
in cold-worked material just detectable by x-rays. 
Fault energy <10 ergs cm-*, dislocation widths ~40 atoms. 
Faults in cold-worked material readily detectable. 
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behaviour of (say) copper and brass, together with the data on resistivity 
and work-hardening, suggest that it may be more convenient to sub- 
divide face-centred cubic materials into three ranges of stacking fault 
energy. Such a classification is shown in the table, but it must be 
emphasized that the alloys in Group II are in critical composition ranges, 
and there may be a gradual transition to Group IT behaviour as either the 
composition or the temperature is altered. Assuming that point defects 
(when present) and narrow dislocations give increases in resistivity of 
only a few per cent, it is seen that the magnitude of the resistance change 
can serve as a rough guide to the stacking fault energy. There is also some 
evidence (Broom and Clothier 1952, Broom 1952 b) that the anisotropy 
of the stacking fault resistance can be detected for Group III materials. 

Quantitative measurements on stacking faults are necessarily made on 
powders, in order to obtain high deformations, and also because line 
shifts in polycrystalline bulk specimens may be caused by residual inter- 
granular stresses (see e.g. Greenough 1952). However, faulting in sufficient 
amounts may be detected from streaks on single crystal oscillation photo- 
graphs, and it is clear that the role of stacking faults in deformation can 
now best be clarified by experiments on single crystals given controlled 
deformations. Low temperature experiments should also be interesting, 
both because the amount of faulting should increase, and because the 
critical alloys of low fault energy are probably the most favourable 
materials for the production of face-centred cubic mechanical twins, as 
observed in copper by Blewitt ef al. (1955). 
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Stark-Splitting in Crystals 
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RARE-EARTH ions in crystals have absorption and fluoresence spectra 
consisting of line multiplets separated by a wave number difference 
~1000 cm-!. The lines within a group may have separations ~100 cm~t 
and are all due to transitions to a multiplet of levels with the same total 
angular momentum (J). The splitting of the levels within a multiplet 
has been examined by Bethe (1929) using the methods of group theory, 
and his treatment of the crystallographic double-groups needed to deal 
with half integral values of J has been extended by Opechowski (1940) 
to include tetrahedral and rhombohedral symmetries. Using a slightly 
different approach Hellwege (1948) has considered in detail the wave 
functions involved in a discussion of the levels in fields having the 
symmetry of any of the 32 crystallographic point groups. These 
theoretical papers do not express the results in a condensed form suitable 
for comparison with experiment ; but Gobrecht (1937) compiled a table 
of the number of levels for a given J in fields of different symmetry, 
and this table has been frequently reproduced (Freed 1942; Yost et al. 
1947; Pringsheim 1949). Unfortunately all these tabulations are in 
error in their treatment of fields of tetragonal symmetry which are 
asserted to produce as many distinct levels as the lowest symmetry 
groups. It is now intended to rectify this error which is not present in 
the original work, and to discuss the splitting of the levels in nearly cubic 
fields, which frequently occur both in pure salts, e.g. alums, and when 
impurity ions requiring charge compensation enter cubic crystals 
(Runciman 1955). 

The 32 crystallographic point groups can be classified into seven systems 
and for the present purpose these will be grouped under four headings 
(a) cubic, (b) hexagonal, containing the hexagonal and rhombohedral 


systems, (¢) tetragonal, and (d) lower symmetry including orthorhombic, 
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monoclinic and triclinic. The point groups in Schoenflies notation under 
these four headings are as follows : 


(a) Cubic :—O,, O, Tg, T), T. 
(b) Hexagonal :—Dey, Dg, Coys Con Co» Dan» Can Daa, Da, Cay, Se, Co: 
(c) Tetragonal :—D,,, Dy, Cay, Cap, Ca, Doa, S,- 


(d) Lower symmetry :—D», Do, Coy, Con, Co, Cy, Se, Cr. 


A recalculation of the splitting for integral J has shown that all point 
groups within one of these four categories produce the same splitting, 
and the number of lines for given J is listed in table 1. The number of 
triply, doubly and singly degenerate levels (ng, n, n,) are written below 
the total number of levels for cubic symmetry. An inspection of correla- 
tion tables for the species of a group and its subgroups (Wilson, Decius 
and Cross 1955) shows that for a cubic field with slight hexagonal distortion 
the triply degenerate levels will split into a doublet and a singlet, and the 
doublet levels will remain doublet ; whereas for a tetragonal distortion 
the triplet levels will again split into a doublet and a singlet, and doublet 
levels will split into two singlets. 


Table 1. Term-splitting for Integral J 


af 0 1 2 3 4 5 6 7 8 General rule 
Cubic ] i 2 3 4 4 6 6 id add 5 every 6 
(Ng Dy 4) (OOL) | (100) | (110) | (201) | (211) | (301) | (312) | (411) | (421) | add (311) every 6 
Hexagonal 1 2 3 5 6 7 9 10 ll add 4 every 3 
Tetragonal 1 2 4 5 7 8 10 it 13 | add 3 every 2 
Lower 
symmetry 1 3 5 tu 9 ll 13 15 LT 2J+1 


Table 2. Term-splitting for Half-integral J 


9/2 | 11/2 | 13/2 | 15/2 | 17/2 General rule 


Cubic 1 1 2 3 3 4 5 5 6 add 2 every 3 
All other 
groups 1 2 3 4 5 6 7 8 9 J+t4 


Levels are counted as doublets, even though the characters for the 
representations are distinct, so long as the characters are complex conju- 
gates of each other, as in this case there is no frequency separation of the 
levels. For orthorhombic and lower symmetry there is complete splitting 
into 2J-+-1 levels for integral J. 

For half-integral J, i.e. for an odd number of electrons, the calculation 
of the splitting involves the double-groups. Bethe’s results have been 
confirmed and table 2 lists the number of levels, either quadruplets or 
doublets, for any of the five cubic point groups. For all other groups 
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there is total splitting into J-++4 doublets, all levels remaining doubly- 
degenerate in any electric field (Kramers 1930). It therefore appears 
as if ions with an even number of electrons are better suited to determine 
the site-symmetry of salts or impurity-activated solids containing 
rare-earth ions. The application of a magnetic field, which in general 


removes all degeneracy, may also help in finding the terms responsible 
for an emission or absorption line. 
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CIX. REVIEWS OF BOOKS 


Eazanding Universes. By Professor E.Scur6.tnczr. (Cambridge University 

Press, 1956.) [93 pp.] Price 17s. 6d. 

ProrEssor Sct RO INGER has produced a most interesting and valuable little 
book. Although the title obviously derives from the modern astronomical 
problem of the universe as a whole, the author is not here concerned with this 
actual problem but only with the purely geometrical aspects of conceptual 
expanding models irrespective of whether their metrics in fact satisfy the field 
equations of general relativity. In other words the book is an exercise in multi- 
dimensional geometry aimed at clarifying those parts relevant to expanding 
universes generally, and anyone wishing to sharpen his wits in this intricate 
realm will find an opportunity in following the quickly moving arguments. 

Beginning with the properties of flat and spherical expanding models, the 
paths of particles and light rays, and the corresponding frequency changes, are 
investigated for them. The properties of waves in general Riemannian space- 
time and in an expanding universe are studied by means of Hamilton-Jacobi 
theory. The book will help to illumine the theory of the actual expanding 
universe in much the way that a treatise on conics might have helped in the 
understanding of Keplerian orbits three centuries ago. 

The subject matter of Professor Schrédinger's book, which is based on a brief 
course of lectures, is so sublime that the most reverent of styles would be appro- 
priate to it, but regrettably the book is marred by rather too many distracting 
informalities and oddities of phrasing, and other such imperfections that could 
easily have been eliminated with great advantage to an otherwise elegant and 
handsome contribution to a highly important and recondite subject. R.A.L. 


The Interference Systems of Crossed Diffraction Gratings: The Theory of Moiré 
Fringes. By J. Gutup. [Pp. vi+152.] (Oxford: Clarendon Press, 1956.) 
Price 25s. 

Tr 1s monograph by Mr. J. Guild of the Light Division of the National Physical 

Laboratory is one of those rare gems—an entirely original book, which breaks 

completely new ground. In this text Mr. Guild develops the theory of the 

production of so called Moiré fringes arising from the crossing of line gratings. 

In particular he discusses the properties, localizations and general characteri- 

stics of these hitherto almost neglected optical effects. As he points out, 

Moiré fringes were first briefly discussed in 1874 by Lord Rayleigh, but the 

subject has remained since then in abeyance. It is only very recently that the 

National Physical Laboratory has vigorously begun to study these fringes with 

the double purpose of (a) studying the characteristics of the gratings them- 

selves and (b) using the Moiré fringes for precision metrological measurement. 

Mr. Guild has supplied a vital part of these new developments. 

This is a small book and highly specialized with six chapters as follows : (1) 
Introduction (2) Fundamentals (3) Gratings in Series (4) Characteristics of 
Moiré Fringes (5) Metrological Applications (6) Grating Defects. 

It is unfortunate that classical optics in university courses now, of necessity, 
occupies much less attention than formerly, owing to the pressure of more 
‘modern’ branches of physics. Thus few undergraduates can be expected to 
profit from this text. On the other hand all postgraduate workers in optics 
and every teacher of optics will welcome this text with plaudits, for the 
author treats the subject in a masterly and comprehensive fashion. Most 
notable is the fact that the author has a sense of the physical feeling for what 
the light waves are doing and his many original and valuable conclusions reflect 
this feeling for the physics rather than the mathematics of the processes he is 
discussing, 
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Mesons and Fields. Vol. 1 (Fields). By 8S. S. Scoweper, H. A. Berar and 
F. p—E Horrman. Vol. II (Mesons). By H. A. Barun and F. pz Horrman. 
(Evanston, Ill.: Row, Peterson & Co.) [Pp. (Vol. I) xvi+449; (Vol. II) 
xiii+446.|] Price 63s. each volume. 


THESE two volumes give a comprehensive account of the properties of mesons 
and of the various theoretical methods which are available for their discussion. 
There are chapters on ~-mesons and on K-mesons, but, naturally enough, the 
main part of the detailed discussion is about pions. 

The first volume deals with what may be called the ‘highbrow’ forms of 
theory ; ice. the general and abstract development of quantum field theory, 
with particular emphasis on the covariant relativistic formulation and the 
renormalization programme. The first section gives the theory of relativistic 
wave equations without field quantization, and includes Feynman’s form of 
positron theory. In Feynman’s scheme positrons are regarded as electrons 
travelling backwards in time, and it is possible to draw the well-known Feynman 
diagrams in which upward pointing lines represent electrons and downward 
pointing lines represent positrons. The discussion of Feynman kernels and 
boundary conditions, which is necessary for this scheme, is fairly subtle; it 
must make heavy reading for any but the specialists. 

The central section of Vol. I is a comprehensive account of the methods of 
field quantization; it includes valuable accounts of Wick’s theorem on 
expressing a product of operators as a sum of creation and annihilation processes, 
and of the Gupta—Bleuler method for eliminating the longitudinal part of the 
electromagnetic field. The remainder of Vol. I is mostly concerned with the 
renormalization procedure. At this point a critical discussion would have been 
very interesting and very valuable ; instead, the authors are content to follow 
the standard references fairly closely. The renormalization idea, originally 
due to Kramers, is as follows: subtracting the infinities arising in field theory 
is equivalent to redefining the mass and charge of the elementary particles. 
Unfortunately the mathematics of this subtraction scheme becomes very in- 
volved; we have to consider expressions which are infinite series of terms, 
each term being a divergent multiple integral. It requires just a Little faith to 
believe that all the subtractions work out as they should. 

The second volume, besides a comprehensive description of the properties 
of mesons, contains an account of the practical phenomenological form of 
theory which is used fairly successfully to give a coherent picture of pion 
phenomena. This material makes Vol. IH an extremely valuable reference 
book for anyone wishing to interpret new results. The simple form of theory 
which is used is mainly based on the pseudoscalar nature of the pion, in variance 
in isotopic space and a pseudovector interaction between pions and nucleons ; 
this approach has been successful in explaining a wide variety of scattering and 


pion production phenomena, in particular by the methods developed by Chew. 


In one section of Vol. II the rigorous field theory methods are-used to analyse 
pion-nucleon scattering by the so-called Tamin—Dancoff method. The analysis 
is complicated, and unfortunately the results are not so satisfactory as those 
coming from the recent phenomenological theory of Low and Chew which is 
briefly mentioned in an Appendix. Here we have an example of the odd 
situation in contemporary tield theory. Since 1946 some very high powered 
mathematical effort has been put into developing the field theory of inter- 
acting elementary particles, but, so far as mesons and nucleons are concerned, 
the practical results are much less than might have been expected. We have 
certainly learned much about. selection rules, but when we have to estimate a 


- erogs-section, the phenomenological theory in many cases gives as satisfactory, 


or more satisfactory, results than the strict application of field theory. This 
does not mean that field theory canbe ignored ; -it remains the only precise way 
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we know for describing elementary particles in accord with special relativity, 
the uncertainty principle and quantum statistics. 

Professor Bethe and his collaborators by giving so thorough an account, 
both of the field theory and the more practical methods, have produced two 
particularly valuable volumes. Aes 


La Résonance Paramagnétique Nucléaire. Edited by P. Griver. (Paris : 
Centre National de la Recherche Scientifique, 1955). [Pp. 198.] Price 
1600 frs. 


Here is a useful book which brings together the contributions to the study of 
nuclear magnetic resonance by workers in France and in Geneva. Although 
covering a good deal of ground the book does not attempt a comprehensive 
survey of the subject, but instead is built mainly round the work of the authors. 
There are eight authors who contribute chapters on their own specialities. 
Coming from so many sources it is not surprising that the chapters are not 
closely integrated and vary in difficulty. Many readers making their first 
serious acquaintance with the subject through this book may find the early 
theoretical chapters rather difficult. The printing is in typescript characters 
and this appears to be the cause of the rather frequent misprints, especially in 
mathematical formulae. 

The Editor, Professor Grivet, sketches the background of the phenomenon 
in an introductory chapter, and points out some of the salient features. In 
the second chapter Dr. Gabillard gives a theoretical account from a macro- 
scopic viewpoint. He develops Bloch’s equations and discusses their solution 
under various conditions. In chapter 3, Dr. Ayant looks at the phenomenon 
from the quantum-—mechanical standpoint and calculates the linewidth and 
relaxation time for liquids and the spectrum for some solids. The theoretical 
concepts would have been made clearer if these two chapters could have been 
dovetailed. 

In chapters 4 and 5 Dr. Gabillard and Dr. Soutif give an account of methods 
detecting nuclear magnetic resonance. The methods described were: the 
crossed-coil induction method, the Q meter method, the bridge method, the 
spin-echo method, other pulse methods, marginal oscillator and super- 
regenerative oscillator methods. Some of these methods are analyzed in 
greater detail and practical circuits are given with values of the required 
components. 

In chapter 6 Dr. Béné, Dr. Denis and Professor Extermann of Geneva 
discuss methods of producing the steady magnetic field. A surprising amount 
of space is devoted to the production of the field by iron-free coils. Electro- 
magnets and their stabilization are discussed, but rather little is said about 
permanent magnets. 

In chapter 7 Dr. Béné tabulates nuclear moment data from all sources and 
Dr. Ayant refers to the shifts and corrections encountered when measuring 
nuclear magnetic moments by the nuclear magnetic resonance method. 
Chapter 8 is concerned with pure quadruple resonance, the theory being given 
by Dr. Ayant and the experimental section, with full circuit details, by Dr. 
Buyle—Bodin. The book closes with a brief review of the theory of nuclear 
moments by Béné, Denis and Extermann. E.R. A 


Progress in Semiconductors. J. Edited by A. F. Gipson, P. AigRatn and 
R. E, Bureuss. (London: Heywood & Co.) [Pp. 220.] Price 50s. 

Tus is a collection of review papers relative to the theory and technology of 

semiconductor devices. The aim in most cases seems to be bibliographical 

rather than expository, but the papers will be convenient as temporary stores 

of information for those active in the field. TMZ 
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Astronomical Optics and related subjects. Edited by Z. Koray. (Amsterdam: 
North Holland Publishing Company, 1956). [Pp. xv+428.] Price 90s. 

Tue Editor and publishers are to be congratulated on this book, which gives an 
account of the proceedings of the Symposium held in the University of 
Manchester in April 1955. 105 guests from five countries attended this 
unusually lively Symposium. The 46 papers here reproduced or summarized 
fall under the following heads: 1. Information theory and optics; 2. Optical 
images and diffraction; 3. Interferometry and coherence problems; 4. 
Electronic devices in Astronomical Optics; 5. Resolution Problems and 
Scintillation ; 6.. Wide-angle optical systems and aspheric surfaces ; 7. Filter 
photography and thin films. Introductory remarks by Professor Kopal and 
concluding remarks by Dr. J. Résch (France) complete the volume. 

K. H. L. 


Elementary Wave Mechanics, 2nd edition. By W. Hurrpyr. (Oxford : 
Clarendon Press. London: Cumberlege.) [Pp. viii+193.] Price 18s. 

Tak second edition of Heitler’s well known book is very welcome. It is about 

half as long again as the first edition ; as it is addressed primarily to chemists 


the expansion is mainly in a section on diatomic molecules and on the chemical 
bond. NSM. 


Surveys in Mechanics : the G. I, Taylor 70th Anniversary Volume. Edited by 
G. K. Barowetor and R. M. Davims. (Cambridge University Press, 1956.) 
[Pp. vii+475.] Price 50s. 


Tats volume, containing ten articles and a biographical note, was prepared to 
commemorate Sir Geoffrey Taylor’s seventieth birthday on the 7th March, 1956. 
Each article surveys the present position in a particular field of fluid or solid 
mechanics in which Taylor has worked sometime during the last forty or fifty 
years. 


The contents of the volume are as follows : 


G.I. Taylor: A Biographical Note, by R. V. Southwell. 

The Mechanics of Quasi-Static Plastic Deformation in Metals, by R. Hill. 

Dislocations in Crystalline Solids, by N. F. Mott. 

Stress Waves in Solids, by R. M. Davies. 

Rotating Fluids, by H. B. Squire. 

The Mechanics of Drops and Bubbles, by W. R. Lane and H. L. Green. 

Wave Generation by Wind, by F. Ursell. abe 

Viscosity Effects in Sound Waves of Finite Amplitude, by M. J. Lighthill. 

Turbulent Diffusion, by G. K. Batchelor and A. A. Townsend. 

Atmospheric Turbulence, by T. H. Ellison. 

The Mechanics of Sailing Ships and Yachts, by K. S. M. Davidson. 

The articles vary considerably in length from about twenty to over one 
hundred pages. Most of the surveys give a completely adequate account of the 
present state of their particular subject. The long article by Lighthill is, how- 
ever, rather more than a straight survey ; it contains some new and interesting 
material on the properties of shock waves. 

The editors and individual authors of this volume are to be congratulated on 
the production of a work which will undoubtedly help all researchers and 
teachers of fluid and solid mechanics to keep abreast with some of the subjects 
which have developed rapidly in recent years. The volume is also a most 
appropriate tribute to Sir Geoffrey’s contributions to research in ee fields 
during the last few decades. J. W.M. 
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Solid State Physics. Volume II. Edited by F. Surrz and D. TurnsuLv. 
(London: Academic Books Ltd.) [Pp. xii+468.] Price $10. 


Tux second volume of this ambitious series is as welcome and useful as the 
first. It contains articles on nuclear magnetic resonance by Pake and by Knight, 
on neutron diffraction by Shull and on specific heats by de Launay. Finally 
there is a massive and comprehensive article by Seitz and Koehler on the 
important subject of displacement of atoms during irradiation which in one 
hundred and fifty pages gives a most detailed survey of this subject. N.F.M. 


Valency and Molecular Structure. By E. Carrmetn and G. W. A. Fow es. 
(London: Butterworths Scientific Publications.) [Pp. xi+256.] Price 
32s. 6d. 

Tr1s book is intended as an introduction to valency theory and molecular 
structure for students beginning an honours course in chemistry. It begins with 
a clear and simple outline of the experimental foundations of the quantum 
theory and its application to isolated atoms, and then goes on to discuss the 
valence-bond and molecular orbital methods for the application of the theory to 
molecules. The discussion contains little mathematical treatment, and will 
undoubtedly be very useful to the many students who find difficulty in mastering 
the elements of theoretical chemistry. 

The second half of the book deals with the nature of the bonding in solid 
structures, the structures of simple inorganic substances, complex compounds, 
and electron-deficient molecules. In this section the scope of the work is much 
more ambitious, and the structures of the most simple and many complex 
inorganic compounds are described and interpreted in terms of electronic 
configurations, special attention being paid to substances of current interest, 
such as ferrocene, the higher hydrides of boron, interhalogen compounds, metal 
carbonyls, and metal alkyls. It is a pity, however, that the discussion is 
restricted throughout to molecular configuration and dimensions ; the energies 
and restoring force constants of covalent bonds provide important information 
which should take its place with the bond lengths to which the authors (and 
chemists in general) pay so much attention. Even when dealing with inter- 
atomic distances, the treatment is hardly critical enough; the structure of 
nitric oxide, for example, is discussed on the basis of ‘ calculated’ N=O and 


N—0 bond lengths with no indication that ‘ calculations’ of bond lengths in ~ 
heteronuclear inorganic molecules rest at the present time on intelligent but 


empirical correlations of observed bond distance and inferred bond order. 
These defects detract somewhat from the value of an otherwise useful introduc- 
tion to structural chemistry. 

The book is well produced, except for the printing in diagrams of structures ; 
this is sometimes almost invisible. AY Gimy 


Radiation Dosimetry. Edited by G. J. Hinn and G. L. Brownetu. (New York: 
Academic Press Inc.) [Pp. xiv-+932.] Price $22. 


Tr 1s book has eighteen chapters each on a different topic and each contributed 
by a recognised authority in the field. It is divided into three main sections ; 
the first entitled “‘ Fundamental Principles of Dosimetry” gives an excellent 
and detailed account of the physical principles involved in the interaction of 
radiation with matter and the application of these principles to dosimetry. 


The second section entitled ‘‘ Radiation Detectors and their Calibration ” — 


surveys the various means of quantitatively measuring radiation and includes 
an excellent and up to date chapter on “ Calorimetric Methods ’°—particularly 
pertinent in the light of the recent recommendations of the International 
Commission on Radiological Protection that radiation dose should be measured 
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as energy absorbed per gm of absorber. The third section—“ Radiation Fields 
and their Dosimetry ”—covers the whole field of dose measurement for electro- 
magnetic radiation with quantum energies from 200 kev to 30 Mev and in addi- 
tion there are chapters on high energy electrons, heavy charged particles, and 
neutrons used as external sources and another on the dosimetry of internally 
administered isotopes. 

The references, given at the end of each chapter, are excellent and up to date. 

This book gives an account of a rapidly expanding subject and collects in 
one volume material which previously could be consulted only in research papers 
and separate books. 

Tt can be highly recommended for all workers in the field. Cols: 


BOOK NOTICES 


Peaceful Uses of Atomic Energy. Proceedings of the International Conference 
in Geneva, August 1955. A United Nations Publication. Available from 
H.M. Stationery Office, P.O. Box 569, London, S.E.1. and other branches. 


Volume 2. “ Physics; Research Reactors.” [Pp. 471.] Price 57s. 

Volume 3. ‘“‘ Power Reactors.” [Pp. 389.] Price 54s. 

Volume 4. “Cross Sections important to Reactor Design.” [Pp. viii+357.] 
Price 54s. 

Volume 5. “ Physics of Reactor Design.”” [Pp. vili+-545.] Price 63s. 

Volume 7. “ Nuclear Chemistry and Effects of Irradiation.” [Pp. 691.] 
Price 70s. 

Volume 8. ‘“‘ Production technology of the materials used for nuclear 
energy.” [Pp. ix-+-627.] Price 70s. 

Volume 9. “* Reactor Technology and Chemical Processing.” [Pp.x+771.] 
Price 70s. 

Volume 11. “ Biological Effects of Radiation.” [Pp. 399.] Price 57s. 

Volume 14. “General Aspects of the Use of Radioactive Isotopes: 
Dosimetry.” [Pp. 305.] Price 45s. A 

Volume 15. “‘ Applications of Radioactive Isotopes and Fission Products 
in Research and Industry.” ([Pp. viii+327.] Price 54s. 

Volume 16. “ Record of the Conference.” [Pp. x+203.] Price 36s. 


Reactor Handbook. Produced by the United States Atomic Energy Commision. 
(London and New York: McGraw-Hill.) 
‘Research Reactors’. [Pp. 422.] Price $6.50. 
“Reactor Handbook: Physics”. [Pp. 709.] Price $12. 
“Reactor Handbook: Engineering”. [Pp. 1075.] Price $15. 
“ Reactor Handbook: Materials”’. [Pp.610.] Price $10.50. 
‘Neutron Cross Sections”. [Pp. 328.] Price $12. 
‘‘ Chemical Processing and Equipment”. [Pp. 302.] Price. $6. 
mati heory. Papers read at a symposium on Information Theory held 
ee te al JE ee oie London, September 1955. (London : Butterworths ; 
New York: Academic Press Inc.) [Pp. xii+401.] Price 70s. 
VUedical Effects of the Atomic Bomb in Japan. Edited by A. W. 
Celene eee 4 WarrREN. (London: McGraw-Hill.) [Pp. 477.] 
Price 60s. | 
Engineering Fluid Mechanics. By CHARLES JAEGER. (Glasgow: Blackie & Son 
Ltd.) [Pp. xviiit+-529.] Price 60s. 
The Fundamental Study of Mathematics. By Soraro NIrro. (Tokyo : 
Maruzen Co, Ltd.) [Pp. 83.] Price not stated, 
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The Diamond Tool Industry in 1955, Issued by the Industrial Diamond 
Information Bureau, 32-34 Holborn Viaduct, London, E.C.1. (pp. 12.) 


Free of charge. 


As in previous years, the Industrial Diamond Information Bureau has compiled 
a survey of the diamond tool industry, dealing with various scientific and 
technical aspects. The subject-matter, mainly based on the Bibliography of 
Industrial Diamond Applications, is sub-divided into 16 main groups ; about 
200 references are quoted. The literature mentioned can be made available to 
enquirers. 


Spheroidal Wave Functions. By J. A. Srrarron, P. M. Morsz, L. J. Cus, 
J.D.C. Lirrtz and F. J. Carparo. (London: Chapman & Hall.) [Pp. xi+ 
613.] Price 100s. 


Christian Theology and Natural Science. By E. L. Mascaty. (London: 
Longmans, Green and Co. Ltd.) [Pp. xxi+323.] Price 25s. 


The Theory of Photons and Electrons. By J. M. Jaucn and F. RouRticu. 
(London : Academic Books Ltd.) [Pp. ix+488.] Price $10. 


Mathematical Theory of Elasticity, second edition. By I. S. SOKOLNIKOFF.— 
(London: McGraw-Hill.) [Pp. xi+476.] Price 71s. 6d. 


Functional Analysis. By F. Rrssz and B. Sz.-Nacy. (Glasgow: Blackie 
and Son Limited.) [Pp. xii+468.] 65s. 


Microscopic Petrography. By E. W. HEtyRIcu. (London: McGraw-Hill.) 
[Pp. 296.] Price 49s. 


Kernmomente. By H. Kopfermann. (Frankfurt : Akademische Verlagsgesell- 
shaft.) [Pp. xvi+462.] Price DM54. 


Rectangular-Polar Conversion Tables. By E. H. Neville. (Cambridge 
University Press.) [Pp. xxix+109.] Price 30s. 
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ERRATA 


The Publishers regret that in the October number of the Philosophical 
Magazine a transposition of subject-titles occurred. 

On p 967 the title should read “Energy Levels of 404’? and not ‘‘ Periodicity 
in Ice-Flowers”’. 

On p. 975 the title should read “Periodicity in Ice-Flowers” and not 
“Energy Levels of *°A”’. 

The Publishers tender their apologies to the Authors concerned and trust: 
that no inconvenience has been caused to them. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


